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1. Introduction.—Since the appearance of Kapteyn’s exhaustive 
investigations of the extended moving clusters in the regions of Orion 
and Scorpio-Centaurus,' the motions of the B-type stars have been con- 
sidered primarily with regard to second order effects such as galactic rota- 
tion and the K-term. However, the recent publication of catalogues of 
radial velocities and spectroscopic parallaxes of B stars suggests the 
possibility of gaining additional information concerning the general state 
of the motions within these two clusters. In Harvard Bulletin 886 were 
published some preliminary results which seemed to call for a more ex- 
tensive study. Accordingly, the present paper summarizes an investi- 
gation of the space motions of 166 stars of classes BO—B9 in the region of 
Scorpio-Centaurus, and 57 in that of Orion. The limits of these areas in 
galactic coérdinates are those fixed by Kapteyn: for the former, / = 
216-360°, b = +30°; and for the latter, 7 = 150-216°, b = +30°. 

2. Computation of the Velocity Components.—Lick Observatory radial 
velocities,? Boss proper motions? and spectroscopic parallaxes by C. J. 
Anger‘ and W. B. Rimmer‘ have been utilized for the computation of the 
space velocities of the stars mentioned above. The corrections of Boss 
and Jenkins and of Oort’ have been applied to the proper motions, which 
are therefore as free as possible from systematic errors. The spectro- 
scopic parallaxes from the two sources have been reduced to a mean zero 
point, those contained in Dr. Anger’s second list being first reduced to the 
scale and zero point of her first paper. From this material, the galactic 
components of the space velocities have been computed in the usual 
manner. The +X axis is here defined as lying toward / = 0°, b = 0°; 
the + Y axis toward ] = 90°, b = 0°; and the +Z axis toward b = 90°. 
A check on the results was provided by the tables of A. Kéhlschutter.*® 
The components have been corrected for a solar motion of 20 km./sec 
toward R. A. 18" 20", Dec. +30°. Unless otherwise stated, all velocities 
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referred to in this paper are so corrected. In addition, it is necessary to 
consider the effects of galactic rotation, and K-term. The former may 
amount to 2 or 3 km./sec., but its influence, if any, will depend directly 
on distance from the galactic center, and should therefore be easy to detect. 
The effect of K-term, on the other hand, cannot be predicted. Correc- 
tions assuming K = +4 km./sec. were therefore computed, and could be 
omitted or applied at will. It may be supposed that the true velocity 
distribution will lie between the distributions found for the corrected and 
uncorrected velocities. 

3. Errors of the Observed Quantities—In a recent paper® Plaskett and 
Pearce find that the Lick radial velocities differ by only —0.4 km./sec. 
from their mean system. The systematic errors must therefore be small. 
They also find + 2.03 km./sec. as the mean probable error of one velocity. 
This and other published results indicate that the probable error may be 
assumed to be +2.0 km./sec. 

The systematic errors of u. and us; may be comparatively large. Since 
they must be of the same order of magnitude as the uncertainties in the 
corrections which were applied, it is not unlikely that they are as great 
as 07004, about 10% of the proper motion in either coérdinate. The 
mean probable error is found from Boss’s catalogue to be close to +0006 
in both coérdinates. 

It is difficult to make a preliminary estimate of the probable size of the 
systematic errors in the spectroscopic absolute magnitudes. It is there- 
fore permissible to present here some results which were derived later in 
the investigation. If the systematic errors in the proper motions are as 
large as indicated above, the stream parallaxes on which the adopted zero 
point is partly based may be in error by an amount equivalent to 0.2 
magnitude in the resulting absolute magnitudes. A comparison of the 
absolute magnitudes of Dr. Anger and Rasmuson,” the latter derived 
from stream parallaxes, shows systematic differences of 0.2 magnitude, 
depending on distance from the convergent of the cluster. Finally, the 
stream velocity of the Scorpio-Centaurus cluster as derived from the cross 
motions alone is but 84 per cent of that found from the radial velocities. 
The discrepancy may be partly accounted for by supposing that the 
assumed absolute magnitudes are, in the mean, too faint. These con- 
siderations, however, do not unite in giving evidence as to the direction 
in which the zero point is in error, but serve to indicate that this error 
may be as large as 0.2 or 0.3 magnitude. The probable error of an ahsolute 
magnitude is certainly no greater than +0.5 magnitude. 

4. Treatment of the Material.—It was desired to present the general 
characteristics of the velocity distributions in graphical form. The di- 
rection of motion and projected velocities in the XY, XZ and YZ galactic 
planes have been computed for each star. The direction in the X Y plane 
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corresponds to the galactic longitude of the direction of motion of the star, 
while the directions in the two remaining planes are defined by 0xz = tan~! 
Vz/Vx, and 6yz = tan~! Vz/Vy. Then for each plane, the sum of the 
projected velocities falling within successive ten degree intervals of the 
corresponding angle @ was computed, the results smoothed by taking 
running means of three ten degree intervals, and the smoothed means 
plotted against the corresponding values of @. This procedure, it will be 
seen, is similar to that employed by Kapteyn! in his first discussion of 
star streaming as exhibited by proper motions. 

5. The Theoretical Velocity Distribution—Before discussing the ob- 
served velocity distribution, it is necessary to consider in detail the effect 
of accidental errors. In the study of proper motions in a limited region 
of the sky, it may be expected that the accidental errors of uw, and uy; 
will be nearly the same. In dealing with space motions, the corresponding 
assumption cannot be made concerning the errors of the three space 
velocity components. The components are computed from the radial 
velocities and cross motions, and in a specific case, one space velocity 
component may be influenced primarily by the errors of the radial ve- 
locities, another by those of the proper motions and parallaxes. Since the 
radial velocities are observed in a manner fundamentally different from 
the cross motions, it cannot be assumed a priori that the errordi spersions 
along the X, Y and Z axes will be the same. It is apparent that the 
general problem is analytically similar to that of the Schwarzschild Ellip- 
soidal Hypothesis, which also supposes unequal velocity dispersions along 
the axes, and a superposed group motion relative to the sun. Considering 
only two dimensions, three cases merit consideration. The first is that 
in which the error dispersion is the same in both coérdinates, the second 
that in which the dispersion perpendicular to the direction of group motion 
is the larger, and in the third the axis of greatest dispersion is oriented 
arbitrarily with respect to the direction of motion of the cluster. The 
necessary formulae may be adapted from those given by Eddington.” 
For the first case, let NV be the total number of stars, V the group velocity, 
r the observed velocity of any star and / related to n, the observed mean 
peculiar motion in either coérdinate by h = 0.5641/y. Then the number 
of stars moving with velocities between 7 and r + dr in directions between 
6 and 6 + dé is given by 


Nh? 


T 


ee? + V2 — 2rV cos 6) rdrd6 


whence it follows that the sum of the velocities of all stars with velocities 
equal to or less than R, moving between 0 and 6 + d@, is 
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R 
Nh? a ae 
ll e her? + J 2rV cos 6 rdr: 
. 0 


It is necessary to set a finite upper limit.R, since the computed curve is 
for comparison with those found by observation, and in practice, the 
inclusion of stars of extremely high velocity produces irregularities out of 
proportion to the number of stars responsible. The function was inte- 
grated numerically with V = 5 km./sec., and a value of / chosen on the 
assumption that the mean of the errors in either coérdinate combined 
with a true mean peculiar motion of +2.0 km./sec. is +3.8 km./sec. 
The integration was performed with two values of the upper limit R. 
The resulting curves have a single symmetrical maximum, and are both 
strongly concentrated in the direction of motion of the cluster; the higher 
the upper limit, the more concentrated the curve. The development of the 
second case is similar to the preceding one. The number of stars with 
velocities between r and r + dr traveling between 6 and 6 + dé is given by 


ns 


T 


— r2(k2 eos? 0 2 sin? LV cos 6 — k2V2 
e r?(k2 cos + h? sin? 0) + 2rk?V cos 6 kV? +drd6 


which leads at once to the following expression for the sum of all velocities 
less than or equal to R lying between 6 and 6 + dé: 


R 
Nhk eV" wf e~ r? (k? cos? @ + h? sin? 0) + 2rk?V cos 6 rdr- 
0 


T 


Here V, 7 and R are defined as before, and k and / are related as described 
above to the observed mean peculiar motions along the X and Y axes. 
Since V is along the X axis, k will be larger than 4. For the integration, 
V was again taken as 5 km./sec., and the observed mean peculiar motions 
along the X and Y axes were supposed to be +2.5 and +5.0 km./sec., 
respectively. The resulting distribution makes apparent the need for 
care in the interpretation of the observations. The function has a pro- 
nounced minimum in the direction of group motion, with maxima distant 
about 40° on both sides. This is in general agreement with the appear- 
ance of the observed distribution discussed in Harvard Bulletin 886, and 
demonstrates that a non-spherical error dispersion is capable of simulating 
the existence of two intermingled clusters or streams. The appearance 
of the velocity distribution in the third case may be predicted qualita- 
tively from its analogue in the Ellipsoidal Hypothesis. If the axis of 
greatest dispersion is inclined at some angle not 90° to the direction of 
motion of the cluster, there will be.two maxima of unequal size, the more 
pronounced being the nearer to the true direction of motion of the cluster. 

6. Comparison with Observation.—The distribution of the velocities of 
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the fifty-seven stars of the Orion region is illustrated in figures la and 10. 
The heavily drawn curves represent velocities corrected for K-term, the 
lighter curves are for velocities not so corrected. In addition, in figure 1), 
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FIGURE 1 


centered on an arbitrary zero point, is included the theoretical distribution 
of the fifty-seven stars, computed on the assumption of a spherical error 
dispersion. Inspection indicates that while the inclusion of the K-term 
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correction effects numerous minor changes in the appearance of the curves, 
the general contours are unaltered. It will also be noted that the dis- 
tributions are somewhat asymmetric, probably due to the fact that the 
error dispersion is not truly spherical. A comparison with the theoretical 
distribution shows that although the forms are not strictly comparable, 
the observed dispersion is no greater than would be anticipated if prac- 
tically all the stars with velocities less than 21 km./sec. belonged to a 
single cluster. The number of velocities greater than 30 km./sec. is six, 
rather more than would be expected, but in all cases except one it may 
be shown that the high velocity is probably due to a large error in the 
parallax, or that the star is a very distant one lying far beyond the cluster. 
It may therefore be concluded that the great bulk of the class B stars 
falling within the Orion region form a single homogeneous cluster. 

In figures 1c to 1f are presented the principal characteristics of the 
velocity distribution of the stars in the Scorpio-Centaurus region. In 
figure lc the dotted outline is a reproduction of the velocity distribution 
of the original ninety-five stars treated in Harvard Bulletin 886. It will 
be seen that its general appearance is well reproduced by the more com- 
plete material. Considering first the X Y velocities corrected for K-term, 
the outstanding feature is the double maximum representing the motions 
of the greater part of the stars. In figure ld this curve is repeated, and, 
superposed on it, the theoretical distribution computed as in the second 
case of Section 5, on the supposition that the error dispersion is not spheri- 
cal. The constants V, 4 and & in the function were chosen purely ac- 
cording to their most probable values. The only deliberate fitting con- 
sisted in assuming that there are only 105 stars in the cluster. The true 
number is larger, but in view of the number of disposable constants in 
such a function, it would have been entirely feasible to fit a curve with the 
right number of stars to the observed distribution. The purpose of the 
computation is served equally well, however, by this demonstration that 
the distribution of purely accidental errors can so well imitate the appear- 
ance of two intermingled clusters or streams. The distributions in the 
XZ and YZ planes are not reproduced here, since they contribute little 
information concerning the state of motion within this region. They 
show, however, a strongly systematic tendency toward positive Z com- 
ponents. Since this is shared by the Orion stars, it may be the result of 
systematic errors, or an error in the assumed value of the Z component of 
the solar motion. In addition to the principal maxima, there appear 
two secondary peaks which will be considered in a later section. How- 
ever, it is clear that by far the greater part of the stars in this region form 
a single cluster, and that the appearance of a double stream is entirely 
due to an ellipsoidal error dispersion. 

7. The Convergents—The convergent of the Orion cluster has been 
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computed by taking the mean values of the X, Y and Z velocity com- 
ponents, respectively. Since all stars in the Scorpio-Centaurus region 
certainly do not belong to a single cluster, it seems that better results 
may be obtained simply from an inspection of the positions and forms of 
the principal maxima of the distribution curves. The elements given 
below are the means of the determinations with and without the K-term 
correction, and have been reduced to the sun for comparison with the 
values derived by Kapteyn. As inspection of the curves does not give the 
velocity, a value of 5 km./sec. was assumed for the Scorpio-Centaurus 
cluster for making the reduction to the sun. 


ORION CLUSTER SCORPIO-CENTAURUS CLUSTER 
EAPTEYN SPACE MOTIONS KAPTEYN SPACE MOTIONS 
R.A. 5° 44” eS 6° 18” T ? 
Dec. —11° —13° —42° —38° 
V 20.0 km./sec. 22.7 km./sec. 18.3 km./sec. 


Part of the discrepancy between the velocities in the Orion cluster can 
be traced to the difference between the K-term correction applied by 
Kapteyn, and that assumed here. 

8. Secondary Maxima.—Two secondary maxima appear in the heavily 
drawn curve of figure lc. In the uncorrected distribution, both have 
disappeared, but there is still evidence of excess motion in the same 
directions. The more prominent peak at longitude 95° is of particular 
interest, since Kapteyn’s elements for the Orion cluster, after being freed 
of solar motion, give a longitude of 108°, and Rasmuson finds 80°. Al- 
though the convergent indicated by the curves of figure la is nearer 130°, 
it must be remembered that the cluster velocities are on the order of 5 
km./sec., and that the direction of the velocity vector is therefore at the 
mercy of systematic errors. Dr. Anger’s space diagrams of the Scorpio- 
Centaurus and Orion clusters show that their outlying members appear 
to intermingle, and suggest that the Orion cluster is not highly con- 
centrated, whence it would not be surprising to find a few of its scattered 
members in the region of space occupied by the Scorpio-Centaurus group. 
However, the physical reality of the secondary maximum must first be 
tested. It was noted in Section 5, that the concentration of a velocity 
distribution with a single maximum should be increased with increase of 
the upper limit of velocities included. This is, of course, only true if there 
be a physical moving group. If the maximum be purely accidental, it 
should be obliterated by the inclusion of stars of higher and higher ve- 
locity. Figures leand If illustrate the result of such anexamination. The 
maximum at longitude 95° increases markedly in concentration, a result 
which should not appear were the original peak purely fortuitous. The 
mean velocity of the stars contributing to the peak is found to be in good 
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agreement with the mean velocity of the stars in the Orion region. Thus, 
on the grounds of the close proximity of the Orion and Scorpio-Centaurus 
clusters, and the agreement in velocity and direction of the motions of the 
stars in this subgroup and in the Orion cluster, it may be concluded that 
far flung members of the latter group have been detected in this manner. 

No identification can be offered of the secondary maximum which 
appears at longitude 195°, although there is evidence that it is not acci- 
dental in character. 
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FIGURE 2 


9. The True Velocity Dispersion—Kapteyn considered that the most 
important results of his discussion of these two clusters were the small 
values (+1.67 and +1.00 km./sec.) of the true mean peculiar motions. 
It is of interest to discover the extent to which the new material checks 
these figures. The true dispersion is best determined separately from 
the radial velocities and from the r components reduced by means of the 
parallaxes to km./sec. It was found that after correcting for the probable 
errors of the radial velocities and proper motions given in Section 3 (the 
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parallax errors play no part here), the true mean peculiar motion in one 
coérdinate as given by the radial velocities was from 1.5 to 5 times that 
derived from 7 components, no matter what method of selection was 
utilized to eliminate stars not members of the clusters. The only solution 
lies in assuming smaller errors in the proper motions, or larger ones in the 
radial velocities. Fortunately, an independent determination of the 
accidental errors of the proper motions is possible. The stars of either 
cluster may be divided into two or more groups according to distance from 
the sun, and for each group, the mean distance computed, and the observed 
mean peculiar motion determined from the + components. Since the 
error in km./sec. produced by a given error in the proper motions will be 
directly proportional to the distance of the star, while the true mean 
peculiar motion may be supposed constant over the cluster, one has two 
or more equations of the form 


4.747 p, |? _ - 
=-2 _ Bb 2 
Mo E tf 


where % is the observed mean peculiar motion in one codrdinate, 7, the 
true mean peculiar motion, r the mean distance and p, the probable error 
of a proper motion component. These equations allow of a solution for 
the two unknowns, 7, and p,. Dividing the Orion cluster into two groups, 
and the Scorpio-Centaurus into three, the true mean motions were found 
to be +1.6 and +0.8 km./sec., respectively, in close agreement with the 
determinations of Kapteyn, and the mean of the two solutions for p, was 
found to be +0"006, the value given in Section 3. This must indicate 
the existence of inherent accidental errors in the measured radial velocities 
which are not brought out by intercomparison of several measurements or 
catalogues. The true mean peculiar motions within the clusters are 
certainly less than +2.0 km./sec. 

10. Star Streaming.—There is no evidence of the two star streams in 
these space motions. However, Wilson and Raymond" in a paper on 
“The Space Motions of 4233 Stars’’ discuss streaming from the viewpoint 
of each spectral class, and for 435 stars of class B determine the axes of 
the velocity ellipsoid. The division of the stars into the two streams is 
illustrated by a diagram in which the numbers of stars moving within 
successive thirty degree intervals of galactic longitude are plotted against 
longitude. It now seems probable that the two resulting maxima at 
longitudes 120° and 330° may be caused by stars of the Orion and Scorpio- 
Centaurus clusters, rather than by star streaming in the ordinary sense. 
Accordingly, the corresponding counts were made for the stars of the 
present paper, and the results embodied in figure 2. Here the upper 
curve is a reproduction: of the diagram of Wilson and Raymond, and the 
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lower curves are derived from the present material. The general agree- 
ment between the two sets of curves is so close that there is good reason 
to suppose that the apparent existence of the two streams is in fact the 
outcome of the presence of the two large moving clusters. 

The writer wishes to express to Dr. Bart J. Bok his appreciation of the 
continued advice and encouragement which have made this paper possible. 
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THE PENETRATION OF IRON METEORITES INTO THE GROUND 
By WILLARD J. FISHER 
HARVARD COLLEGE OBSERVATORY 


Communicated February 7, 1933 


Treatises on exterior ballistics show that objects moving under gravity 
in a resisting medium of uniform density ultimately descend with a velocity 
such that the resistance to motion balances the weight: an idea said to 
be traceable to Huygens. Schiaparelli extended the idea to meteors, and 
showed that, neglecting the attraction of the earth, the ultimate velocity 
approached is little affected by the initial or extra-atmospheric velocity, 
and little affected by the law of variation of air density with height [pro- 
vided only that the mass of the air immediately disturbed is large in com- 
parison to the mass of the object].! For large meteors the ultimate 
velocity does depend on the cosmic velocity. His analysis also neglected 
the effect of ablation by melting in the violent blast of hot air impinging 
on the front of the meteor. ; 

With meteorites of ordinary size the striking velocity at the ground, 
being independent of the extra-terrestrial velocity, should depend only 
on the mass and dimensions of the object and on the density of the air at 
the ground. For use of this principle it is necessary to take recourse to 
a formula for spheres? 
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v? = P/0.014 R? (1) 


in which P is the weight, R the radius and the units are meter-kilogram- 
second, engineering system. Air resistance is supposed proportional to 
v? and to air density. The known inaccuracy of this law of resistance in 
the neighborhood of the velocity of sound is regarded as unimportant for 
this application. 

In 1839 and 1840 a committee of the French artillery, consisting of 
Piobert, Morin and Didion, carried on a series of experiments at Metz 
and vicinity, part of which were devoted to determining the law of pene- 
tration of round shot into various materials of fortification. Beside the 
official, extended reports of this committee, their theoretical and practical 
results on penetration are presented in Didion’s Traité de Balistique, 
Paris, 1848, Section VI, and have been copied by Crantz, from whom they 
are quoted in part by Borgstrém.* The differential equation for the re- 
sistance of targets is based on a study of the forms and dimensions of 
penetration funnels and of projectiles dug out of the targets. The ex- 
pression for depth of penetration of cast iron round shot is 


E = K-2R:-log[1 + (v/u)?] (2) 


in which E is the penetration, R is the radius of the projectile and u and 
K are empirical constants, the density of the projectile being contained 
in K. The units are meter-kilogram-second, engineering. Tables of 
these constants are given for various targets. 

A slight modification of the constants enables the centimeter-gram- 
second absolute system to be used, and permits a substitution of m/* 
instead of radius of the projectile. 

In the case of a meteorite of modest size, for which a terminal velocity 
is calculable, the terminal velocity in the air path is the same as the striking 
velocity of the penetration into the ground; hence the velocity may be 
eliminated between two equations, and there results 


1 bg 
a 4% is 
sm Ook log E + abd m (3) 


in which s is the penetration into the target, m is the mass of the meteorite, 
g is the acceleration of gravitation, d is the density of air, k is a form factor 
which converts m’* into an effective resisting cross-section. & contains 
the density of the meteorite, b and a are empirical constants, functions of 
“the ground.” The units supposed are c.g.s. absolute. The formula 
means that the penetration of a meteorite of given shape and material 
into a given kind of ground is proportional to the cube root of the mass of 
the meteorite multiplied by a factor which increases with the cube root 
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of the mass, but less rapidly. But there is no assumed dependence of 
penetration upon cosmic velocity. 

As the expression consists of products and powers, the obvious way to 
test the rough data given by the records of meteorite falls is by plotting 
on logarithmic paper. If penetration were strictly proportional to cube 

TABLE 1 


ACTUAL PENETRATIONS OF IRON METEORITES 


REPORTED REDUCED 
METEORITE MASS DEPTH KG. CM. NOTES 
Avée Hex 1.28 kg. 30+ cm. 1.2 30 Broke apple branch 5 


cm. thick. 
Boguslavka Hex 199 kg. 130 cm. 199 130 Sandy ground. 
Boguslavka Hex 57 kg. 200 cm. 57 200 Clay. 
Braunau Hex 23.625 kg. 100 cm. 23.6 100 Vertical hole. 
Cabin Creek MO 107 Ib $ ft. 48.5 91 Broke limb of a pine. 
Estherville MS 437 Ib. 8 ft.; 6 ft. 198.0 244; 182 Passed first through 6 
ft. water into clay. 


Estherville MS 170 Ib. 4.5 ft. 77.0 Dry gravelly soil. 
Estherville MS 92.51b. 5 ft. 41.9 5 

Hassi Jekna FO 11/4 kg. 80 cm. 1.25 80 Desert sand, by Arabs. 
Hraschina MO 71 Ib. 3 Klafter 39.8 569 Recently plowed field. 
Mazapil MO 3.95 kg. 30cm. 3.95 30 Corral. 

Nedagolla NPA 10+1b. !/2m. 4.5 50 

N’Goureyma BrO 37'/2 kg. 1m. 37.5 100 Clay. 

Pitts O 421/.0z. 8 in. 1.2 20 

Pitts O 57 oz. 16 in. 1.6 41 Fresh plowed land. 
Pitts O 30 oz. 7 in. 0.85 18 

Quesa FO 10.67 kg. 40cm. 10.7 40 Obliquely. 

Rowton MO 73/4 Ib. 18 in. 3.5 46 Turf field. 

Treysa MO 63 kg. 160 cm. 63 160 In a wood. 


Hex = hexadrite; MO = medium octahedrite; MS = mesosiderite; FO = fine 
octahedrite; NPA = nickel-poor ataxite; BO = brecciated octahedrite; O = octa- 
hedrite. (No nickel-rich ataxites have been known to fall.) 

Where there are no notes on soil, there is no information. 

In some cases, the masses were diminished artificially before weighing. The measure- 
ments of penetrations are usually rough, and sometimes they are estimates. 

The penetrations of Avée, Cabin Creek and Estherville (437 Ib.) must be affected 
slightly by the obstacles met above ground. 

The penetration of Estherville (437 lb.) is given diversely by C. U. Shepard and J. 
Lawrence Smith. 

Falls lacking penetration data are omitted. In plotting, masses less than 1 kg. are 


omitted. 


root of mass, the points plotted should lie on one or more straight lines, 
slope = '/3, the intercept of each line depending on the kind of ground 
hit. The variability in form and density of meteorites, the differences 
in ground encountered, and the usual badness of the measurements of 
penetration, with occasional badness of the measurements of mass, produce 
a scatter of points plotted sufficient to prevent complete certainty. 
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The data collected from the literature of meteorite falls are naturally 
grouped by density into two classes, for irons (and stony-irons) and for 
stones. In the present paper only the data for irons and stony-irons are 
tabulated. The sources are those indicated in Prior’s Catalogue and his 
Appendix, but Farrington’s Meteorites of North America, Flight’s Chapter 
in the History of Meteorites, and Chladni’s Feuermeteore, and certain 
encyclopedias, have been used where possible. 

The data for these meteorites are shown in table 1, with notes on the 
ground penetrated—unfortunately few—and other points of interest. 
The figure, somewhat reduced, shows the same data plotted on logarithmic 
paper. Abscissas are kilograms; ordinates, centimeters. The dash line 
has the slope 1/3, for comparison. The plotted points lie in a band be- 
tween two lines, the lower of which is nearly parallel to the comparison 
line. (The two highest points are plotted for completeness; they have 
probably little relation to facts. The measurements for Hassi Jekna are 
probably unreliable, and the penetration of “drei Klafter,” or three fath- 
oms, of the Hraschina iron, is incredible.) 

For the purpose of extending the discussion to the great iron meteorites, 
the lower bounding line will be taken to represent the lower limit of the 
penetration of iron meteorites into the ground, and its equation will be 
taken to be s = 18.6 m”. 

“The ground” may be represented by clay, in Didion’s tables, for which, 
and for cast iron round shot, we take ‘“Terre argileuse de Saint-Julien,”’ 
K = 19.9, uw = 112 and 


penetration E; = 19.9-2R-logi[1 + (v/m)?]. (4) 
For limestone, ‘‘Roches de calcaire oolitique des Genivaux,’’ K = 3.05, 
u = 258 and 
E, = 3.05-2R-logi[1 + (v/ue)?]. (5) 
With these two formulas, and with a value of striking velocity computed 


for terminal fall, the relative penetrations of cast iron round shot of equal 
masses may be roughly computed for falls into clay and limestone by 


E, _ 19.9 logi{1 + (v/112)*] 


E, 3.05 logi[1 + (v/258)?] 





(6) 


M.K.S. units. Equation (1) is used for computing v. ; is computed 
with E, = s = 186m’. 

Table 2 shows the masses, minimum penetrations into clay, terminal 
velocities and minimum penetrations into odlitic limestone, of certain 
large iron meteorites, treated as cast iron round shot. 

The scatter of the points in the figure does not admit a determination 
of the effect of the bracketed factor in equation (3). Table 2 shows that 
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TABLE 2 


PENETRATIONS OF GREAT IRON METEORITES (COMPUTED FOR CLAY AND LIMESTONE) 


TERMINAL 
WEIGHT, CLAY, VELOCITY, LIMESTONE, 


METEORITE KG. CM. cm./SEC. CM. REMARKS 
Bacubirito 27000 = 558 807 51 Lies in a hollow in rock; 
its crater? 
Bendego 5000 = 318 608 27 
Bethany: 
Two in Windhoek ae on a * Three of Bethany se- 
Mukesop 178 105 349 - lected from over 50. 
Cape York (1) 36500 + 617 848 57 a Cope bee >on 
Cape York (2) 3000 + 268 560 22 ree iar oehig cent 
Cape York (3) 500 + 148 415 if, ee 
Cape York, Savik 3500+ 282 574 23 j San NE Eaten er 
fell on snowbanks. 

Chupaderos (1) 14114 449 723 40 

Chupaderos (2) 6767 352 641 30 

Cranbourne 3500 282 574 23 

El Morito 11000 + 414 695 36 

Hoba West 88000 + 827 982 80 Weight before corrosion 
(after L. J. Spencer). 

Willamette 13500 = 443 719 39 Base up, and level with 
the ground. 


the minimum (computed) penetration into limestone is such that even 
with a considerable increase in penetration the meteorites would still be 
easy to find without digging. Even the greatest of them all, Hoba West, 
has actually penetrated the Kalahari superficial limestone to only its 
own thickness, and lies with its flat base level with the ground. But the 
minimum (computed) penetrations into clay show that such large irons 
would produce deep craters in ‘‘the ground’’; not so large, however, but 
that in the course of centuries the holes would be effaced by the wash of 
the weather. 

Herein may lie an explanation of the well-known fact, that the large 
iron meteorites are so generally found in rocky uplands. There only 
would they be stopped within easy reach of the prospector. Buried a 
meter or more deep, ‘‘the ground’’ of prairies and valleys may contain 
many great iron meteorites, or the rust of their corrosion. Had the fifty- 
one Bethany irons fallen into the Nile delta, there they would still be, and 
there they would stay. 

Similar plotting of the numerous data for stony meteorites brings out 
little of interest. The stones are specifically so much lighter than the 
irons, and are usually so small, that the top soil has an inordinate effect on 
their penetration, it being wet or dry, plowed land or grassy. Sometimes 
this results in reversal of the usual rule of penetration; the smaller of the 
Ellemeet or St. Michel or Parnallee pairs, etc., penetrated farther than the 
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larger. Furthermore, the penetrations have often been evidently not 
measured, but estimated in easy multiples of inches, feet or decimeters. 
Even the weights are frequently merely easy multiples of the pound and 
the kilogram. 

A plot on logarithmic paper shows a wedge-shaped cloud of points, 
narrowing upward toward large weights, with an average slope of 1/3 or 
more—certainly not less. The effect of moisture in the deep subsoil is 
probably well shown by the three Tilden stones, all of which fell on parched 
clay. The larger the stone, the deeper the penetration, but in ratios not 
explainable by mass alone; apparently the going was easier in the deeper, 
moister clay. Lines connecting the Stalldalen points, or the Drake Creek 


cm 


—_— 





10 


Logarithmic plot of the masses (kilograms) and earth penetrations 
(centimeters) of iron meteorites. 


points, or two of the New Concord points, show a slope not much different 
from '/;. At present the data on the penetration of stony meteorites into 
the ground seem unpromising, even for the derivation of qualitative results. 


1 Ahnighito and Hoba West are the largest meteorites known; but neither is large 
in comparison with the air mass which it would immediately disturb, even in vertical 
fall. So that the word large is used in two senses, both implied but not definitely dis- 
tinguished by Schiaparelli. The proviso in brackets is introduced by the writer. 

2 Quoted from C. Crantz, Compendium der theoretischen dusseren Ballistik, Leipzig, 
1896, p. 64, by L. H. Borgstrém, Die Meteoriten von Hvittis und Marjalahti, Helsingfors, 
1903, p. 7. The book is not accessible in Cambridge, and the constant 0.014 has not 
been verified. 

3 Crantz, Compendium, p. 362; Borgstrém, Hvittis und Marjalahti, pp. 6, 10. 
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Communicated January 26, 1933 


In the course of a comparative investigation of the elements of the 
third group in the Brown laboratories, we have had occasion to prepare 
trimethyl gallium. In our initial attempt, we employed the familiar 
Grignard reagent, but obtained a compound which proved to be an ether- 
ate, from which trimethyl gallium could not be prepared in the pure state. 
We were, however, able to substitute ammonia for ether and thus obtain 
an ammoniate. In the meantime, Dennis and Patnode* prepared tri- 
ethyl gallium by the action of diethyl mercury on metallic gallium at 
elevated temperatures (160°). By means of the Grignard reagent they 
obtained only an etherate. Very recently Renwanz* describes the prepa- 
ration of trimethyl gallium etherate by means of the Grignard reagent. 
Below, we communicate our results on the preparation and properties of 
trimethyl gallium as well as those of trimethyl gallium etherate and tri- 
methyl gallium ammine. In a second paper, we will communicate our 
results on chlorine derivatives of trimethyl gallium. 

Trimethyl Gallium.—(1) Preparation. We prepared trimethy] gallium 
by the action of dimethyl zinc on gallium trichloride. Approximately 
20 g. of gallium trichloride were treated with dimethyl zinc in an appara- 
tus of the form outlined in figure 1. The dimethyl zinc was prepared 
by the action of methyl iodide upon the familiar zinc copper couple; 
the reaction mixture was heated to boiling for 40 hours, when the methyl 
iodide had disappeared and the contents of the reaction vessel were com- 
pletely solidified. On heating the reaction mixture, dimethyl zinc began 
to be formed at a temperature of 80° and the reaction was completed at 
150°. The dimethyl zinc was condensed and collected under an atmos- 
phere of nitrogen in an ampule cooled in liquid ammonia. The product 
was almost completely free from methyl iodide, a trace of the latter being 
removed by the addition of a very small quantity of sodium to the di- 
methyl zinc at a temperature of 40°.4 The product was redistilled and 
condensed into ampules of known capacity for use as required. These 
ampules, of the form shown at A in figure 1, were provided with a small 
capillary B. When it was desired to withdraw dimethyl zinc, the ampule 
was attached to the apparatus as shown in figure 1 by means of deKhotin- 
sky cement, the capillary projecting through the stopcock C. The ap- 
paratus being exhausted, on turning the stopcock the capillary was broken 
and dimethyl zinc vapor could be introduced into the apparatus as required. 
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Gallium trichloride, prepared by the action of chlorine on metallic 
gallium, was sublimed several times and introduced directly into the re- 
action tube D through tube E, which was then sealed off. A second tube 
F was joined to D and this tube was in turn joined to the ampule 4A, 
containing dimethyl zinc as already described. The apparatus was attached 
to a manometer at G and to a pump and a source of dry nitrogen through 
a trap H, which was cooled in a bath of liquid ammonia. 

Dimethyl zinc vapor reacts vigorously, and even violently, with gal- 
lium trichloride at room temperatures. The reaction is controlled by 
means of the stopcock C which limits the rate of flow of vapor. As the 
reaction proceeds, it becomes less vigorous and, when a third or more of 
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FIGURE 1 


the dimethyl zinc has been introduced into the reaction tube, it is neces- 
sary to heat the reaction mixture to temperatures ranging from 80° to 
120°. When approximately one half of the dimethyl zinc has been intro- 
duced, trimethyl gallium begins to be formed. As the tube D is heated, 
tube F is cooled in liquid ammonia and the gallane, together with unreacted 
dimethyl zinc, is condensed. The reaction tube is cooled in ice water. 
dimethyl zinc is again introduced through stopcock C and the liquid in 
F is recondensed in the reaction tube. The reaction mixture is then 
heated to approximately 120°. The course of the reaction can be followed 
by means of the manometer. The volatile products in the reaction tube 
are again condensed in F, dimethyl zinc is condensed on the material in D 
as before and this process is repeated until the dimethyl zinc has been 








294 CHEMISTRY: KRAUS AND TOONDER Proc. N. A. S. 


used up. The volatile materials are heated and recondensed in D several 
times to insure complete disappearance of the dimethyl zinc. 

The total reaction takes place according to the equation: 

2GaCl; + 3(CHs3)eZn = 2Ga(CHs)3 + 3ZnCl, 

and is practically quantitative. In the initial stages of the reaction, the 
gallium is only partially methylated. Dimethyl zinc, however, reacts 
completely with the formation of zinc chloride. As the trimethyl gallium 
is formed in the course of the reaction, its presence is evidenced on con- 
densing the product in F at liquid ammonia temperatures. Dimethyl 
zinc freezes at —40° while trimethyl gallium freezes at —19°. At liquid 
ammonia temperatures, the gallane separates out in crystalline form in 
the presence of considerable quantities of dimethyl zinc. 

When the reaction is completed, the product is condensed in F and the 
tube is sealed off from the remainder of the apparatus. The product 
may be transferred from the tube through the attached capillary as al- 
ready described. The crude product was distilled several times and sam- 
ples were condensed into weighed, fragile, glass bulbs for analysis. The 
remainder of the product was filled into a small ampule provided with a 
capillary for use as needed. 

(2) Amalysis.—For analysis, trimethyl gallium was oxidized to GazO3 
in weighed Pyrex test tubes. The compound, sealed in a weighed, fragile, 
glass bulb, is introduced into a test tube which is then cooled to a tempera- 
ture of —76° in freezing ammonia and the fragile bulb is broken. The 
compound is so reactive toward air that it ignites and oxidizes with evolu- 
tion of light in a bath of freezing ammonia. No product is lost in this 
process and, when the initial process of oxidation is over, sulphuric and 
nitric acids are added to complete oxidation. The acids are fumed off 
and the tube is weighed in the usual way. 

Analysis: substance, 0.2554, 0.1858; GayO3, 0.2090, 0.1519; per cent 
Ga found, 60.87, 60.81, mean 60.84; Ga required for (CH;)3Ga, 60.72. 

The molecular weight of trimethyl gallium vapor was determined by 
introducing the vapor into a bulb of 347.4-cc. capacity at known pressures 
(approximately 15 cm.). The result of several determinations gave a 
mean value of 115.1 against the theoretical value of 114.8. 

(3) Properties —Trimethy] gallium is a volatile liquid boiling at 55.7 + 
0.2° under a pressure of 762 mm. and freezing at —19°. Its vapor pres- 
sure is given by the equation logioPmm, = —1705/T + 8.07. At —33.5°, 
0° and 25°, the vapor pressure is 8.9, 66.6 and 221.8 mm., respectively. 
The molal heat of vaporization according to the vapor pressure equation 
is 7800 calories and the specific heat of vaporization is 68 calories. Tri- 
methyl gallium is extremely reactive toward oxygen. The solid reacts 
pyrophorically with air at —76°; a trace of the vapor combines explosively 
with air at ordinary temperatures. 
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Trimethyl gallium combines with ether to form an etherate which boils 
at 98.3°. It forms an ammoniate melting at 31°, which is fairly volatile 
at room temperatures. It reacts with hydrogen chloride in the pure state 
or in solution in ether, one or more methyl groups being substituted by 
chlorine. 

Trimethyl Gallium Monoetherate.—(1) Preparation. Approximately 30 g. 
of gallium trichloride were dissolved in 150 cc. of ether and added slowly 
to twice the equivalent amount of methyl magnesium iodide in ether solu- 
tion under an atmosphere of nitrogen. The reaction was not vigorous, 
although there was a noticeable evolution of heat. The apparatus em- 
ployed is shown in figure 2. The gallium trichloride, after preparation, 
was sublimed into the tube A, which was then detached and connected to 


B f° 











O 


H 














FIGURE 2 


a tube containing ether (not shown in the figure), and 150 cc. of ether were 
condensed on the trichloride. The siphon tube B was then introduced 
through the top of the condensation column at C and the solution of gal- 
lium trichloride was blown into the Grignard reagent contained in the 
flask D. The siphon tube was then removed and replaced by a ther- 
mometer. The apparatus was connected with a source of nitrogen at 
E, a pump at F and a manometer at G. The reaction mixture was re- 
fluxed for 40 hours, when about two-thirds of the ether was distilled off. 
At this point, a white solid began to appear. After again refluxing the 
mixture for 12 hours, most of the remaining ether was distilled on a bath 
at 170°. The residue in the reaction flask was then heated to the boiling 
point for an additional’4 hours. The ether which had collected in H was 
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distilled into J under reduced pressure and flask J was detached and re- 
placed. The system was then evacuated to 3 cm. and, with an ice bath 
surrounding H and a liquid ammonia bath surrounding K, the reaction 
flask D was gradually heated to 260°. The liquid distilled into the re- 
ceivers H and K.. The crude product in H was distilled under a pressure 
of 5 cm. into receiver K, cooled in liquid ammonia. The product was 
finally purified by distillation through a fractionating column. The liquid 
was collected in ampules which were provided with capillaries as already 
described. From these ampules, the liquid was transferred by distilla- 
tion as needed. Based on the gallium trichloride used, a yield of 90% 
was obtained. 

(2) Analysis —The compound was analyzed for gallium by oxidizing 
known samples with fuming sulphuric and nitric acids in a Pyrex test 
tube and weighing as Ga,O;. The samples were weighed in sealed fragile 
bulbs similar to those described in the preceding section. The bulb was 
broken in a test tube cooled in liquid ammonia and 2 cc. of fuming sul- 
phuric and 5 cc. of nitric acid were added dropwise. When the reaction 
at liquid ammonia temperature had subsided, the tube was cautiously 
warmed to room temperatures, being re-cooled if necessary. After fum- 
ing off the acids, the tube was weighed. The results are tabulated below. 


TABLE 1 


ANALYSIS OF TRIMETHYL GALLIUM MONOETHERATE 
SAMPLE GazOs % Ga 
0.4345 0.2153 36.86 
0.2994 0.1502 37.32 
0.2931 0.1500 38.07 
0.2972 0.1499 37.52 
0.4180 0.2085 37.10 


The gallium content of the various samples varied from 36.86 to 38.07; 
the gallium content of the monoetherate, (CH3)sGa-(C2Hs)20, is 36.91%. 
The variation in the gallium content of the samples analyzed is due to 
the fact that the etherate is largely dissociated in the vapor phase and 
measurable separation of the two components occurs during distillation. 

(3) Properties —Trimethyl gallium monoetherate is a colorless liquid 
which fails to freeze at —76° and boils at 98.3°. It is readily soluble in 
liquid ammonia. The vapor pressure is given approximately by the equa- 
tion: logiwPmm. = —2440/T + 9.45. The vapor pressure is 0.18, 3.2 
and 18.3 mm. at —33.5°, 0° and 25°, respectively. From the* vapor 
pressure curve, the molal heat of vaporization is calculated to be 11,200 
calories. 

The vapor density of trimethyl gallium monoetherate at room tempera- 
tures was determined by weighing a known volume (of vapor) at a pres- 
sure somewhat below the saturation pressure. The density varied some- 
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what with the sample employed, the molecular weight ranging from 107 
to 115. The calculated molecular weight for the monoetherate is 189, 
whence it follows that the vapor is dissociated to the extent of about 70% 
at 25°. 

Trimethyl gallium etherate is oxidized slowly by atmospheric oxygen 
if the exposed surface is not too large. It is hydrolyzed by an aqueous 
solution of potassium hydroxide. At room temperatures only one methyl 
group is substituted while, at 100°, two methyl groups are replaced by 
hydroxyl. In two experiments, at 35°, the methane was collected and 
the gallium residue was analyzed for gallium. In two other experiments 
at 100°, the methane was collected but the gallium was not determined. 
In this case, however, the gallium content was known fairly closely from 
the amount of etherate used. The results are as follows: 

(4) Determination of Methane, 35°, m. mols. compound, 1.145, 1.224; 
m. mols. CH, 1.16, 1.31: 100°, m. mols. compound, 1.745, 1.28; m. mols. 
CHy, 3.35, 2.55: ratio CH,/Ga; (35°) 1.01, 1.07; (100°) 1.92, 1.99. 

The salts formed upon hydrolysis are soluble in water but have not been 
studied further. 

Trimethyl Gallium Monoammine.—When liquid ammonia is condensed 
on trimethyl gallium monoetherate, it dissolves readily and, upon evapora- 
tion of the ammonia, there is left behind a white solid containing one mole- 
cule of ammonia. The compound is sufficiently volatile at ordinary tem- 
peratures to permit of its transfer as a vapor. It is convenient as a start- 
ing material for the preparation of other compounds. 

(1) Preparation.—A sample of trimethyl gallium monoetherate, sealed 
in a suitable fragile bulb, is introduced into a reaction tube, and cooled in 
a bath of liquid ammonia. After breaking the fragile bulb, liquid ammonia 
is condensed in the reaction tube to dissolve the etherate. The ammonia 
is allowed to evaporate and the tube is pumped for some time through 
a liquid ammonia trap to remove all traces of ether. The ammine tends 
to form a supersaturated solution at room temperatures and it is usually 
necessary to induce crystallization by chilling. The product can be 
transferred to a suitable weighed tube by carrying the vapors at very low 
pressures through a tube cooled in liquid ammonia. 

(2) Physical Properties—Trimethyl gallium monoammine, (CHs)3- 
Ga-NHs, is a white, crystalline solid which melts at 31°. It sublimes in 
a vacuum at room temperatures. It is readily soluble in ether and liquefies 
in ammonia vapor at atmospheric pressure and room temperature. It 
is fairly stable toward oxidation but is quite sensitive to water vapor. 

(3) Hydrolysis and Analysis.—Trimethyl gallium monoammine dis- 
solves readily in an aqueous solution of potassium hydroxide, under which 
condition it is hydrolyzed with evolution of one molecule, each, of methane 
and ammonia. Samples of the ammine, sealed in weighed, fragile bulbs, 
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were treated with 7.5 N potassium hydroxide. The gases produced in 
the reaction were carried through standard sulphuric acid to absorb the 
ammonia and the remaining gas was collected over mercury by means of 
a Toepler pump. The residue remaining in the reaction tube was ana- 
lyzed for gallium. 

(4) Analysis for Nitrogen: substance, 0.1619, 0.2141; cc. 0.09367 N 
H2SO,, 13.15, 17.38; % N found, 10.66, 10.65, mean 10.65; % N required 
for (CH;)3Ga-NHs, 10.63. 

(5) Analysis for Gallium: substance, 0.1619, 0.2141; GayO3, 0.1152, 
0.1525; % Gafound: 52.93, 52.98, mean 52.95; % Ga required for (CHs)s- 
Ga-NHs, 52.89. 

(6) Methane Determination: m. mols (CH3);3Ga:-NHs3, 1.23, 1.63; cc. 
of gas, 27.4, 56.9; mol. wt., 16.3, 16.3; m. mols. CHy, 1.22, 1.64. 

Hydrolysis evidently occurs according to the equation: 

(CH3)3Ga:N Hs + KOH = (CH3)2GaOK + CH, + NHs3. 


1 DuPont FELLOW in Chemistry. 

2 Dennis and Patnode, Jour. Amer. Chem. Soc., 54, 182 (1932). 
> Renwanz, Ber., 65, 1308 (1932). 

4 Hein, Z. anorg. allgem. Chem., 141, 212 (1924). 


CHLORINATION PRODUCTS OF TRIMETHYL GALLIUM 
By CHARLES A. KRAUS AND FRANK E. TOONDER! 
CHEMICAL LABORATORY, BROWN UNIVERSITY 


Communicated January 31, 1933 


The fact that trimethyl gallium hydrolyzes in an aqueous solution of 
potassium hydroxide with the formation of one molecule of methane at 
ordinary temperatures and two molecules at higher temperatures? indi- 
cates that the methyl groups are readily replaceable by more electronega- 
tive elements or groups. This behavior is characteristic of the elements 
of the third group. It was thought, therefore, that chlorination products 
of trimethyl gallium might be prepared by treating the gallane with hy- 
drogen chloride under controlled conditions. Trimethyl gallium ammine, 
treated with an excess of hydrogen chloride at temperatures up to 160°, 
is completely chlorinated, the three methyl groups going to form methane, 
while ammonium tetrachlorogallate remains behind in crystalline form. 
This reaction was employed for determining the methyl content of tri- 
methyl gallium. 

The reaction actually takes place in several stages. When trimethyl 
gallium ammine is treated with one molecule of hydrogen chloride in ether, 
finely divided ammonium chloride is precipitated, which then reacts with 
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the gallane to form dimethyl] gallium chloride monoammine (CH3)2GaCl-- 
NHs, and one molecule of methane, which is evolved. On treating the 
ammine in ether with two molecules of hydrogen chloride, dimethyl 
gallium chloride and a molecule of ammonium chloride are formed; 
on treating the ammine under the same conditions with three molecules 
of hydrogen chloride, a second methyl group is substituted and the final 
products are methyl gallium dichloride and ammonium chloride. The 
dichloride forms a monoammine as well as a pentammine. The following 
equations describe the reactions: 


(1) (CH;);Ga-NH; + HCl = (CH;);Ga + NHC 

(2) (CH3)3Ga ob NH.Cl = (CH3)2GaCl-NH3; + CH, 
(3) (CH3)eGaCl-NH3 + HCl = (CH3)2GaCl + NH.Cl 
(4) (CH3)2GaCl a HC] = CH;GaCl. + CH, 

(5) CH;GaCk + HCl = GaChk + CH, 

(6) GaClk + NH,Cl = NH,GaCh 

(7) (CH3)3Ga-NHs3 ++ 4HCl1 = NH,GaCh +- 3CH, 


We have been able to identify the ammoniates of the various chlorine 
derivatives of trimethyl gallium. 

Ammonium Tetrachlorogallate—(1) Preparation. Weighed samples of 
trimethyl gallium ammine were prepared in the weighed reaction tube 
A of the accompanying figure by the method described in an earlier paper.” 
The reaction tube was connected to the system at B by means of deKhotin- 
sky cement and was disconnected for weighing as needed. Hydrogen 
chloride was passed from a generator attached at C into the calibrated 
receiver D, the pressure being determined by means of a manometer 
attached at E. On cooling the reaction tube A in liquid air and opening 
stopcock F, the hydrogen chloride was condensed in the reaction tube. 
Any desired quantity could be introduced by this method. On allowing 
the reaction tube to warm up to room temperatures, methane was evolved. 
The methane was determined by cooling the reaction tube in liquid air 
and pumping the methane through G into a receiver by means of a Toepler 
pump. On addition of two molecules of hydrogen chloride, approximately 
one molecule of methane was obtained. Excess hydrogen chloride was then 
introduced and the reaction mixture heated to 160°. A crystalline solid 
resulted and a total of three molecules of methane were recovered. The 
product in the reaction tube appeared to be homogeneous, did not sublime 
in a vacuum at a temperature of 200° and was evidently ammonium 
tetrachlorogallate. This was confirmed by the weight of the product in 
the reaction tube. The results of two preparations follow: (CH;)sGa-NHs, 
0.2108, 0.4468; mol. wt. gas, 16.1, 16.1; m. mols. CH, found, 4.66, 10.06; 
m. mols. CH, calculated, 4.80, 10.17; weight product found, 0.3671, 0.7787; 
weight required for NH,GaCh, 0.3672, 0.7783. 
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The weight of the final product corresponds closely with that required 
for ammonium tetrachlorogallate. The methane recovered was slightly 
low in one case due, without doubt, to the difficulty of recovering the 
methane completely from the tube which, at the time, was cooled in liquid 
air. 

(2) Analysis of Ammonium Tetrachlorogallate——Weighed samples of 
ammonium tetrachlorogallate were dissolved in water, acidified with nitric 
acid and treated with a slight excess of silver nitrate to precipitate silver 
chloride. The precipitate was filtered, washed and weighed. 

(3) Chlorine, substance, 0.3672, 0.7787; AgCl, 0.9192, 1.9370; % Cl 
found, 61.90, 61.54, mean 61.72; % Cl required for NH,GaCh, 61.76. 
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The filtrate was treated with hydrochloric acid to precipitate the trace 
of silver present and the gallium was precipitated as hydroxide with 
tannin in the presence of acetic acid.* After washing, the precipitate was 
ignited and weighed as usual. 

(4) Gallium, substance, 0.3672, 0.7787; GazO3, 0.1505, 0.3195; % 
Ga found, 30.49, 30.52, mean 30.50; % Ga required for NHiGaCh, 30.37. 

Ammonium tetrachlorogallate can be heated to 200° without melting 
and without sublimation. It is readily soluble in water without apparent 
hydrolysis. It seems to resemble the.corresponding compounds of boron, 
aluminum and thallium. 

Dimethyl Gallium Chloride Monoammine and Diammine—(1) Prepa- 
ration. Dimethyl gallium chloride monoammine was prepared by treat- 


FIGURE 1 
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ing trimethyl gallium ammine with two molecules of hydrogen chloride 
in ether solution. The apparatus employed was that already described. 
Initial reaction takes place between the hydrogen chloride and ammonia 
with the formation of ammonium chloride, which appears as a fine suspen- 
sion in the ether solution when the reaction tube warms to room tempera- 
ture, after having been cooled in liquid air for introduction of the hydro- 
gen chloride. Reaction occurs between the suspended ammonium chloride 
and the trimethyl gallium in solution, methane being evolved. At room 
temperatures, reaction is completed at the end of 48 hours, the end-point 
being indicated by the complete disappearance of ammonium chloride. 
The ammonia, resulting from the disappearance of ammonium chloride, 
combines with the dimethyl gallium chloride as ammine. The mono- 
ammine absorbs a second molecule of ammonia at ordinary temperatures 
to form a diammine. To demonstrate that the reaction proceeds through 
the initial formation of ammonium chloride, a reaction was carried out 
in which ammonium chloride was precipitated in ether by addition of equiva- 
lent amounts of ammonia and hydrogen chloride and then allowing this 
solid to react with a known sample of trimethyl gallium etherate. The 
reaction proceeded in precisely the same way as when the ammine was 
treated with hydrogen chloride. Trimethyl gallium also reacts with 
coarsely crystalline ammonium chloride, although the reaction is much 
slower than with the finely divided precipitate. No reaction occurs be- 
tween ammonium chloride and trimethyl gallium in liquid ammonia. In 
the absence of ammonia, hydrogen chloride reacts with the etherate in 
ether solution, the products being dimethyl] gallium chloride and methane. 

Direct analysis was not made of the final product, but the weights of 
the initial materials and of products were determined. The results of 
three such determinations are given in the following table. In no. 1, 
trimethyl gallium ammine served as starting material, while in nos. 2 
and 3, the etherate was used. The product of reaction with etherate was 
treated with ammonia to convert to the ammine. 


TABLE 1 
DIMETHYL GALLIUM CHLORIDE PREPARATION 
no. 1 no. 2 no. 3 

wT. M. MOLS. wT. M. MOLS. wT. M. MOLS. 
(CH;);Ga:NHs; 0.7417 5.627 1.0278 5.45* 0.4728 2.503* 
HCl 126.1lcc. 5.63 122 ce. 5.45 55.9cce. 2.49 
CH, 126.4 cc. 5.64 124 cc. 5.53 56.lcce. 2.50 
(CHs)2GaCl-NH; O:Siue O08: eeces ‘the 0.3811 2.50 
(CH3)2GaCl-2NH; 0.95388 5.63 0.9268 5.47 0.4231 2.50 


* Trimethyl gallium monoetherate. 


(2) Properties—Dimethyl gallium chloride monoammine is a white 
solid melting at 54°. -It is not appreciably volatile in a vacuum at 25°. 
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Dimethyl gallium chloride diammine is a white solid melting at 112°. It 
loses one molecule of ammonia when pumped at 60° At 25° the dis- 
sociation pressure is between 0.3 and 0.5 mm., and at 40° it is 1.0 mm. 
The compound is stable in air but is hydrolyzed by water vapor. The 
diammine is readily soluble in liquid ammonia and somewhat soluble in 
alcohol and acetone. It is insoluble in ether, in which respect it differs 
from the monoammine, which is readily soluble. Thus, the diammine 
may be precipitated from an ether solution of the monoammine by means 
of ammonia vapor. 

Methyl Gallium Dichloride-—(1) Preparation. When trimethyl gal- 
lium ammine is treated with three molecules of hydrogen chloride in 
ether solution, methyl gallium dichloride and ammonium chloride are 
formed. Initially, ammonium chloride is formed, which reacts to form 
dimethyl gallium chloride monoammine and this, in turn, reacts with hy- 
drogen chloride to form ammonium chloride and methyl gallium dichloride 
according to equations (1), (2), (3) and (4). When the etherate is treated 
with two molecules of hydrogen chloride in ether in absence of ammonia, 
methyl gallium dichloride is formed. 

These reactions were carried out according to the procedure already 
described. Direct analysis was not made of the compounds, but the 
weights of the materials were determined. Methyl gallium dichloride 
forms a monoammine and a pentammine. The results of several prepara- 
tions are given in the following table. The starting substance in the first 
two experiments was trimethyl gallium monoammine and in the third, 
trimethyl gallium etherate. In this reaction, ammonia was added after 
completion of halogenation. 


TABLE 2 


METHYL GALLIUM DICHLORIDE PREPARATION 


no. 1 No. 2 No. 3 
wt. M. MOLS. WT. M. MOLS. WT. M. MOLS. 
(CH;);Ga-NH: 0.3296 2.501 04000: S778. 6a cs 4.62" 
HCl 167.lcc. 7.46 .252.1 cc. 11.25 207.6cc. 9.24** 
CH, 111.0cc. 4.96 166.6cc. 7.44 206.0cc. 9.20 
5 SS Sh ae a ema sc gry eS i Se aoisks ee 1.99 
CH;GaCk + NH,Cl 0.5282 2.53 ibe: | Slee: Ser,” Sierras sae 


CH;GaCl.NH; + NH,Cl 0.5714 2.53 OGG tt os a es rats 
CH;GaCl.5-NH; + NH,Cl 0.7325 2.48 1.1064 3.76 1.11382 4.62 


* The etherate was used. ** The etherate requires 2 molecules of HCl. 


Treated with three or two equivalents of hydrogen chloride, respec- 
tively, the ammine or the etherate is converted to methyl gallium dichlo- 
ride with formation of two molecules of methane. 

(2) Properties.—Methyl gallium dichloride is a white solid which does 
not melt or sublime in a vacuum up to 75°. It is readily soluble in ether. 
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With ammonia it forms a monoammine- and a pentammine. Methyl 
gallium dichloride monoammine does not melt up to 80°. Its vapor pres- 
sure is below 0.01 mm. at 25°. The pentammine is insoluble in liquid 
ammonia and has a dissociation pressure of about 25 mm. at room tem- 
perature. The ammonia is readily pumped off at 80°. 


1 DuPont FELLow in Chemistry. 
2 These PROCEEDINGS, 19, 292-298 (1983). 
3 Moser and Brukl, Monatshefte f. Chem., 50, 186 (1928). 


THE THREE-ELECTRON BOND IN CHLORINE DIOXIDE* 


By L. O. Brockway 
GaTES CHEMICAL LABORATORY, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated February 13, 1933 


The stability of certain molecules! containing an odd number of electrons 
has been explained on the basis of the formation of three-electron bonds. 
Evidence has now been obtained which substantiates the validity of a 
similar explanation for the stability of ClO2, and in particular for its lack 
of tendency to polymerize to ClQ,. 

The possible electronic structures which may be written for ClO: in 
accordance with the quantum mechanical rules for the formation of 
electron-pair bonds*® are represented as follows: 


Oc ere: i! 
- + 
<O:CI: 0: II 
4- — 
‘OO: Gio: III 


Only the electrons in the incomplete shells are indicated. The unsym- 
metrical arrangement ClI—O—O is not considered, since it is highly im- 
probable for chemical reasons. The formal charges indicated for the atoms 
result from splitting shared electron pairs. 

‘The choice of the correct structure among these could be made by a 
theoretical calculation of the energy of each. If the energy for a single 
one lay considerably below the energies calculated for the other two, that 
structure would closely represent the molecule in its normal state. If the 
separate energy calculations gave nearly the same result, the true value 
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of the energy would be obtained only by considering a combination of all 
the structures. Normal chlorine and oxygen are rather near to each other 
on the electronegativity scale,* and the electronegativity of positive chlo- 
rine is greater than that of neutral chlorine, and may be nearly equal to 
that of oxygen, so that the transfer of one electron from oxygen to chlorine 
in ClO, may well involve only a small change in the energy of the molecule. 
Structures I, II and III above would then correspond to energy levels 
near to one another. 

The criterion for the formation of a three-electron bond is that two 
structures, in one of which there is an unshared pair of electrons on one 
atom and a single electron on the other, and in the other the single electron 
and the pair are interchanged, have nearly the same energy. It has been 
shown! that in this case a lower (and hence more nearly correct) value for 
the energy is obtained by taking a linear combination of the wave func- 
tions representing the two structures than corresponds to the use of either 
function alone. This extra resonance energy which serves to stabilize 
the molecule may be ascribed to the formation of a bond with the three 
electrons in the same way that the interchange of two single electrons on 
two atoms affords a resonance energy which is the principal energy of the 
electron-pair bond. The formation of the two kinds of bonds may be 
indicated diagrammatically in the same manner, that is, by placing two 
dots between the atoms for the electron-pair bond and three dots for the 
three electron bond. 

Since structures I and II or I and III taken together meet the conditions 
for the formation of a three-electron bond, a better representation of 
ClO, is given by the two formulas: 


:0:ClT0: IV 
:0'7C1:0: V 


In both IV and V one oxygen atom is bonded to chlorine by an electron- 
pair bond and a three-electron bond, and the other by an electron-pair 
bond alone. Since the calculation for the helium molecule-ion He.* to 
which reference has already been made indicates that the energy and 
interatomic separation of a three-electron bond is about equal to that of a 
one-electron bond and about half that of an electron-pair bond, the combi- 
nation indicated here in ClO, should be equivalent to one and one-half 
single bonds or about midway between a single and a double bond. - 

It was pointed out in a discussion‘ of resonance between several Lewis 
electronic structures that the observed bond distances between two atoms 
corresponds to that of the strongest type of bond involved in the resonating 
structure and not to a mean (if the elementary structures are of equal 
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importance). Since the structures IV and V are indistinguishable, they 
have precisely the same energy and therefore contribute equally to the 
resonating structure. On this basis the observed chlorine-oxygen sepa- 
ration will be the same for both bonds and should correspond to the com- 
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FIGURE 1 


bined electron-pair-three-electron bond distance. The predicated value 
of 1.57 A for this distance is based upon the properties of the combined 
bond mentioned above and on the table of covalent radii.‘ The single 
and double bond distances for CI—O are 1.65 A and 1.48 A, respectively. 
I have determined the bond distance in ClO, with the aid of electron- 




















TABLE 1 


MODEL II MODEL III MODEL IV MODEL V MODEL VI 


MODEL I 


SIN 0/2 


150° 125° 110° 90° 50° 


180° 


r 
OBSERVED 


1.338 
. 589 
1.650 


4.63 


4.25 1.228 


7.72 
11.20 
14.15 
17.05 


1.257 
1.592 
1.597 
1.574 


4.35 
7.92 
10.80 
13.90 
17.40 


1.330 
7.80 1.569 


4.60 
10.65 1.575 
14.30 
17.10 


4.75 1.37 


1.364 
1.518 
1.648 
1.568 
1.587 
1.580 
0.037 


4.7. 


0.276 
0.396 


0.538 
0.703 


0.875 


First min. 


1 


.90 
11.15 
14.30 
17.30 


~ 


1.551 
1.657 
1.601 
1.55 


7.60 1.527 


7.55 
Lis 
13.85 
17.45 


First max. 


1.645 
1.574 


11.12 
13.90 
17.40 


Second min. 





1.619 


1.619 


Second max. 


t 
ve) 


1.583 
1.587 
0.008 


1.556 
1.580 
0.020 


1.583 
1.582 
0.032 


Third min. 


1.608 
0.026 


1.590 
0.039 


Mean value 


Average deviation 
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= 1.58 + 0.03 A. 


Final value Cl—O 
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diffraction photographs of the gas. The 
method was the same as that used in the 
study of the hexafluorides.’ Theoretical 
intensity curves for six molecular models 
corresponding to a range of bond angles 
from 50° to 180° are reproduced in figure 1. 
It will be noted that the positions of the 
maxima and minima change very slightly 
throughout the range. This is explained 
by the large scattering power of chlorine 
relative to that of oxygen; the interference 
effects observed arise almost exclusively 
from the chlorine-oxygen separation so 
that the effect of the change of the 
oxygen-oxygen separation with variation 
of the bond angle is scarcely detectable. 
In table 1 are given the results of the 
observations as interpreted with the use 
of each of the theoretical models. 
The first column contains the anc be 
values for the observed minima and max- 
ima, in which 6 is the scattering angle and 
\ the wave-length associated with the elec- 
tron beam. ‘Each value is the average of 
six measurements on three different photo- 
graphs. The remainder of the table shows 
the interatomic distance (a = Cl-O separa- 
tion) as calculated for each minimum and 
maximum of all the models with the aid of 
the equation 

sin 6/2 

“Ee 





x = 4na 


The mean values for a do not include those 
calculated for the first minimum since the 
position of the first minimum on the 
photograph is always affected by the fogged 
area around the central image. 

The most probable model is considered 
to be the one giving the most consistent 
a-values. Since in this case the deviation 
in the least consistent set is only a little 
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larger than the experimental error (which is + 1%) no very definite con- 
clusion in regard to the bond angle is afforded by these results. An angle 
near 120° is expected from theoretical considerations. The uncertainty 
in the value of the angle does not affect the determination of the Cl-O 
separation, however, if the extremely improbable Model VI is excluded. 
The most probable value for the Cl—O distance is thus found to be 1.58 
+ 0.03 A. 

The excellent agreement between the observed value and the value 
predicted for the structure discussed above (1.57 A) substantiates this 
electronic structure for ClO., and provides positive evidence for the 
existence of the three-electron bond in this molecule. 

I wish to thank Professor Linus Pauling for many invaluable suggestions 
and criticisms. 

* Contribution from the Gates Chemical Laboratory, California Institute of Tech- 
nology, No. 358. 

1 Linus Pauling, J. A. C. S., 53, 3223 (1931); J. Chem. Phys., 1, 56 (1933). 

2 Linus Pauling, J. A. C. S., 53, 1367 (1931). 

3 Linus Pauling, [bid., 54, 3570 (1932). 

4 Linus Pauling, Proc. Nat. Acad. Sci., 18, 293 (1932); Ibid., 18, 498 (1932). 

5 L. O. Brockway and Linus Pauling, Jdid., 19, 68 (1933). 


THE NEUTRON, ATOM BUILDING AND A NUCLEAR 
EXCLUSION PRINCIPLE 


By WILLIAM DRAPER HARKINS 
GEORGE HERBERT JONES CHEMICAL LABORATORY, UNIVERSITY OF CHICAGO 


Read before the Academy, Monday, November 14, 1932* 


I. Introduction: Discovery of the Neutron.—The recent experiments 
of Curie and Joliot,® based on those of Bothe and Becker,’ have been 
shown by Chadwick’ to indicate the existence of a new fundamental unit 
of structure, the neutron, which is undoubtedly of great importance in 
the building up of atomic nuclei. 

Neutrons were found to be emitted by beryllium when it is bombarded 
by fast a-particles. That the beryllium nucleus consists of a neutron, 
or a complete condensed hydrogen atom in a nuclear state was an integral 
part of the hydrogen-helium theory of nuclear constitution developed in 
1915 by Harkins and Wilson.' The prediction of the existence of a neutron 
in the beryllium nucleus arose from a curious numerical discrepancy which 
will be described below. The ordinary beryllium nucleus was considered 
as a carbon (12) nucleus in which one helium nucleus is replaced by a 
neutron. It is now found that if the a-particle enters the beryllium 
nucleus the neutron leaves it. 
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According to the theory light nuclei or even charge consist in general of 
a-particles alone. That is, the composition of each atom is that of a whole 
number of helium atoms as follows: 


Element He Be & O Ne Mg Si $s 
Atomic number 2 + 6 8 10 12 14 16 
Atomic weight predicted 4 8 12 16 20 24 28 32 
Atomic weight determined 4 9 12 16 20 24 28 32 
Composition of atom He HesxH He; He, Hes Hes Her Hes 


The only exception to the general theory was found to be beryllium 
and for this it was necessary to assume the presence of a complete con- 
densed nuclear H particle, that is, proton plus electron, or neutron in the 
nucleus. The formula of the beryllium nucleus was written a2(7f) in 
a paper? submitted early in 1919. Here m represents a proton, 8, an 
electron and the parentheses indicated the union into a neutron, or par- 
ticle of zero charge. Paired neutrons (7282) were also assumed to exist 
in many atomic nuclei. 

Following this the existence of neutrons, or atoms of zero atomic number 
and zero nuclear charge as separate particles was postulated by Harkins? 
(April 12, 1920) and by Rutherford‘ (June 3, 1920). 

The writer was led to this prediction by (a) the apparent existence of a 
neutron in the beryllium nucleus, (b) of a double neutron (7262)° in all of 
the ordinary members of the neon?*—uranium series, and by (c) the 
apparent necessity of the existence of neutrons to explain the formation 
of nuclei of high positive charge, as uranium, which exert a high repulsion 
toward a-particles or protons, but through whose region of repulsion 
neutrons could easily pass. 

A few quotations from this 1920 paper present some of these ideas. 


1. Attention is called to the suggestion that atoms of zero atomic number may 
exist and be of importance in atom building. Such atoms might have masses 4, 3, 2 
and 1, and possibly other values, and they would contain no non-nuclear electrons, so 
they would have no chemical, and almost none of the ordinary physical properties, 
aside from mass 

2. Suppose that several a-particles are approaching each other. Their mutual 
self-repulsion due to their net positive charges, tends to keep them apart unless they 
are driven together so strongly that they approach sufficiently close to enable the 
attractive effects of the mutual couples to become of significance. As the total net 
positive charge on the nucleus, once it is formed, increases, the repulsion due to this 
cause may be assumed to increase to such an extent that a-particles are no longer able 
to pass through the region of repulsion, and to attach themselves, though electrically 
neutral particles. ... could easily pass into and through this region. 

3. It is not improbable that some isotropic atoms are formed by the addition of the 
group (78)°. , 

4. The atoms of the meta-neon-uranium series seem to have nuclei which consist 
of a-particles mostly, but which contain in addition one other group. Since the nature 
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of this group has not been determined experimentally, its exact size is uncertain. How- 
ever, the difference between a member of the a- or helium series and a corresponding 
member of the meta-neon-uranium series consists in the presence in nuclei of the latter 
type of a group consisting of 2 positive and 2 negative electrons, of the formula (7262)°, 
where 7 represents a positive and 6 a negative electron and 0 indicates that the net 
charge on the group is zero. Jf this group exists by itself, it is a particle belonging to an 
element of zero atomic number. 


While the question of stability is not discussed in this paper, it may be 
stated that at that time the writer did not consider any neutron of higher 
mass than 2 to exist more than momentarily except in a moderately 
heavy nucleus, but that the first quadruple neutron present in other 
nuclei is that in argon 40, whose nuclear formula was written a(78), or 
ag( Be). 

II. Neuton, the New Element.—One of the objects of this paper is to 
point out that the neutrons in our universe constitute a new element of 
atomic number zero, which, on account of the neutral character of the 
atomic nuclei of which it consists, may be called neuton. 

This element is undoubtedly of fundamental importance in the formation 
of other elements. Since its atoms are very light and, unlike those of 
hydrogen, are not held to the earth in the form of chemical compounds, 
it should be found in space with a moderately high relative concentration, 
that is, as compared with the other elements in space. 

Since its effective thermal collisions are mostly with nuclei, rather than 
with atoms, the (at present theoretical) study of the distribution of thermal 
energy between neuton and any gases, liquids or solids with which it is 
mixed is of considerable interest. 

III. Emission of y-Rays by Nuclei Excited by Neutrons, and Nuclear 
Energy Levels.—Neutrons seem to be remarkably efficient in the excita- 
tion of atomic nuclei with which they unite. When nitrogen, by the im- 
pact of a neutron is converted into boron and helium, mass may be lost 
(energy equivalent to 1.4 million electron-volts provided Aston’s mass data 
are entirely correct) and converted into y-rays or it may be nearly con- 
served. In the disintegrations with capture of the neutron which have 
thus far been investigated, kinetic energy has either been conserved, or 
a part of it has disappeared and its energy has appeared in the form of 
-rays. 

While examining half a dozen photographs of the disintegration of nitro- 
gen obtained by Gans, Newson and the writer,’ two remarkably good 
(A) photographs were found which represent capture of the neutron. 
These are marked (HA) in table 1, where they are compared with results 
from three of the best photographs obtained by Feather® (FA), together 
with four (FB) from slightly imperfect, and five less perfect (FC) photo- 


graphs. 
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TABLE 1 


AMOUNT OF ENERGY (IN 10° ELECTRON-VOLTS) WHICH DISAPPEARS INTO THE REACTION 
AND Is PRESUMABLY CONVERTED INTO y-Rays 
N“ + ni —> nit — > BU + Het 
1 1.4 (HA) i (FB) 1.3 (FC) 
2a 2.3 (HA) . (FB) 2.0 (FB) 2.3 (FC) 


8 (FA) .5 (FC) 
3 3.8 (FA) 3.8 (FA) 3.9 (FB) 
4 5.0 (FC) .8 (FC) 


If the rest masses alone of N'‘ and n! are considered the rest mass of the 
newly formed N® is 15.0147 + 0.004, which undoubtedly represents an 
excited state of this nucleus. In addition there has been 0.002 or more 
mass units of kinetic energy to dispose of in the cases which have thus far 
been investigated. The excess energy is presumably emitted as a y-ray. 

The simplest explanation of the constancy of the lowest value, 1.4 x 
10° e.-v. as found in the table is that this value is due to about this amount 
of error in Aston’s mass values. 

On the basis of the last paragraph the amount of energy converted into 
y-rays in the disintegration process, as given in table 1), strongly suggests 
differences between definite energy levels, although not quite so strongly 
as table 1. It is obvious that more accurate mass data must be secured 
before it can be determined whether mass disappears, is conserved, or is 
gained in this nuclear reaction. 


TABLE 1B 
ENERGY WHICH DISAPPEARS IN 10° ELECTRON-VOLTS 


0 0 
0.6 0.6 
: | 
3 i 2.4 2.5 
4 . 4.4 


The question arises as to whether photographs of such events taken in 
the Wilson Cloud Chamber, might be used to show from just what point 
on the tracks the photon is emitted. 

A photon of energy 1.5 X 10° e.-v. has a momentum of 2.1 X 107g. 
cm. sec., while either the B'! or the He* atom has a very much larger 
momentum, e.g., 4.63 and 5.41 X 10~ for the boron nuclei represented 
by 1 (1.5A) and 3 (3.9A) in table 1. In the former, the momentum of the 
photon is 1.9 per cent, and in the latter 3.8 per cent of that of the boron 
nucleus, while the percentage value probably rises even higher for other 
events. It is probably not impossible to detect the position of emission 
of a y-ray of the maximum momentum provided it occurs in both a favor- 
able position and direction. 

The photographs of the disintegration of nitrogen by neutrons indicate 
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that the N” nucleus moves less than 0.15 mm. before it disintegrates. 
This gives an upper limit for its life as 10—” seconds. 

In the reaction, O* + n! —> O!’—> C!3 + He‘ as studied by Feather, 
mass is very nearly conserved, so only the kinetic energy of the neutron 
is (partly) available for transformation into radiation. 

It should be noted that the nitrogen disintegrations, for which 1.5 x 
10° e.-v. of energy in the form of mass may be available for the trans- 
formation, occur much more often than those of oxygen, for which only 
the kinetic energy of the neutron is used to disrupt the nucleus. 

The present status of the experimental work on the disintegrative syn- 
thesis of atoms or disintegration by capture may be summarized as follows: 


Periodic System of 
the Atomic Species 


§ 
5 
| 
r 4 
2 
g 
3 


32 
3 
x» 
29 
28 
27 
26 
25 
24 
23 
22 
2i 
20 
(9 
18 
17 
16 
15 
14 
13 
12 
i" 
10 
9 
8 
7 
6 
5 
4 
3 


o-n 


be Zou SS FaHand<xBhr>SELSZINSS LISLE IEGHL-LISSBISE 
7 RANQSRSRSARHABABHHRRQSIVSI VQ VIZGASSILSAIS 
Atomic Number 
FIGURE 1 
The number of extra protons or of hemi-alpha groups is given on the x-axis and the 
number of extra neutrons on the y-axis. The sum of these gives the total number of 
neutrons in the nucleus. 


1. a-Particle of double charge and small excess mass. Kinetic energy 
seems to be transformed into mass (possibly into y-radiation). 

2. Proton of single charge and high excess mass. Excess mass seems 
to be converted almost entirely or entirely into kinetic energy, or partly 
into kinetic energy and some other form of energy (as y-rays). 

3. Neutron of zero charge and high excess mass. Excess mass may 
disappear and presumably escapes as y-rays. Kinetic energy may also be, 
and usually is, converted into y-rays. 

Much work is needed to show other relations and to give a more general 
picture of disintegrative syntheses. Thus far when an a-particle has gone 
into a heavier nucleus a proton or a neutron has come out, and when a 
proton or a neutron has gone in an a-particle has come out. 
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IV. Structure of the Nucleus—Two general theories of the constitution 
of atomic nuclei have been proposed: 

(1) The nucleus is built up largely of a-particles or groups (Harkins 
and Wilson!), and 

(2) The Harkins-Masson'’™" nuclear formula, which may be written 


(pre) z (pe), (1) 
(np) zn; (2) 


(2) ® 


nz (4) 
az/2(Me)1/2 (5) 


represents the constitution of any nucleus. The possibility that this 
formula might be made the basis of a theory of the nucleus was pointed 
out by the writer to Heisenberg" in 1929, and he has used it as the basis 
of his recent nuclear theory. (See also the earlier paper of Iwenenko.!”) 

In this paper p or 7 represents a proton; e or 8, an electron; nor w,a 
neutron; yw, a double neutron m; a/2 or 7, a hydrogen nucleus of mass 2 
or the corresponding hemi-alpha group; a, an a-particle or group; Z, 
the atomic number, the number of hemi-alpha groups or the number of 
n or extra protons; and J, the isotopic number, or number of extra neutrons. 
The total number of neutrons is N = Z + I. 

Form (5) of the formula indicates that in general, with few exceptions 
the hemi-alpha (n) groups are present as a-particles. 

Since every neutron in a nucleus is supposed to contain an electron, 
and since on the basis of formula (1b) no electron is present except as a 
neutron, it is obvious that any relation which exists concerning the number 
of electrons holds for the number of neutrons as well. 

The most important nuclear relations found by the writer are expressed 
below in terms of the Harkins-Masson formula. 

1. No complex nucleus contains more than two protons to one electron, 
that is, not more than one proton can combine with a single neutron. 

2. The general stability of atoms nuclei is determined by two variables, 
Z, the number of extra protons, and N/P the ratio of the number of neu- 
trons to the total number of protons. A mathematically equivalent state- 
ment is that the stability depends upon Z, and N/Z, or the ratio of neutrons 
to extra protons. Let 
N 


ee Bi okeisinee 


be taken as the dependent variable. Then all stable atoms lie in a band 
(Fig. 2) which is broad for low values of Z, and narrows as Z increases. 
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In the radioactive region 6-changes start if the top of the band is too closely 
approached. That is where y or \ approaches the maximum limit, and 
a-changes begin if y or \ becomes too small, or approaches the minimum 
limit. The fundamental nature of this relation pointed out in 1921 by 
the writer has just been recognized by Heisenberg, who uses it as the most 
important in his theory. 
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Thus for any relation which is valid for y there is also a mathematically 
equivalent relation which holds for }. 

3. The number of neutrons is even for almost every known atomic 
nucleus. Thus 99.2 per cent of the atoms in the earth’s crust contain 
an even number of neutrons. This is shown by figure 3, which gives 
the atomic percentages of all of the known atomic species with the excep- 
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Abundance of atomic species as a function of the neutronic number 
(total number of nuclear neutrons). 


tion of those too rare to be visible on the plot. The height of the peaks 
in this plot is remarkable if it is taken into account that that for oxygen 
would be about 1000 feet high on a linear scale if the points for 3 and 5 
neutrons have the heights given in ‘the figure. 

4. The hemi-alpha (n) particles are in general associated in pairs, 
that is, they exist as a-particles or groups. It is obvious that an odd 
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hemi-alpha group, such as that in nitrogen cannot be part of an a-particle, 
but such odd or relatively free hemi-alpha groups seldom occur. 

It may be noted that on this basis the number of a-groups in any nucleus 
is one-half the atomic number (Formula 5). 

5. Relations 3 and 4 indicate that the total number (P) of protons in 
a nucleus is usually even. Thus P = 2Z + I = N + J is usually even. 

V. A Nuclear Exclusion Principle—An important nuclear relation 
becomes apparent upon the adoption of the system of classification ex- 
hibited below. 

All atomic nuclei may be divided into four classes, each of which should 
be represented by 25 per cent of all atoms, provided inside the region of 
stability the numbers of neutrons and protons which combine are de- 
termined entirely by chance. 


TABLE 2 
TOTAL EXTRA TOTAL PERCENTAGE OF 
NEUTRONS NEUTRONS EXTRA PROTONS MATERIAL 
N I PROTONS P EARTH METEORITES 
Class I Even Even Even Even 87.4 95.4 
Class II Even Odd Odd Odd 10.8 S31 
Class ITI Odd Odd Even Odd 1.8 2.5 
Class IV Odd Even Odd Even 0.007 0.0 


Nuclei of class IV may be divided into sub-classes (a) and (6). Only 
4 atomic species of class IVa are stable, and these represent rare species. 
No atoms of class IVb are stable. 

Since atoms of class IV are 13,000 times more rare than those of class I, 
it is apparent that nuclei for which both the number of neutrons and of 
extra (or eta) protons, are odd, are almost excluded from existence. 

However, the relation becomes more remarkable when it is found that 
all of the 0.007 per cent of material which belongs to this class, follows 
the relation which follows. 

All of the species of class IV which are known to exist (H’, Li’, B” and 
N"™) have an isotopic number zero, their maximum atomic number is 7, 
and maximum mass number 14 and they belong to subdivision a of this 
class. 

Thus a nucleus az/2), which consists of a-particles alone, will take 
up a single hemi-alpha (n) particle to form a stable nucleus if the number 
of a-particles is three or less, but will not if the number is larger than an 
unknown lower limit whose value at present seems to be three. 

However, such nuclei will in general take up a neutron, and for larger 
values of Z/2, a double neutron (m2) or yu. 

The above relation is the basis for the assumption that the order of 
binding is p, n, n, p, to form the a-particle, and then n, p, n, p, n, p, n, 


p, n, p, n, p, n, n, P, p, n,n, p, , N,N, p, p, Nn, Nn, p, p, n,n, p. (Fig. 1.) 
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VI. Second Exclusion Principle-——No nucleus which contains an odd 
extra proton, that is an odd hemi-alpha (n) group is stable if the number 
of extra neutrons is even (class IVb), but is stable if I is odd (class II). 
This causes the nucleus to shift into class II. 

This principle indicates that for such a nucleus to be stable the neutron 
of the odd hemi-alpha (n) group must be paired with the odd extra neutron. 

Thus for such a nucleus to be stable the total number of its neutrons 
must be even, not odd. 

The existence of this relation gives additional evidence for the pairing 
of neutrons. 

However, since the spin of such a nucleus has values such as 3/2 or 5/2 
in units 4/27, the pairing is due in a primary sense to something other 
than spin. 

These relations seem to indicate the existence in atomic nuclei of the 


a Qa 
following groups: ; np, a/2 orn; m2 Or uw; Mp, 3” Or 7M; Nop. OF a. 


VII. Binding Energy of Nuclear Groups.—It should be of considerable 
interest to show the connection between the binding energy of nuclear 
groups and the relations presented in this paper. However, this would 
lead to too lengthy a discussion to be presented here. It is of interest 
to note that the binding energy in millions of electron-volts of two hemi- 
alpha (n) groups is 23, while that of a proton and electron to form a neutron 
is about one, and of a neutron and proton to form a hemi-alpha group is 
also about one. The addition of 7 to as (carbon nucleus) to give nitrogen 
gives 8.5, while the addition of a second 7 to form oxygen gives 20 million 
electron-volts. Thus the constituents of an odd hemi-alpha (n) group are 
bound much less firmly than those in an a-group. This and other facts 
agree remarkably well with the general point of view presented in this 
paper. 

It is not improbable that a theory intermediate between the two (1 and 
2) presented in Section 4 may in the end prove the most satisfactory. 

VIII. The Emission of Neutrons as Related to the Harkins-Masson 
Formula.—It is of considerable interest that the three atomic species 
which have been found to emit neutrons on bombardment of their nuclei 
with a-particles all have the nuclear formula 


with Z = 3, 4, 5, and the isotopic number J = 1. Thatis, the Harkins- 
Masson formula indicates the existence of one extra neutron in each of 
these three nuclei. : 

The relations are shown by figure 4, which gives the atomic number on 
the x-axis, and (J) the number of extra neutrons on the y-axis. 
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In each case an a-particle adds itself to the Li,’, or Be,® or B,!! nucleus, 
thereby increasing the nuclear charge and atomic number by 2, and form- 
ing B,!', C,!* or Ni%. The isotopic number of all of the above species 
is one as represented by the subscript. However, each of the elements 
B, C and N has a stable isotope of zero isotopic number. Since the 
newly formed nuclei B,"', C,!* or N,® possess a large amount of energy 
due to the kinetic energy of the a-particle, and the mass relations are not 
unfavorable, a neutron is emitted and the lower isotope is produced. 

This figure indicates that neutrons should also be emitted by a-ray 
bombardment of C;!4%, O,!7, Ne, and Si,”. 

On account of the higher nuclear charge somewhat faster a-particles 
are needed than with the lighter atoms. These species are not convenient 
sources for neutrons because each of them is a relatively rare isotope, 
though less rare for magnesium and silicon. 

However, with F,;”, Na,”, Al?” and P,*!, there is no known stable 
isotope of isotopic number zero, so the emission of a proton and not of a 
neutron is indicated, though it is not entirely precluded. 
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FIGURE 4 

Emission of neutrons by nuclear collisions of a-particles with cer- 
tain nuclei. 


Thus the emission of neutrons, instead of protons from Li,’ and Be," is 
due to the existence of two stable species of the extremely rare class IV, 
Beo” and No". 

Either figure 4 or the relations of Section 6 indicate that the energy of 
emission of the neutron should increase in the order: neutrons from Li’, 
B" and Be’, which is what is found by experiment. 


* Revised February 3, 1933. 
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BIOLOGICAL SIGNIFICANCE OF PROTECTIVE MECHANISMS 
INHERENT IN THE MYOCARDIUM 


By Cart V. WELLER 
DEPARTMENT OF PATHOLOGY, UNIVERSITY OF MICHIGAN 


Read before the Academy Wednesday, November 16, 1932 


To the clinician who sees the effects of disease in the living patient 
and to the pathologist who observes the evidences of disturbed function 
and altered morphology in the tissues removed by the surgeon or ob- 
tained at autopsy, it appears that all too often the most important struc- 
tures of the body are those which are most vulnerable. Since the serious 
morbific lesions and those which are lethal are the most impressive, it 
is easy to understand the basis for this assumption. Evidence can be 
adduced, however, that to some degree such necessary structures as the 
myocardium have acquired certain inherent protective mechanisms which 
are not equally operative in many other tissues of the body. Some of 
these protective factors have obvious mechanistic explanations, while 
others are due to inhibiting actions which are without adequate explana- 
tion at the present time. 

One of the most obvious of the protective mechanisms of the myo- 
cardium is dependent upon the pattern of the arterial blood supply. Not 
only is there a variable anastomosis between the branches of the right 
and left coronary arteries, but the number of these connections probably 
increases as age advances, so that the heart becomes better prepared to 
withstand occlusive coronary disease in that very period of life when 
such lesions are prone to occur. Also there is found in the myocardium 
a peculiar intermingling of arterioles derived from different arterial 
branches. Thus it is often seen that the necrosis of infarction following 
occlusive coronary disease will alternate in distribution with patches of 
intact muscle. This is especially well seen when such lesions reach the 
stage of fibrosis. Because of the interdigitation of two regionally inde- 
pendent arterial systems the myocardium is able to withstand the occlu- 
sion of one of the smaller coronary branches without danger of rupture 
or aneurysm formation. If the same amount of necrotic muscle or of 
scar tissue were concentrated in a single area death would be inevitable. 
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The myocardium enjoys comparative immunity from metastatic malig- 
nant disease. It is probable that no biological significance, in the sense 
of evolutionary modification, should be attached to the ability of the lungs 
to retain neoplasm cells which might otherwise find their way into the 
myocardium. This function of the lungs is obviously of minor import 
compared to that of aeration of the blood. However, even in those cases 
in which hematogenous skeletal and visceral metastases are present the 
myocardium almost always escapes. That this is not due to purely 
mechanical factors is indicated by the fact that with certain neoplasms, 
as, for instance, the melanoblastomas, the myocardium may contain a 
large number of secondary tumors. 

In infectious diseases similar phenomena may be present. Of all of the 
internal organs the aorta is most frequently the site of macroscopically 
and microscopically demonstrable syphilis. While both congenital and 
acquired syphilis of the myocardium are by no means rare, it is frequently 
noted at autopsy that acquired syphilitic aortitis with easily demonstrable 
spirochetes may extend to, and involve, the aortic valves without gross 
or microscopical evidence of invasion of the myocardium. It is by no 
means clear why there should not be an invasion of the heart muscle similar 
to that of the first portion of the aortic arch. It is true that the examina- 
tion of many blocks from such hearts will usually show active syphilitic 
foci. The limited extent of such foci, and particularly the paucity. of 
spiral organisms in them, is very often in sharp contrast to the condition 
in the aorta. That this is not an essential, or an unvarying characteristic 
of heart muscle is shown by the extraordinary involvement of the myo- 
cardium in many cases of congenital syphilis in which spirochetes may 
be found in tangled masses throughout the interstitial tissue. There are 
also examples of diffuse and active acquired syphilis of the myocardium 
with many spiral organisms, but such are comparatively rare. It is evi- 
dent that some protective mechanism must be operative which in many 
hearts restricts syphilitic disease to minimal, and more or less latent, foci. 

A striking example of the principle under discussion occurs in human 
trichinosis. Death from massive infestation with trichinae may be due to 
circulatory failure in the fourth to seventh week of the disease. In such 
cases the embryos in enormous numbers will be found in an early stage 
of encystment in the skeletal muscles. In the myocardium, however, 
not a single embryo can be,found, although a characteristic active intersti- 
tial myocarditis with numerous areas of lymphocyte and eosinophile cell 
infiltration is present. Necrotic heart muscle fibres also occur. While 
some have believed that this myocardial damage with its subsequent 
inflammatory reaction is of toxic origin, there is experimental evidence 
that the embryos reach the heart muscle during the period from the eighth 
to the twelfth day after ingestion just as they do the voluntary muscle. 
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By the fourteenth day inflammatory foci are present but the embryos 
can now be found only with great difficulty, if at all. Either they have 
been killed and have disintegrated, or they have passed back into the 
circulation. As far as size is concerned, heart muscle fibres would be 
adequate for encystment of embryos of trichina. Yet something pre- 
vents this process. The protection is not altogether adequate, however, 
for an interstitial myocarditis results. 

From such diverse examples as these one must conclude that the myo- 
cardium possesses certain protective mechanisms which seemingly are 
inherent. It is reasonable to consider that these mechanisms are of deeper 
biological significance than their immediate objective manifestations, and 
that they have been developed and perpetuated in order to maintain, 
insofar as possible, the integrity of the myocardium, so essential for the 
preservation of the organism. It is probable that many of the varied 
responses of body tissues to extrinsic pathogenic factors have similar sig- 
nificance. Conversely, the recognition of these inherent protective mecha- 
nisms, particularly as they occur in different forms of life, may be of 
use in interpreting generic relationships. 


METACHRONISM IN CILIATED EPITHELIUM 
By LEONARD G. WORLEY 
ZOOLOGICAL LABORATORIES, HARVARD UNIVERSITY 


Communicated January 18, 1933 


Cilia are remarkable not only because they beat, but because they pos- 
sess the ability to act in sequence—metachronously. This metachro- 
nism or coérdinated rhythm results in the production of waves of activity 
which pass over the ciliary field. The direction toward which such propa- 
gated waves travel may be the same as the direction of the effective stroke 
of the individual cilia, or the reverse. In any given case, however, the 
same relation of the effective stroke to the wave of propagation is constantly 
maintained. Thus, in the intestine of lamellibranchs, both are always 
directed away from the stomach and toward the anus, while in the palate 
of Rana pipiens the two are always in direct opposition. In those few 
cases, notably among the Protozoa and some of the Coelenterates, where 
a change in the direction of the effective stroke may occur there appears 
to be a coincident reversal in the direction of movement of the wave propa- 
gation. The cirri of the gills of Nephthys beat at a right angle to the wave 
of propagation and a similar relation exists in the lateral cilia of the gills 
of lamellibranchs. 
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Beating and metachronism are intimately associated and one is rarely 
observed without the other. These activities of the ciliary system are, 
however, not mutually dependent. If a ciliated protozoan or a ciliated 
epithelium is placed in a dilute solution of chloroform, ether or chloretone, 
metachronism soon disappears although the cilia continue to beat for a 
time in an uncoérdinated manner. An anaesthetic disrupts coérdinated 
ciliary movement without simultaneously causing a cessation of beat- 
ing. A comparable result obtains when the region of a ciliated epithelial 
cell occupied by the ciliary cone is cut. When the distal fifth of an iso- 
lated tall columnar ciliated cell from the intestine of Anadonia is cut off 
with a microdissection needle, the cilia lying above no longer beat in 
sequence, but behave independently or synchronously. Beating and 
metachronism are quite definitely distinct elements of ciliary activity. 

Metachronism itself involves two activities which are separable, for 
an epithelium in a weak anaesthetic solution loses its coérdination in at 
least two stages. Metachronism between cells (pluricellular metachro- 
nism) disappears almost immediately when a group of cells is so treated, 
while the cilia of many individual cells continue to beat in a codrdinated 
fashion either metachronously or synchronously. Not until several 
minutes later does unicellular metachronism, the metachronism of indi- 
vidual cells, wholly cease. This effect is most clearly demonstrated if a 
few of the filaments of the gill of Crepidula are subjected to treatment 
with a saturated solution of ether or a 20 per cent solution of magnesium 
sulphate in sea water. It is possible therefore for adjacent cells in an 
epithelium to behave independently. Approximately the opposite has 
been noted in the gill epithelia of Ostrea. Here, under certain conditions, 
all of the cilia of each cell may beat simultaneously, but the cilia of alter- 
nate cells act in sequence. When such a state of affairs appears, the first 
cell of the series beats all of its cilia at once, then, after an appreciable 
‘delay, the next cell acts in exactly the same manner, then the next and so 
on. In this case, pluricellular codérdination is maintained although uni- 
cellular metachronism is lacking. These observations serve to illustrate 
the fact that pluricellular and unicellular metachronism are independent, 
either being capable of expression without the other. 

The basal body is now generally believed to be the motor center for the 
ciliary beat since isolated cilia usually beat if the basal body is present. 
It might be supposed that a neuroid impulse, if such it may be called, 
brings about a metachronic behavior of the cilia of an epithelium by 
stimulating the basal corpuscles in succession. The results presented 
above, together with the observation that isolated individual epithelial 
cells are capable of displaying coérdinated ciliary movement, have led 
me to the conclusion that the individual cells govern their own metachro- 
nism and that some other influence, possibly the neuroid impulse, normally 
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is concerned only with effecting a sequential response of these cells. It 
appears to me that there are two different mechanisms at work here. The 
operation of one of these results in a pluricellular coérdination due to a 
stimulation of the cells in order; the activity of the other, an organiza- 
tion present within each cell, effects a typical response within the cell. 
Coérdinated beating of the cell cilia may be considered the characteristic 
response. 

Another conclusion, based upon the observations just described, should 
be stated. Many authors believe that the metachronism between cells 
is determined within the epithelium and is not dependent upon an extra- 
epithelial influence. If this be the case, this determining factor must 
travel from cell to cell in a horizontal direction with respect to the epi- 
thelium. The synchronous beating of all of the cilia of each cell described 
above for Ostrea, suggests that all of the cytoplasm of a cell may be simul- 
taneously affected. The cell apparently transfuses the impulse with great 
rapidity. Through the epithelium as a whole, however, the progress of 
the impulse is slow. These considerations lend support to the conclu- 
sion that the impulse does not flow through the tissue at a constant rate 
in the manner that the propagation wave crosses the surface. The latter 
probably cannot be considered an external expression of the velocity of 
the impulse through a single cell. This impulse possibly might best be 
thought of as one that progresses in jerks, just as is the case in a true 
nervous impulse that pauses at each synapse. 

Lucas (1931) has suggested that the cell wall acts like a synapse be- 
cause he found that the speed of the propagation wave was influenced by 
the number of cell boundaries per unit length over which it passed. The 
greater the number of cell walls per unit length, the slower was the velocity 
of the wave. Lucas (1932) was, however, unable to observe any delay 
in the progress of the wave across a cell wall. It is my suggestion that the 
delay of the impulse at each cell boundary is of just sufficient length to 
allow it to be overtaken by the wave traversing the surface. The sluggish 
transmission in a ciliated epithelium has never been explained, although 
this is one of the most characteristic features of ciliary activity. Can this 
slow transmission be due to the large number of synapses in the form of cell 
walls that the impulse is required to cross? True nervous impulses progress 
rapidly because of the relatively small number of synapses in their paths. 

The experiments which have led to these conclusions have been carried 
out ynder the kind direction of Professor G. H. Parker. A more detailed 
account of this work will appear in a later publication. 


A. M. Lucas, “‘An Investigation of the Nervous System as a Possible Factor in the 
Regulation of the Ciliary Activity of the Lamellibranch Gill,” Journ. Morph., 51, 
147-194 (1931). 

A. M. Lucas, ‘‘Coérdination of Ciliary Movement. II. The Effect of Temperature 
upon the Ciliary Wave Length,”’ Journ. Morph., 53, 265-276 (19382). 
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THE INTRACELLULAR FIBRE SYSTEMS OF PARAMECIUM 


By LEonarD G. WorRLEY 
ZOOLOGICAL LABORATORIES, HARVARD UNIVERSITY 


Communicated January 18, 1933 


Within the present century, at least two sets of intracellular fibres, 
considered of significance in the codrdination of the beat of the cilia of 
Paramecium, have been described. Rees (1922) has emphasized the im- 
portance in this respect of a group of fibres, each of which has its origin at 
a basal granule in the ectoplasm and then follows a course through the 
endoplasm toward the nuclear region where all of the fibres are said to 
converge and from what has frequently been called a ‘‘motorium.’’ These 
fibres have been termed by some authors the ‘‘fibres of Rees’ and they 
constitute, together with the motorium, the so-called ‘‘neuromotor ap- 
paratus” of Paramecium. This ‘‘neuromotor apparatus’ has been demon- 
strated in sections of Paramecium stained with iron haematoxylen. Rees’s 
figure of this system has been reproduced in several publications. 

An ectoplasmic system of fibres has also been noted in fixed Paramecium 
material. Apparently Schuberg (1905) was the first to observe that the 
basal bodies of the cilia were directly connected by means of longitudinal 
ectoplasmic fibres. These fibrillae, which run the length of the cell, were 
disclosed in specimens killed with osmic acid and treated with dahlia, 
tannin and iron haematoxylen. This investigator pictures the basal cor- 
puscles arranged in longitudinal rows and shows a fibre connecting all 
of the corpuscles of each row. 

Silver impregnation methods have been employed recently by Klein 
(1928; 1929) and Gelei (1932) in revealing the structure of this same system 
of fibres. These authors have described the apparatus in considerably 
greater detail than was done by Schuberg and have called the rather com- 
plicated organization of fibres which they have demonstrated in a large 
number of ciliates, the “‘silver-line system.’’ The longitudinal fibrillae 
have been referred to as the ‘‘direct connecting fibres’ by Klein (1928), 
but according to this writer, there also exist ‘indirect connecting fibres’’ 
in the ectoplasm which connect transversely the longitudinal rows. . The 
direct connecting fibres and the cross fibres divide the ectoplasm into small 
rectangular or hexagonal areas or “‘Feldchen’’ which over the greater por- 
tion of the surface of Paramecium are arranged somewhat like bricks in 
a wall. The paths of the longitudinal fibres are straight while the indirect 
connecting fibres must follow a crooked course around the ‘‘Feldchen.” 
The cross fibres meet the primary connecting fibres in the intervening 
spaces between basal bodies. They are consequently in communication 
with the basal bodies only indirectly. 
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It is noteworthy that these fibre organizations, both the neuromotor 
mechanism of Rees and the silver-line system of Klein and of Gelei, have 
been disclosed only in fixed and stained or impregnated material. As 
far as I am aware, neither sets of fibres have ever been observed in living 
cells. For this reason and others, the existence of such an apparatus as 
Rees has described has been doubted by Alverdes (1923) and Jacobsen 
(1931) and its rdle as a neuromotor device has been questioned by Hofker 
(1928) and Parker (1929). The function of the ‘‘neuromotor apparatus,”’ 
if it exists at all, remains to be ascertained. The silver impregnation 














FIGURE 1 


Longitudinal or direct connecting fibres of the aboral surface of Para- 
mecium showing cilia and their basal bodies and the attachment of the latter 
to the longitudinal fibres. Photograph from a fresh preparation. 


method employed by advocates of the silver-line system might be criti- 
cized because there is no positive assurance that the structures seen in 
such preparations are not, in part at least, of an artificial nature. The 
drying process used by Klein cannot help but result in some distortion of 
the material so treated. Klein (1929) has, however, ascribed a coérdina- 
tion and conductile function to this set of fibres. It occurs to me that 
it would be appropriate to possess unquestionable proof of the existence 
of such fibres before a function is attributed to them. With this considera- 
tion in mind, the following observations are presented. 

It has been found possible by employing microdissection methods to 
remove the ectosarc from Paramecium individuals and to flatten it against 
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the cover glass within a moist chamber. Two microdissection needles, 
with their tips bent at right angles to the main shafts, are usually utilized 
in this process. These are inserted through the ectoplasm and into the 
endoplasm of an individual Paramecium specimen at as nearly as possible 
the same point. They are then moved apart slowly and a large rent out 
of which the endoplasm pours is consequently formed. The surface ten- 
sion of the film of water on the cover slip aids in flattening out a large 
area of ectoplasm in a plane parallel with the glass. With the endoplasm 








FIGURE 2 


Longitudinal or direct connecting fibres (horizontal) and cross or indi- 
rect connecting fibres (vertical) near the cytostome (lower margin) of 
Paramecium, Photograph from a fresh preparation. 


thus disposed of, observations of the ectoplasmic structure not possible 
otherwise are facilitated. 

In living, unstained ectoplasmic sheets secured in the manner just de- 
scribed, the longitudinal or direct connecting fibres of the silver-line system 
can be clearly seen extending antero-posteriorly, joining the basal bodies 
in the manner described by Schuberg and Klein. On the dorsal or aboral 
surface the basal corpuscles of two adjacent rows are usually alternate 
in position or nearly so (Fig. 1). Cilia, basal granules and longitudinal 
fibres are visible. The margin of such a sheet of ectoplasm frequently 
becomes ‘‘frayed out” and single isolated fibres bearing their basal bodies 
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are not uncommonly seen. Such isolated fibres resemble in their appear- 
ance strands of loosely strung beads. 

The longitudinal fibres of the dorsal side appear not to be held together 
by any definite structures other than the undifferentiated ectoplasm and 
they break away from each other and float about in the fluid surrounding 
the preparation. The indirect connecting fibres or cross fibres are sup- 
posed to be of approximately the same diameter and structure as the 
longitudinal connecting fibres and they should be expected to be seen with 
almost equal clarity. They are, however, invisible in my preparations. 
Indirect connecting fibres have never been observed in living material 
in areas where basal corpuscles do not lie opposite each other on adjacent 
rows. In the vicinity of the cytostome (Fig. 2), however, basal granules 
are arranged in horizontal as well as longitudinal rows and in the inter- 
vening spaces between basal granules, cross fibres pass, radiating in several 
directions from the pharynx. These are fairly straight fibres and they 
appear to be continuous across the longitudinal rows. Brown (1930) 
faintly indicates perpendicular fibres of this kind around the cytopharynx 
in the figures accompanying his work. In figure 2, the proximal margin 
of the cytostome is visible at the lower edge of the photograph. The longi- 
tudinal or direct connecting fibres extend horizontally in the picture while 
the cross fibres run vertically, meeting the primary fibres perpendicularly 
or nearly so. Isolated fibres with their attached basal bodies appear 
in the upper right-hand corner of the figure. 

No evidence for the existence of the fibres of Rees has ever been ob- 
tained in my studies, but I believe that it may: be safely stated that the 
direct connecting fibres and the indirect connecting fibres in the vicinity 
of the cytostome in Paramecium are not artifacts. 
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PALAEMONETES 
By F. A. Brown, Jr. 


Communicated January 27, 1933 


tracts the white pigment. 


blood. 








THE CONTROLLING MECHANISM OF CHROMATOPHORES IN 


ZOOLOGICAL LABORATORIES, HARVARD UNIVERSITY 


The common prawn, Palaemonetes vulgaris, possesses four kinds of 
pigment which are active in adapting the animal to its environment. 
These four pigments are red, yellow, white and blue; this report is con- 
cerned with only the first three. The red and yellow pigments are located 
in one type of chromatophore; the white pigment in another. 
chromatophores are far less numerous than those containing the red and 
yellow pigments. Perkins (1928) discovered that injection of a sea-water 
extract of the eyestalks of Palaemonetes would produce contraction of 
their red and yellow chromatophores. Koller (1928) working with the 
sand-shrimp, Crangon vulgaris, confirmed Perkins’ discovery and, in addi- 
tion, located a second hormone from the rostral region of Crangon. 
second hormone caused the chromatophores to expand. Although my 
work has been done principally upon Palaemonetes, several differences 
in the chromatophoral responses of Crangon and Palaemonetes give me 
reason to believe that the mechanisms which control the chromatophores 
in these two forms are unlike. Crangon differs from Palaemonetes in 
having chromatophores containing melanin. As the following results 
show I can confirm Perkins’ statement that the hormone formed in the 
eyestalks of Palaemonetes contracts the red and yellow pigments. 
can add further that there is in this form a second hormone which con- 





The red chromatophores of Palaemonetes are capable of expanding 
or contracting in accordance with the stimuli received through the eyes. 
The completely contracted condition is found in those animals that are 
upon a white background. Complete darkness also causes contraction 
of the pigment. This condition may be produced artificially by injecting 
into an animal a sea-water extract of the eyestalks or of the central nervous 
organs from other Palaemonetes. Removal of one eyestalk quite appre- 
ciably decreases the rate of contraction of the red pigment. 
fore concluded that the contraction is brought about by a substance found 
principally in the eyestalks but in smaller amounts in the central nervous 
organs. This material is carried through the body of the animal by the 


The fully expanded condition of the red chromatophores is seen in ani- 
mals that are upon a black background. Experimentally, this is induced 
by removal of the eyestalks from the animal, or by anesthesia. Stoppage 
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of the circulation in a region of the body results in an expansion of the red 
pigment in the part thus deprived of flow. If such a drastic method as 
isolating a piece of the integument with its contained chromatophores 
in sea-water is used, the red pigment expands at a normal rate for the 
first half to three-quarters of an hour and then the rate becomes lower 
than that seen in normal black adaptation. This decrease in rate is due 
to a change in the chromatophore itself rather than to an antagonizing 
hormone, for if the chromatophore is more than two-thirds expanded 
when isolated, complete expansion will follow. All experiments of inject- 
ing sea-water extracts of various body parts, occluding arteries supplying 
different regions, or destroying parts of the body, have thus far revealed 
no evidence in favor of a red-expanding substance that is introduced into 
the blood stream at any time. Hence it is concluded that no specific 
substance outside of the chromatophore is responsible for producing ex- 
pansion. 

The white reflecting chromatophores undergo the same degrees of ex- 
pansion or contraction as do the red ones. They are fully contracted in 
animals upon a blue or black background. Darkness also causes contrac- 
tion and injection of sea-water extract of eyestalks brings about a con- 
traction. Thermal or electrical stimulation of the eye-stubs in animals 
without eyestalks causes contraction of the white pigment. These facts 
indicate that contraction of the white chromatophores is accomplished by 
a hormone that originates in part in the eyestalks and in part in some 
other region of the body. 

White chromatophores most widely expanded are seen in animals upon 
a white background in the light. Much the same condition is seen upon 
a yellow background. Occasionally white chromatophores are found 
contracted on a white background. Also, these chromatophores may be 
found in any state of expansion or contraction in animals without eye- 
stalks, though the response to background is no longer displayed. The 
white chromatophores are less predictable in their behavior than are 
the red ones, though the vast majority of animals exhibit the conditions 
described in this and the preceding paragraph. Because of the great 
similarity between the expanding white and the expanding red chromato- 
phores, both in rates of change and in the experimental procedures that 
bring about the change, it seems that here also no specific expanding 
hormone is present. 

Thus far nothing has been said of the ever present yellow pigment which 
is located in the same chromatophore as the red. Without exception 
this pigment responds in the same way as does the red, the only difference 
being that the yellow expands more rapidly and contracts more slowly 
than the red. There appears to be no reason to seek a controlling sub- 
stance for it separate from that which controls the red. 
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The behavior of the red and of the white chromatophores is such that 
their contractions cannot be regarded as due to the same substance. Upon 
a black background the white is contracted and the red fully expanded. 
Upon a white background the condition is reversed, the white being fully 
expanded and the red contracted. The possibility that the hormone that 
contracts the white chromatophores causes expansion of the red, is elimi- 
nated when it is seen that darkness, a blue background or injection of eye- 
stalk extract, produces a contraction of the red and of the white chromato- 
phores simultaneously. Similar evidence shows that the substance that 
brings about contraction of the red chromatophores cannot be necessary 
for complete expansion of the white ones. Animals upon a red background, 
or animals without eyestalks, often show both types of chromatophores 
expanded at the same time. 

These facts require for their explanation the assumption of at least 
two different substances carried in the blood to affect the chromatophores. 
One substance contracts the red chromatophores and the other contracts 
the white ones. Expansion in both cases is not dependent upon a hor- 
mone but is a property inherent in the chromatophore itself. 


Koller, G., ‘‘Versuche tiber die inkretorischen Vorgange beim Garneelen-farbwechsel,”’ 
Zeit. vergl. Physiol., 8, 601-612 (1928). 

Perkins, E. B., ‘‘Color Changes in Crustaceans, Especially in Palaemonetes,” Jour. 
Exp. Zoél., 50, 71-105 (1928). 


DEGENERATION OF MELANOPHORES IN FUNDULUS 
By J. M. ODIORNE 
ZOOLOGICAL LABORATORIES, HARVARD UNIVERSITY 


Communicated January 21, 1933 


The killifish, Fundulus heteroclitus, like many other fishes, assumes a 
light shade in light surroundings, and a dark one in dark surroundings. 
The changes in the melanophores which are chiefly instrumental in bring- 
ing this about have been the subject of many investigations. 

During the course of such an investigation it was noticed that fishes 
kept in light surroundings over long periods of time became unable to adapt 
themselves to dark surroundings to the usual degree. Closer examina- 
tion showed that this inability was the result of degeneration of some of 
the melanophores, followed by loss of the pigment contained in them. 
The experiments which have yielded the results briefly set forth here 
were undertaken at the suggestion of Professor G. H. Parker and carried 
on under his supervision. 
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Using a binocular microscope equipped with an ocular micrometer 
ruled in squares, the melanophores in a particular region of a particular 
scale on the back of a given fish were counted. The area of this region 
was approximately 0.75 sq. mm. and it was found possible to locate the 
same region again without serious error. Thus the melanophores in such 
a region were counted at intervals and any change noted. 

The melanophores of fishes that are placed in white surroundings soon 
contract and remain so. After a few days, examination shows that the 
melanophores are degenerating, for clumps of pigment granules are visi- 
ble, and counts reveal a decrease in the number of melanophores. As 
this process goes on the pigment granules form large aggregations which 
accumulate at the edges of the scales, so that the boundaries are empha- 
sized. The pigment granules are later lost at the surface of the body. 

The following table shows the decrease in the number of melanophores 
seen in fishes kept on a white background before a window and at room 
temperature. The counts for 23 fishes are here added together. 


Total number of melanophores at beginning of 


experiment 1835 
After 7 days 1714 
After 11 days 1621 
After 18 days 1478 
After 25 days 1283 
After 28 days 1224 


If black surroundings are substituted such a reduction fails to appear, 
though all other conditions remain the same. Instead the number of 
melanophores remains nearly constant. 

If fishes which have been for some time in white surroundings are trans- 
ferred to black surroundings the decrease is checked, and a gradual in- 
crease begins. 

At a temperature of 10°C. fishes kept in white surroundings show only 
a slight reduction in the number of melanophores over a period of several 
weeks though the melanophores contract and the fishes assume a shade 
similar to that of specimens kept at room temperature. But if at the end 
of the period the temperature is allowed to rise to room temperature, de- 
generation quickly sets in. Experiments are now going on which it is 
hoped will yield further information on this matter. 

Ordinarily the experimental fishes have not been fed, but experiments 
in which controls are fed indicate that inanition is not responsible for the 
degeneration. In winter Fundulus lives for months in the laboratory 
without food and it does not seem likely that lack of food would cause 
degeneration of melanophores within a week when the animals are freshly 
caught. Furthermore, if inanition were the principal factor it is difficult 
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to understand why the melanophores do not degenerate when the fishes 
are kept in black surroundings. Murisier (1920) found that when no 
food was supplied to trout after hatching, the melanophores degenerated 
provided the fishes were kept on a white background. He attributed 
this degeneration to starvation, although if the young trout were kept 
on a black background the melanophores actually increased in number. 

It seems more probable that the difference in the degree of expansion 
of the melanophores, occasioned by the different backgrounds, are respon- 
sible for the results. Gradually an hypothesis has taken form which pos- 
tulates that continued expansion of the pigment cells favors an increase 
in their number and in the amount of pigment contained in them, while 
continued contraction opposes such an increase. The development of 
this hypothesis is connected with the names of Keeble and Gamble (1904), 
Franz (1910) and Babak (1913), but a discussion of their work is out- 
side the scope of this paper. Other experiments of Murisier seem to sup- 
port such a theory, as do experiments which I have performed on young 
paradise-fishes (Macropodus opercularis). The degeneration of melano- 
phores in Fundulus described in this paper leads me to add that continued 
contraction not only checks any increase in the number of melanophores 
but is likely to bring about a decrease in their number. It would appear 
that the experiments of Murisier cited above may be interpreted in such 
a way as to support this idea. 

Kuntz (1915) mentioned a decrease in the number of melanophores in 
Paralichthys albiguttus when kept on a white background, followed by 
increase when kept on a black background. Within eleven days a reduc- 
tion of about 30 per cent occurred, but this seemed to be the limit. In 
Fundulus a longer time is required to produce such a decrease, but a re- 
duction of 50 per cent is common, and in occasional cases it is much greater. 
Kuntz attributed his results to the influence of the background, but made 
no attempt at any hypothesis. 

Recently I have learned of a similar investigation which has been con- 
ducted by Dr. F. B. Sumner and Mr. N. A. Wells. These investigators 
have kindly sent me a copy of the manuscript of a paper which is to appear 
in an early issue of the Journal of Experimental Zodlogy. The guppy, 
Lebistes reticulatus, has been used in their experiments. This fish shows 
reduction in the amount of melanin and in the number of melanophores 
when under the influence of a white background, while under the influence 
of a black background the melanin increases in amount and the melano- 
phores in number. My results are in substantial agreement with theirs. 
In Fundulus the process of ejection of melanin, after degeneration of the 
melanophores in which it was contained, is similar to that which they 
describe for Lebistes, and Fundulus also possesses ‘“‘iridosomes,”” which 
frequently become very conspicuous in fishes whose melanophores have 
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been long contracted or are degenerating. It is not surprising that the 
results obtained from these fishes should be in agreement, for Fundulus 
and Lebistes are related forms, but it is interesting that they have been 
obtained simultaneously and independently, though with different methods 
of attack. 
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GROUPS IN WHICH EITHER ALL THE OPERATORS OR ALL THE 
SUBGROUPS OF THE SAME ORDER ARE CONJUGATE 


By G. A. MILLER 


DEPARTMENT OF MATHEMATICS, UNIVERSITY OF ILLINOIS 


Communicated February 10, 1933 


It is obvious that a necessary and sufficient condition that all the sub- 
groups of the same order in an abelian group are conjugate thereunder 
for every such order is that the group is cyclic but for non-abelian groups 
it is somewhat difficult to find a necessary and sufficient general condition 
that all the subgroups of the same order contained therein are conjugate 
under the group for every such order. One of the main objects of the 
present paper is to prove that whenever all of the Sylow subgroups of a 
given group G are cyclic then all of its subgroups of the same order must 
be conjugate thereunder for every such order. While this condition is 
sufficient for the conjugacy of all the subgroups of the same order con- 
tained in G it is not always necessary as results from the fact that all the 
subgroups of the same order contained in the icosahedral group are con- 
jugate under it. 

Suppose that all the subgroups of the same order contained in G are 
conjugate under it for every such order and that G involves a subgroup 
H of order p”, » being an odd prime number. If H were non-cyclic and 
involved operators of order p? then G would involve both a cyclic and a 
non-cyclic subgroup of order p? since every non-cyclic group of order 
p” contains a non-cyclic invariant subgroup of order p?. Hence it results 
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that every subgroup of order p” contained in G is either cyclic or involves 
only operators of order p. If G involves a subgroup of order 2” which is 
not cyclic but involves operators of order 4 then this subgroup must be 
the quaternion group and it must also be a Sylow subgroup of G since this 
subgroup of order 2” can involve only one subgroup of order 2 and a Sylow 
subgroup of G must have the same property. The group of order 120 
which has a (2, 1) isomorphism with the icosahedral group but does not 
involve this group is an instance of such a G. Hence it results that if 
all the subgroups of the same order contained in a group are conjugate there- 
under for every such order then every prime power subgroup must satisfy 
one of the following conditions: it is cyclic, involves only operators of prime 
order or it is the quaternion group. 

Suppose now that every Sylow subgroup of G is cyclic and that the 
order of G is g = py"py”...p%, where pi, po, ..., p, are distinct prime 
numbers in descending order of magnitude. It is well-known that G is 
solvable and that it contains one and only one subgroup of order pj'd%’. . . 
pn, where a = 1, 2, ...,¥% — 1. In the special case when \ = 2 it is 
easy to see that all the subgroups of order p7'p}’, m1 = m, mz = me, contained 
in G must be conjugate under G since G involves only one subgroup of 
order p{' and hence all of these subgroups involve the same subgroup of 
order p7' and these subgroups of order p}’ must be conjugate under G 
since its Sylow subgroup of order p$” is cyclic and hence all of its sub- 
groups whose common order is an arbitrary divisor of p}" are conjugate 
under it. 

Mathematical induction may now be employed to prove that all the 
subgroups of the same order contained in G are conjugate under G for 
every such order since it may be assumed that this is true for the subgroups 
whose orders involve all the prime factors except the smallest which divides 
the order of such a subgroup. Each of these subgroups is therefore in- 
variant under a cyclic group whose order is a power of this smallest factor 
equal to the highest power thereof which divides the order of the subgroup 
in question. Hence the following theorem has been established: Jf all 
the Sylow subgroups of a given group are cyclic then all of its subgroups of 
the same order are conjugate under 1t for every such order. 

In what follows G will be used to represent a group such that all the 
operators of the same order contained in G are conjugate under it for every 
such order. If such a G involves a subgroup of index 2 and of odd order 
none of its operators can be transformed into its inverse under this sub- 
group. Hence such a G contains an operator of order 2 which is commu- 
tative only with its powers under G. It therefore results that all the 
operators of G which are not found in this subgroup must be of order 2 
and hence this operator of order 2 transforms all the operators of the given 
subgroup into their iriverses. This implies that this subgroup is abelian 
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and hence of order 3 since otherwise its operators of the same order could 
not all be conjugate under G for every such order. 

If G would involve a subgroup of index 2 and of even order all of its 
operators whose orders are powers of 2 and which do not appear in this 
subgroup would be of a larger order than any such operator appearing in 
this subgroup since otherwise all of these operators could not appear 
in a single set of conjugates under G. Such an operator would have to be 
transformed into its inverse by an operator appearing in this subgroup. 
Hence we have arrived at a contradiction by assuming that G contains 
a subgroup of index 2 and of even order and the following theorem has 
been established by the arguments of this and of the preceding paragraph: 
The symmetric groups of degrees 2 and 3 are the only ones which have the 
properties that each of them contains a subgroup of index 2 and that all of 
their operators of the same order are conjugate thereunder for every such order. 

If a group contains a cyclic Sylow subgroup of composite order then its 
operators of this order cannot all appear in the same set of conjugates 
under the group since a generator of such a group is transformed into an- 
other generator thereof by an operator whose order divides the order of 
this cyclic group. In particular, none of the Sylow subgroups of G can 
be cyclic and of composite order. It is also obvious that G cannot involve 
an invariant subgroup of odd index since the quotient group corresponding 
thereto would be of odd order and hence some operators of this quotient 
group could not correspond to their inverses under G, which is contrary 
to the hypothesis. In particular, if a group contains an invariant Sylow 
subgroup of order 2” then it must involve operators of the same order 
which are not conjugate under it, and if it is solvable it contains such 
operators except in the two cases noted at the close of the preceding 
paragraph. 

Suppose that G is a transitive group of degree p, p being a prime number. 
When ~ > 3 it is obvious that G cannot be either metacyclic or a subgroup 
of the metacyclic group. Hence it must be at least doubly transitive 
and therefore of even order. Its operators of order p generate a simple 
group which must also be doubly transitive. As this does not involve 
any cyclic permutations of order p — 1 its operators of order p cannot all 
be conjugate under it. Since the corresponding quotient group is known 
to be cyclic! it results that G would involve an invariant subgroup whose 
index would be either 2 or an odd number. This is contrary to the results 
obtained in the preceding paragraph and hence it has been proved that 
the only two transitive groups of prime degree in which all the operators of 
the same order are conjugate for every such order are the symmetric groups of 
degrees 2 and 3, respectively. These are clearly also the only two non- 
perfect groups which have this property. 

From what precedes it results that if there exists a group G besides the 
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two just noted which has the property that all of its operators of the same 
order are conjugate under it then it must be perfect and have an order 
which is divisible by 4. Moreover, it cannot involve a subgroup of prime 
index since it could otherwise be represented as a transitive group of prime 
degree, which has been proved to be impossible. If such a perfect group 
were not simple and of composite order its quotient group would have 
this property. Hence it results that the group of lowest order which 
could have the property in question would be simple and of composite 
order. If it contained an operator of order p”, p being a prime number, 
the generators of the cyclic subgroup generated by this operator would 
be transformed under the group according to a group of order p”~?. In 
particular, if a Sylow subgroup of such a group were abelian then all of 
its operators would be of order ». No such simple group seems to be 
known. 
1G. A. Miller, Bull. Amer. Math. Soc., 4, 141 (1896). 


AN EXTREMELY SIMPLE METHOD OF PERIODOGRAM 
ANALYSIS 


By DINSMORE ALTER 
DEPARTMENT OF PHYSICS AND ASTRONOMY, UNIVERSITY OF KANSAS 


Communicated January 23, 1933 


Despite the many advantages of the periodogram, when properly 
interpreted, over all other forms of analysis, in the search for hidden cycles 
or periodicities, it has never been popular on account of the tremendous 
amount of computation involved. 

There will be issued shortly a study of all forms of the periodogram 
which have been used and a critical discussion of the advantages and 
disadvantages of each, by Doctor E. S. Pearson and the writer. Among 
these forms the correlation periodogram seems to be the best when 
thoroughness of analysis throughout a wide range of types of data and 
ease of calculation are considered. 

In this article is presented a new form of periodogram as powerful as 
the correlation form and so simple that it requires little more than onie- 
third as long to compute. In addition to this, all calculations are made 
on an ordinary adding machine. The accompanying rather long example 
was computed in twelve hours. It is hoped that with this great saving 
in time many investigators will make use of the periodogram method. 
In the work of astronomers and meteorologists this is especially important. 
Given n data x, %2....x, with x; representing any one of the data. 
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We shall consider the values given as departures from the mean, so that 
=x; = 0. Denote by / the time interval between two data as, for example, 
between x, and x,,;. Square brackets will be used to denote absolute 
values. Form all the sums x; + x;,; and then sum the absolute values 
of all these sums for a given fixed value of /. 


Define: 


n—l 
M, = pm [2; + X41] 


i=1 n—l 


If after the time / there is a partial repetition in the data, M; is increased. 
Plot the M/;, as ordinates and the various trial values of / as abscissae. 
The resulting curve has maxima for those values of / for which there is a 
partial repetition in the data. 

In order to investigate the variations of MM, and their usefulness in 
determining hidden cycles or periodicities, we shall compute the relation- 
ship between M, and the correlation coefficient. 








Define 
n 2 n—l 2 n 2 
xi xy Xi 
a i v oe ae 7 
ox =), — 50x =D) Sage =e). 

i= 1 1 j=, n —1 isp n—l 

n—l Xj X41 

1% 

7, = (As yy Rie 


i=, (n — l)o\x-ox 
Using Peter’s equation, thus assuming normal distribution of the x’s: 
€4 = 0.8453 M, 


and using the standard equation 








yA2 
& = 0.6745 \=4 j = (0.6745 V ox a o3x + 201X'02x'7;. 


Equating these 
0.6745 


M, = 0.8453 Vox + 02x + 2oyx-o2"-7. 





Our population may be considered as the m data so that ox is a fixed 
value. ox and oox are, therefore, approximations to ox from samples of 
size m — 1. Unless / is considerably larger than '/2n they will each be 
rather close approximations to ox. 


Write 


Ox =oxt+A 
Define: 
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2 2 2 
o;x + ox A 
i laces aie = ox + Aox + re 


OyX'O2X = ox + Ao\x. 


Neglecting A? we then write 


__0.6745 
M, = ‘* 
M, V9 oof Vi+r. 








Differentiating the logarithms we obtain 


5M, b0'x 1 &, 


M, o'x ee a Ue 


At / = 0 and/ = n/2, o’x takes the forced value ox and éc’x = 0 for 
these values. If 1] < m/2, 50’x is very small since o’x differs from ox due 
only to inclusion of some of the x’s twice. The largest value of d0’x, 
within this region, is to be expected at / = 1/,m where half the data are 
used twice. 

For values of / > !1/2n, some of the data are omitted from o’x; i.e., those 
between x; and x,_;. Therefore, in this region 60’ can become large as 
l approaches n. 

In the region where 50’x/o’x is small, we have the ratio of accidental 
errors of M; to M, equaling one-half the ratio of accidental errors of 
r,to (1+ 7). Therefore, in this region the new periodogram will be less 
affected by accidental errors than is the correlation periodogram. The 
reverse will be true as / approaches n. 

The last statement is easily shown by an extreme case m = 10,000, 
1 = 9500. 17, is given quite accurately by 500 pairs of data. But suppose 
that the first 500 data and the last 500 are subject to large oscillations with 
respect to those of the intermediate data. oix and ox will each be far 
larger than ox and, therefore 50’x will be very large. As an example of 
such data we may take the rainfall of the eastern United States, which 
had large oscillations in annual values between 1820 and 1827 and oscilla- 
tions larger than the average from 1912 to 1919. 

Few periodograms are computed to values of / much larger than one- 
half n. 


In the accompanying example are shown superposed periodograms of 
R Scuti. The characteristic feature of these periodograms is the very 
definite beat pattern given by interference of large, nearly equal periodicities 
in its luminosity. m = 228. The time unit for lags is 25 days and the 
calculations have been pushed to the extreme value of / = 168. The 
two forms of periodogram follow each other with surprising exactitude. 
For lags less than about 3000 days the new periodogram seems to be the 
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more symmetrical in form. For greater lags the correlation form is the 
more symmetrical. 
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Having computed the periodogram, it may be best in investigating the 
accidental frequency distribution of the highest peaks, to compute 7; and 
consider the ratio 7;/or; for these exceptional values of /. This procedure 
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is necessary due to the weakness that only the approximate frequency 
distribution m,; is known except at the value / = !/2 n. 

If the data have not been reduced to a zero mean, the shape of the 
periodogram is not affected but the whole is raised by twice the mean 
value. The modification of the probability formulae for such data is a 
routine task. 


SOME REMARKS CONCERNING SCHRODINGER'S WAVE 
EQUATION 


By GrorGcE D. BIRKHOFF 
DEPARTMENT OF MATHEMATICS, HARVARD UNIVERSITY 


Communicated January 30, 1933 


In the present note I propose to approach the wave equation of Schré- 
dinger by a method which, although closely related to methods used by 
de Broglie and Brillouin,! Schrédinger? and Dirac,’ is distinct from these 
and has the advantage of fixing the position of the wave equation from a 
purely mathematical point of view: namely, the wave equation is the 
simplest form of linear partial differential equation involving a parameter, 
for which the classical process determining asymptotic series solutions 
gives a ‘‘multiplier equation” identical with the Hamilton-Jacobi equation, 
while the Cauchy characteristics of the multiplier equation (along which 
the asymptotic ‘‘wave packet’ is always propagated) then become the 
dynamical trajectories of the corresponding Hamiltonian system. 

Let L(y, \) = 0 be any linear (ordinary or partial) differential equation 
in y with n 2 1 independent variables x, ... x,, and involving a parame- 
ter \. This parameter is to be thought of as large in absolute value. 
In the case which interests us certain of the coefficients of y and its de- 
rivatives in the above equation become large as || increases. Under 
these circumstances it has been found in many cases that the operation 
of differentiation of y or of its derivatives with respect to any variable 
x; is asymptotically equivalent to multiplication by \0x;/0S, where S 
is a definite function of x, ... x,. Thence one is led to an asymptotic 
series solution, 


y~ (oy +24 (1) 


where v, v1, ... do not involve 4. It is with certain general facts con- 
cerning this classical, formal process and its relation to the Schrédinger 
wave equation that the present note deals. If we introduce the modified 
differential operators 
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the original equation in y may be written as 


[1] (1] [2] 
L(y, 5, ...2 2, is 4 oe ee »)=0 (2) 
Ox OX, Ox?” 





where LZ may be expanded in the form of an infinite series in 1/), 
L 
L=lbt+—+... (2’) 


in the case which we are considering. Here Ly may be called the “principal 
part” of Zo, and may be more explicitly written as follows 


[2] 
Loy) = w+ ree t+ ety id 


; 3 
Ox; oe i (8) 


where we assume 
fi; = Shin Sie = Sing = jet = Shik = Seg = Enji, ete. 


If now the formal series (1) be substituted in the equation (2), and the 


coefficients of e*/X* on the left be equated to 0 for k = 0, 1, ..., the first 
equation determines S and may be written 

P(m,...4 - sees 28) m0, (4) 

” Ox,’ OX» 
where 
oS oS 
= + iat i#—-— +. 5 
=f res 2 bis a (5) 


is a polynomial of degree m in OS/0x, ... OS/0x, if m is the order of the 
equation in y. Conversely any polynomial P determines a corresponding 
Ly. We shall term (4) the “multiplier equation” for L(y) = 0. This is 
a differential equation of the first order in S which does not contain S and 
is in general non-linear. 
When we proceed with the further similar equations for k = 1, 2, 

we find that the second equation determines v in (1), the third determines 
%, etc. We shall be interested primarily in the equation for v. Let us 
observe first that we have, to terms in 1/4, 

 . os (S., 1 (2. ~ ) 

Oar mn” ten ee 


os 1 oS 
. *([8 432 |s a et =m) 
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a ~ & “(2542 2 [2 +5 2 | +228 25Sn), te 
Ox; Ox; Ox; NOx;JLOx; Ox; d Ox; Ox; 
Here we have used an obvious operational notation. On substituting 


these expressions in equation (2), the terms in 1/\ which involve y, dis- 
appear because of the multiplier equation (4), and there remains the 


following equation for k = 








OP ov or 8S 
—+- a( )es = 0, (6) 
i dy; Ox; 4,7 Oy; Oy; Ox; Ox; 
where we have written y; = OS/0x; fori = 1, 
Suppose now that we take any complete solution S(m, ... %,; 4, 
. Cy-1) + C, of the multiplier equation. Then x, ... Xp, 1, ... Vn 


may be regarded as independent except for the equation of condition 
P = 0. If we write 





ee. Te 
—=-=— ¢=1,...8n 7 
. ( ) (7) 
and let denote the coefficient of v) in (6), the equation (6) becomes 
a + dv) = 0, whence v = xe /™. (6’) 
Hence the equation (6) specifies how v varies along any curve x; = x;(r), 
(¢ = 1,... ), in m, ... %_ Space representing a solution of (7). 
But we have also for any 7 along such a curve 
dy; OS day _ eS oF 





dr F' Ox;Ou;dr F Ox; Ox; dy; 
Furthermore, from the ee equation (4) we infer (in x, ... x, space) 
oP oS 


ee = 0, 
= = ae Oy; Ox; OX; 





so that we deduce 
YS 


= a (= 1,... 2). (8) 


The equations (7), (8) are the 2m ordinary differential equations in 
X1, ... Xqy Vy --- Yq Of Hamiltonian type which define the Cauchy char- 
acteristics of the multiplier equation (4), provided we restrict attention to 
those solutions (involving 2 — 1 arbitrary constants) for which P = 0. 

The general significance of the above results is clear: Consider a regular 
transversal surface T of m — 1 dimensions so that, for the given function 
S, the corresponding characteristics intersect T in one and the same sense 
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throughout. The value of 1 may then be given arbitrarily on the trans- 
versal surface 7, and is then determined elsewhere by integrating along 
these Cauchy characteristics. 

For k # 0, the equations analogous to (6’) are of the form 


At + dy + = 0 (= 1, 3,3; ...), (6") 


where ®, depends linearly and homogeneously upon v%, 1%, ... v,~1. Hence 
each coefficient v, is determined up to a constant multiple of v by its 
value on T. 

If in particular v, v1, ... are taken as 0 in TJ except throughout some 
region V of 7, and as vanishing along the boundary of this region, to- 
gether with their partial derivatives, these functions may be taken to be 
0 except within the corresponding tube of characteristics (on which they 
are determined to the extent indicated above), and vanishing similarly 
on the boundary of the tube. Such a region V changes position as the 
“time’’ + changes, and the series solution thus gives rise to what we shall 
term an “asymptotic wave packet.” 

Let us call any function v which satisfies the linear partial differential 
equation (6) an “amplitude function,’ associated with the corresponding 
“phase function” S. 

If v9 is any assigned amplitude function of S whatsoever, then the integral 


[= rd vedx x2. . .dX, (9) 
V(r) 


does not vary with r. Here V(r) is the region into which the arbitrary region 
V at r = O%s carried along the Cauchy characteristic in ‘time’ r. Conversely 
if this integral is independent of 7 for every V, v9 will be an amplitude function 
of S. 

To prove this result we observe that the integral J may be written 


fae os Xn) Jdx.. dX, 
- 


where the transformation from V to V/r is x; = f(x, ... x»), (@=1,... m), 
and J is the functional determinant 0(x, ... x,)/O(x1, ... %,). Hence 
the derivative of J at 7 = 0 (when x; = x;,) is 


Ov OP 4 
2-— — + m— |mdx...dxy. 

Tie Ox; Oy; "de tiie is 
But for Ar small we have 


lea +. aoe 5 hy Aye 
Ov; 


and, in consequence, 
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2 2 2 
r=14(5 OP 4 p_OP oS) art. 











+ Ox. OV; 47 Oy} Oy; Ox; Ox; 
whence 
dJ 


o'P Or oS 
dr sai » i + z; . 


T=0 i Ox; oy; i,j oy; Oy; Ox; Ox; 














Hence the integral written above vanishes because of the condition 
(6) upon v%. Consequently dJ/dr = 0 for + = 0, and likewise for any 7, 
so that J is necessarily constant. Conversely if the integral J is inde- 
pendent of 7 for all regions V, it is clear that v) must be an amplitude 
function. This completes our proof of the above italicized statement. 





Suppose now that we take m + 1 variables /, 1, ... x,, with 
os os os 
P zoe H “ee n> e..2:¢ & 
eg (= *m x, vg 


so that the ‘multiplier equation’ takes the form of the Hamilton-Jacobi 
partial differential equation with energy H. The corresponding “principal 
linear equation’ L)(y) = 0 is then the usual Schrédinger wave equation 
2ri OW 2mt O 2ni O 
ome we bt EN hy, Bg i = 0 10 
ie ag (« ag h on) ¥ (20) 





provided we take } = 2mi/h. Furthermore the corresponding Hamil- 


tonian equations are 
os 
dt ot 
1 





— dr 
together with 
a oe (i n) (11) 


so that 0S/0t = const., and r = ¢ + const. 

The Schrédinger wave equation is therefore merely the principal equation 
which has the usual Hamilton-Jacobi partial differential equation H as its 
multiplier equation. 

Evidently ¢ = f is a transversal surface T in this special case and the 
n + 1 dimensional element of volume is didx,...dx,. Furthermore we 
may write 


S = 65 + Sm, ... Sa; Gr ++» Cent) Goan 


since y, = 0S/dt appears only linearly. Here the Cauchy characteristics 
become the dynamical trajectories of the corresponding Hamiltonian 
system (11). 
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If we write to a first approximation 


y poe eu, y as ee, 

(v real, \ pure imaginary), we see that along any part of an asymptotic 
wave packet, the integral /Yydidx,...dx, reduces to fu*didx,...dx, 
which by our general result remains constant. It follows then that 
i vydx. . .dx, also remains constant. To establish this we begin with 
the m + 1 dimensional element of volume, taking limits 4 and #4, for ¢ with 
t; — t) = At, and then let At tend to zero. Hence we arrive at the follow- 
ing result: 

In the special case of the Schrédinger wave equation in y, the “asymptotic 
wave packets” follow the corresponding dynamical trajectories, while the 
squared amplitude integral [ydx;...dx, remains constant over any part 
of the packet, at least to the order of terms in 1/i.° 

In conclusion, it is desirable to note what happens under any change 
of independent variables in the general case 


mm =film,...x,) (= 1,... 0), 
It is to be observed first of all that because of the identities 
Olly i” dx; 0! oly Ox, Ox, O?!y 1 O%x, tly 


ae =) = — rae "9 t ” 
Ox; j Ox; Ox; Ox; Ox; kil Ox; Ox; Ox, Ox) r kb Ox; Ox; Ox, ss 














the components of 


1 
L=bt+yhit.. 


do not remain individually invariant in the equation (1), and in particular 
the principle part Lo is not carried over into the new principal part by 
the ordinary rules. In fact the coefficients £ in the principal part transform 
by the rules valid for the attached Hamilton-Jacobi equation. Hence 
Schrédinger’s wave equation in the form (10) is only maintained (in 
general) under a linear transformation of the independent variables. 

This fact indicates that any coordinate system from which we start is to 
be regarded as a privileged absolute system of reference for the Schrédinger 
wave equation, up to an arbitrary linear transformation. 


1 L. de Broglie and L. Brillouin, Selected Papers in Wave Mechanics, 19-54 (translation, 
London, 1928). 

2 E. Schrédinger, Collected Papers on Wave Mechanics, 1-30 (translation, London, 
1928). 

*P. L. M. Dirac, Quantum Mechanics, 119-123 (Oxford, 1930). 

4 See my paper, Trans. Am. Math. Soc., 9, 219-231 (1908). 

5 Dirac asserts incorrectly (loc. cit., p. 121) that the amplitude is constant along a 
‘dynamical trajectory. He overlooks the terms in the second partial derivatives in 
the equation corresponding to (6) in the special case. 
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ON THE MEASURE OF THE NON-SPECIAL GEODESICS ON A 
SURFACE OF CONSTANT NEGATIVE CURVATURE 


By Gustav A. HEDLUND 
DEPARTMENT OF MATHEMATICS, BRYN MAwr COLLEGE 


Communicated February 7, 1933 


Let E denote the unit circle in the complex z-plane and let W be its 
interior. Let the metric be defined as 


ds* = Ads 

(1 — 22)? 
A hyperbolic geometry in W is thereby defined, the straight lines of which 
are segments of circles orthogonal to E. These will be called NE (non- 
euclidean) straight lines, and that part of a NE straight line with one end- 
point in W and the other end-point on E will be called a NE ray. « 

As is well known, there exists within V a regular non-euclidean polygon 
S of 4p sides, p > 1, with sides and vertices all equidistant from the origin 
0, and with angle 7/2p at each vertex. Associated with this is a group T 
of linear fractional transformations which carry E into itself and such 
that under I, S is carried into a set of regions which fill up Y without 
lacunae. Considering congruent points as identical there is thereby 
defined a closed orientable surface of genus » and of constant negative 
curvature. This surface will be denoted by S, as only points of the 
polygon S need be considered, provided the sides of S taken in certain 
pairs be considered identical. 

The direction of a NE ray emanating from a point P(x, y) of S, will 
be measured by the angle ¢ formed with the NE ray from P to the point 
z = 1, ¢ to be taken in the counterclockwise sense and such that 0 S ¢ S 2z. 
Let M denote the three dimensional region (x, y, ¢) obtained by associating 
with each point (x, y) of S, the values of ¢,0 S ¢ S 2x. If the point 
(x, y, 0) is considered identical with (x, y, 27) the region M is a closed 
manifold. 

A set of points, F, of M, will be said to approximate all points of M with 
degree of approximation e > 0, if (x, y, ¢) being any point of M, there is a 
point (x’, y’, ¢’) of F such that 


V (x! — x)? + (y’ — 9)? + @' — 6)? <e 





A directed NE ray evidently determines a set of points, G, of M. A 
directed NE ray will be said to be non-special if the set of points, G, ap- 
proximates all points of M with degree of approximation e > 0, ¢ arbitrarily 
small. If a directed NE ray is non-special, the same ray, but oppositely 
directed, is also non-special. Ifa NE ray is not non-special, it is special. 
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THEOREM I. Given any point P(x, y) of S, there are non-special NE 
rays, (x, y, >), with this point as initial point. Moreover, the set of values 
of » determining such non-special rays is everywhere dense in the interval 
05 ¢ 8 2n. 

Let those points (x, y, ¢) of M which determine non-special NE rays 
be called non-special points. From a conjecture of Birkhoff! it would 
appear that the set of non-special points of M should be a point set of the 
measure of M. The object of this paper is, first of all, to outline a proof 
that this is the case, then to extend the results to complete geodesics 
which are made up of pairs of non-special rays. 

Considerable use is made of the work of J. Nielsen? and it will be neces- 
sary to assume his notation and results as given in the paper cited. With 
each point of E, except for a denumerable set D, Nielsen associates a 
unique infinite symbol. With each point of D there is associated two 
infinite symbols, which may be considered equivalent. There are two 
processes which give the infinite symbol of a point Q of E, only one of 
which needs to be considered in this outline. First E is divided into a 
set of 4p intervals, each of which is given a symbol. These constitute the 
intervals of the first stage. Then these intervals are subdivided to obtain 
the intervals of the second stage, each of which has a symbol. By a 
continuation of this process E is divided into more and more subintervals 
and the point Q becomes a common point of an infinite sequence of in- 
tervals belonging to successive stages. The symbol associated with Q is 
the sequence of single symbols of these intervals containing Q. 

Let those points of E for which the corresponding symbol contains every 
finite admissable sequence of single symbols be called non-special points 
of E. Let the other points of E be special. 

THEOREM II. The NE ray PQ, P in S and Q on E, is special or non- 
special according as Q is a special or non-special point of E. 

THEOREM III. The non-special points of E constitute a point set of linear 
measure equal to the circumference of E, that is, of measure 2r. The proof 
of this theorem requires three lemmas. 

Lemma 1. Given any interval 6 of any stage, let the set of subintervals of 
5 obtained in the fourth following stage bec. Among the symbols of the set o 
occur all single symbols. Moreover, any given single symbol occurs such 
that the symbol of the interval of the preceding stage, in which it is contained, 
is not either of the adjacent symbols going around a polygon of the network n. 

Lemma 2. There exists a constant k > 0, such that the ratio of the length 
of any interval of any stage to the length of the interval in which 1t 1s contained 
in the preceding stage 1s greater than k. . 

Lemma 3. Given a finite symbol, a, there exists a constant c > 0, such 
that if 5 is any interval of E, there exisis a finite set of subintervals of 5, the 
points of which subintervals all have the symbol a in their symbols and such 
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that the sum of the lengths of these subintervals is greater than cé, where c 
depends on a, but not on 6. 

With the aid of this last lemma, Theorem III is readily proved. The 
proof of this theorem amounts to proving that the special points of E 
constitute a set of measure zero. A special point has the property that 
its symbol omits at least one finite symbol. The finite symbols constitute 
a denumerable set 


G1, M2, G3, .. 


Let that set of points of E, each of which has the property that its symbol 
omits the finite symbol a;, be denoted by A;. The special points of E 
are all included in the totality of point sets A;,i = 1, 2, .... Hence the 
statement of the theorem holds if each of the sets A;, 7 = 1, 2, ..., is 
proved to be a set of measure zero. The set A; is of measure zero if the 
complementary set B; is of measure 27. But the set B; is made up of a 
set of intervals on £, since those points whose symbols contain the symbol 
a; in a specified position constitute a set of intervals. 

Now consider the symbol a; and let c; be the corresponding c of Lemma 3. 
A proof is desired that the sum of the lengths of those intervals whose 
points have the property that their symbols contain a; is equal to 27. 
Those points whose symbols have a; at their beginning occupy an interval 
5, = S, of E. From Lemma 3 there are intervals having the desired 
property and lying in the interval 2x — S,, the sum of whose lengths is 
c(2r — S\). Let Se = S; + c;(2r — S,), and there are intervals of the 
desired property such that the sum of their lengths is S,. By a repetition 
of this process we obtain intervals having the desired property such that 
the sum of their lengths is 


Se = S,-1 + c;(2r a S,-1)- 


But S, is monotonic increasing and cannot exceed 27, hence it approaches 
a limit U, and we have 


U = U+¢,(2r — U) 


or U = 2x. The desired theorem is proved. 

A point Q of E determines with each point P(x, y) of S, a NE ray and 
hence a direction ¢ at P. The set of values of @ thus determined by Q 
constitute a function ¢9(x, y) which is analytic for (x,y) in S. If Qisa 
special point of EZ, from Theorem II, the points (x, y, ¢o(x, y)) of M are 
all special, and hence the special points of M can be considered as lying on 
a set of analytic surfaces. From Theorem III, given « > 0, the special 
points of E can be enclosed in a set of intervals 


0,0’, a — 1, 2, ae (1) 
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such that the sum of the lengths of these is less than «. But then the 
special points of M will be enclosed in a denumerable set of volumes 
bounded by the pairs of surfaces defined by 


goi(x, y), go'i(x, y), t — 1, 2, core (2) 


From the fact that the sum of the intervals (1) is less than e, it can be 
proved that the sum of the volumes (2) is less than ae, where a does not 
depend on «. Hence the special points of M constitute a point set of 
three-dimensional measure zero. We thus have the theorem: 

THEOREM IV. The non-special points of M constitute a point set of the 
measure of M. 

A complete geodesic, g, with end-points A and B on E, will be said to 
be non-special if, P being any point of g, the two rays PA and PB are 
non-special rays. Evidently, if this is true for one point of g, it is true 
for any point of g. 

Let the set of points M be that subset of the non-special points of M 
such that if (x, y, ¢) is a point of M, hence non-special, then (x, y, ¢’), 
¢’ = ¢ (mod n), is also a non-special point. The following theorem is then 
readily proved: 

Tueorem V. The set of points M is of the measure of M. 


1G. D. Birkhoff, ‘Dynamical Systems,’’ Amer. Math. Soc. Colloquium Publications, 
9, p. 248. 

2 J. Nielsen, ‘‘Untersuchungen zur Topologie der geschlossenen zweiseitigen Flachen,”’ 
Acta Mathematica, 50, 191-230 (1927). 

















VoL. 19, 1933 GEOPHYSICS: J. C. MERRIAM 349 


SYMPOSIUM ON CLIMATIC CYCLES 


NATIONAL ACADEMY OF SCIENCES, WASHINGTON, D. C., Aprit 26, 1932 


I. INTRODUCTORY REMARKS 


By Joun C. MERRIAM 


CHAIRMAN, COMMITTEE ON LONG-RANGE WEATHER FORECASTING, NATIONAL 
ACADEMY OF SCIENCES 


CARNEGIE INSTITUTION OF WASHINGTON 


The conference on climatic cycles is an outgrowth of discussion on ques- 
tions touching possibility of a periodic or cyclic element in earth climate. 
It is assumed that any information making it possible to foresee conditions 
of climate or of our physical surroundings is of importance. 

As the basis of regular routine in living we have already worked out 
with great accuracy cyles of day and night with varying length of the day. 
From history and from study of the earth’s movements we are acquainted 
with variation of the seasons. Is it possible that there are dependable 
periodicities which affect human living but are still unknown? 

The presence of periodicities or cyclic movements would not by them- 
selves guarantee the nature of what is commonly called weather. Fluctua- 
tions of temperature, of humidity, movements of the air and many other 
earth conditions we have not yet learned to predict with complete accuracy 
on the basis of the periodic influences which cause seasons. Other factors 
related to varying local conditions may affect the swing of general influences 
in such way as to erase their effect, at least temporarily. 

It is recognized that influences upon weather and upon the climate 
probably derive from a number of sources: We are, however, no longer 
greatly concerned about effects arising from the internal mass of the earth. 
In general the influences are seen as mainly external, and principally origi- 
nating from the sun. But with all that has been learned there is yet much 
to know regarding variation in solar and other radiation effects. We have 
also much to learn concerning the way in which the earth receives and, 
so to speak, utilizes the influences coming from the outside. 

In studies of periodicity we are concerned with the factors which have 
to do with reception of radiation, as also with the nature of fluctuations in 
the source of that radiation. This symposium represents largely an at- 
tempt to present from several points of view the facts regarding radiation 
and its influence upon the earth. 

The paper by Dr. Douglass, on climatic cycles as illustrated in tree rings, 
illustrates an interesting influence of radiation upon one of the most deli- 
cate instruments on the earth, namely, a growing organism recording in its 
tissues the development of its structure from period to period. The in- 
























GEOPHYSICS: A. E. DOUGLASS 





350 Proc. N. A. S. 





fluence of climate upon tree growth is an extremely important problem, and 
represents one of the most interesting means of registration of such factors. 
That there is change in plant life due to climate one cannot avoid seeing 
in development of vegetation at this time of the year. 

The paper by Dr. Walter S. Adams, on studies of the sun-spot cycle, 
presents the evidence relating to variation in the surface of the sun through 
observation of the sun itself, and a study of the elements involved in these 
changes. 

The paper by Dr. Abbot is based upon study of radiation as it is measured 
physically upon the surface of the earth. It stands between research on 
the sun and study of the sun’s influence upon the plant. 

The paper by Dr. Bowman presents evidence from history of the earth 
preserved in the geological record. One of the most interesting phenom- 
ena of time is the fact that such slight variation of temperature is noted 
through a period in earth history to be measured in many hundreds of 
millions of years, and yet in that time there have been fluctuations which 
marked periods of warm climate in northern latitudes, and glacial climates 
in what have in other times been relatively warm regions. 

The symposium is an attempt to present the results of these types of 
research as basis for discussion regarding possible further extension of our 
knowledge of the earth and of the influences which impinge upon life. 


II. EVIDENCES OF CYCLES IN TREE RING RECORDS 
By A. E. DouGLass 
UNIVERSITY OF ARIZONA, TUCSON, ARIZONA 


This paper describes some preliminary results in a search for a theory 
of climatic change, based on cycles. The search has been warranted by 
two possible advantages in approaching the subject: first, tree ring rec- 
ords, climatic in character, reaching a maximum length of three thousand 
years and accurately dated, whose climatic interpretation is based on 
ecological considerations and a limited number of actual tests; second, 
a cyclogram method of analysis of plotted curves, whose flexibility and 
speed of operation has made these studies possible. It is not intended 
in this paper to deal with the interpretation of these cycles beyond pointing 
out their apparent relationship to solar variations. 

The cyclogram is a pattern in which simple cycles are automatically 
separated from one another by a process which causes them to appear as 
rows of dots pointing in different directions. Each cyclogram includes the 
complete series of data and the cycle may appear anywhere within it: 
Thus it is independent of the beginning and ending of the sequence. Its 
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character is denoted by the regularity and conspicuousness of the align- 
ment of dots. We are therefore free to detect a cycle whether it is simple 
or composite, constant or variable in length, long or short in length, long 
or short in duration, subject to reversals of phast or subject to fractioniz- 
ing. The method also gives a first approximation to amplitude, often 
sufficient. By the use of some simple light-sensitive instrument the ampli- 
tudes can be measured quantitatively, or after the cycle is well delineated, 
its various elements may be determined by simple mathematical means. 
Thus in comparison with other methods the length of a cycle is determined 
in a few minutes with a precision approximating a least square solution 
(this was found in a test on variable stars reported to the American Astro- 
nomical Society in 1919); the beginning and ending of a cycle and com- 
position and possible variations appear at once, as never possible in pre- 
vious methods: the amplitude, however, is less precise, having been esti- 
mated on a scale of four, yet nearly every maximum must be present and 
there must be enough of them to give the impression of a straight line in 
the pattern. 

The capacity of this method may be seen in figure 1, which shows a 
multiple standard for calibration and its cyclogram. We can easily detect 
four of the five cycles of the original; in the analyzing instrument the 
fifth is seen as easily as these four. 

Discontinuities in a cycle become at once apparent as illustrated in 
figure 2 which shows the first cyclogram of sun-spot numbers, made in 1914. 
The bending of the horizontal fringes points out the discontinuities listed 
by Turner in 1913. It is evident that the solar cycle is not a precise and 
invariable period, yet we call it acycle. Therefore a cycle may be defined 
as the recurrence of similar conditions at similar intervals. The beginning 
and end of duration of a cycle may be observed with precision, which is 
sometimes of first importance, as will be illustrated in connection with the 
Hellmann cycle. (Figure 2 originally appeared in Astrophysical Journal, 
Vol. XLI, page 174, April, 1915.) 

Reliability Tests.—Natural sequences shown in tree growth may be dis- 
tinguished from accidental or ‘‘scrambled”’ sequences using the same values, 
if the sequence is long enough. In attempting to select the genuine out of 
three curves, two of which were accidental—this is called the ‘‘triangle’’ 
test—curves of 60 to 75 years in length could not be distinguished. Curves 
of 125 or more years met with better success. Some eighteen tests have 
been made on group curves of trees averaging about 175 years in length. 
Successful selection was made in every case. Over eighty per cent were 
easy. Two methods of selection were used; namely, smoothed curves 
and cyclogram analysis. Smoothing an accidental curve tends to reduce 
it to a straight line, while the natural sequence retains its larger variation. 
This is shown in figure 3. 
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PLATE 1 


1. Calibration curve (negative) and its cyclogram (cyclograms best 
seen by slanting view). 

2. Cyclogram of sun-spot numbers 1755-1914. 

3. Triangle test by smoothing: natural sequences on left, accidental 
sequences on right. 

4. Triangle test by cyclograms: natural sequence top row, others 
accidental. 

5. Independent cycle analyses of western pine records (a and b), mean 
(c and e), sun-spot cycle ratios (d). 
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The cyclogram test as shown in figure 4 consists of three cyclograms in 
four different settings. One of the three horizontal rows is genuine, the 
others are accidental. The genuine is readily distinguished by its clearer 
background denoting fewer short cycles. 

Thus we see that natural sequences tend to exhibit longer cycles than 
accidental sequences. It will be noticed that there is no attempt here to 
distinguish the genuine by its exact cycle length. Eventually that may 
become the easiest method. These tests show that there is something in 
the curves besides accidents. 

Western Zone Cycles.—Admitting the complexity of cycles in tree ring 
records and the probability that some are accidental, the question is, can 
we by quantity methods find some order in the apparent confusion. Pre- 
liminary results have been obtained by analysis of a large number of western 
tree records. ; 

More than 52,000 rings in 305 trees, mostly pine from regions between 
the Rocky Mountains and the coast, have been measured and analyzed. 
This was done in 1926 and before the close of that year it was noted that a 
relation seemed to exist between the cycles observed and the sunspot cycle. 
Fearing that the personal equation had not been completely eliminated, 84 
curves were replotted in a scale unknown to me and the analyzing process 
carried through a second time before publication. The second result 
agreeing with the previous one is shown in figures 5a and 6. Mean cor- 
relation coefficient + 0.53 + 0.06. 

The cycles were summated as shown in 5c and e compared with simple 
fractions of 34 years which is closely three times the sun-spot cycled. 
Speaking generally, the cycles found in tree ring records seem to be very 
nearly simple fractions of small multiples of the sun-spot cycle. Similar 
expressions have been reached by Abbot, Clayton and C. E. P. Brooks. 

Hellmann Cycle—The most important fractional part of. the sun-spot 
cycle found in trees is the '/2 value, 5.5 years to 5.7 years. It is believed 
that this was first presented in its relation to the sun-spot cycle by Hell- 
mann, of the Prussian Meteorological Bureau, in 1906, in a study of 
North German rainfall. It is therefore called by his name. His curve 
is shown in figure 6. 

It is readily found at other places. A California rainfall curve expressly 
presented by an author to show the absence of cycles is shown in figure 7. 
On being smoothed it seems to show this Hellmann cycle. Horton in 1903 
mentions this short cycle in California rain. Its length and association 
with the sun-spot cycle was published by the present writer in 1909. It 
has been plotted recently by Nicholson, who, however, does not attach 
any special significance to it. It is in use by Gorton at the Scripps Institu- 
tion of Oceanography. It appears in many centuries of sequoia growth. 

About 1914 the Hellmann cycle was recognized as a marked character- 
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PLATE 2 
6. The Hellmann relation. 
7. Similar relation in San Diego rainfall. 
8. Hellmann cycle in early Arizona trees. 
9. Same in Grand Canyon trees, 1700-1920, isolated maxima. 
10. Hellmann effect in 80 European trees. 


istic of Arizona trees in the fifteenth century. Two different curves of 
tree growth are shown in figure 8; the maxima are marked to bring out 
the Hellmann cycle and a suggested position of the sun-spot cycle is placed 
below. This cycle was readily seen in the cyclograph in early analytical 
work. It was further noticed that from 1650 to near 1725 the eleven- 
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year cycle was absent. About 1725 it returned and after 1800 began again 
to show the double crested character. This failure near 1700 caused much 
doubt of the possible relation to sun-spot records but an account of it 
was published in 1919. Three years later Maunder’s description of the 
great dearth of sun-spots from 1645 to 1715 was published and offered an 
explanation of the absence of the cycle in Arizona trees. It was the ability 
of the cyclogram to show beginning and ending of a cycle that drew at- 
tention to this coincidence. The story may be seen in ordinary curves as 
published in 1919 (Climatic Cycles and Tree Growth, page 103; Carnegie 
Institution. The relationship of Arizona and California tree growth to 
southern California rainfall and temperature and to the sun-spot cycle was 
given in the same volume on page 104). The Hellmann cycle shows in 
the means of the Arizona groups including 58 trees, but it is necessary first 
to remove a 19-year cycle. The Grand Canyon group, however, shows 
the cycle very clearly in “isolated maxima” cyclogram of figure 9. 

In north Europe certain groups of trees show an excellent sun-spot cycle 
with a tendency to introduce an added maximum at time of sun-spot mini- 
mum. In a qualitative result in 80 European trees in which the number 
of occurrences of maximum growth in each year were plotted, a strong 
maximum corresponded to the sun-spot maximum and a lesser one appears 
at sun-spot minimum, as shown in figure 10. 

Sun-Spot Cycle——European trees near the Baltic show increased growth 
during sun-spot maxima. This is illustrated in the two sections shown in 
figures 11 and 12 The sun-spot maxima, 1830, 1837, 1848, 1860, 1871, 
1884, 1895 and 1905 are marked with dots. On plotting the ring growth 
of 57 trees of north Europe, mostly pines, a curve is obtained which 
closely resembles the curve of sun-spot numbers as shown in figure 13.* 
Integration of this tree growth on a 11.4-year cycle is almost identical with 
the similar integration of the sun-spot cycle. It has been suggested that 
these effects are due to thinning the forest with German regularity.. Ona 
recent visit there I was assured that the forest is thinned every three years 
or so for instruction purposes in the University. 

The Arizona trees have a tendency to give a cycle of about 14 years. 
The cyclogram of these trees in figure 14 shows the Hellmann relation in 
the left quarter, which is the early 125 years of its length. In the central 
part, near 1700, it shows an approximate 10-year cycle during the dearth 
of sun-spots, and in the right-hand third something close to 21 years. 

Early analysis of the Flagstaff 500-year tree records (1909) gave a cycle 
of about 21 years since 1700. This has remained very prominent in the 
last 200 years of Arizona tree growth. In various curves of cycle inten- 

* Figure 10 has appeared in Climatic Cycles and Tree Growth (Carnegie Institution: 


I, 1919; II, 1928) I, page 78: Figure 11, Jbid., I, Plate 8A: figure 12, II, Plate 8: figure 
13, I, page 77: figure 14, II, Plate 9, no. 7. 
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11. Eberswalde (German) pine showing large growth at sun-spot maxima. 
12. Swedish spruce showing same. 
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sity (periodograms) a considerable crest from 19 to 22 has been evident. 
There has always been uncertainty regarding the real length of this cycle. 
A series of integrations at 18, 19, 20, 21 and 22 years in figure 15 give 20 
years as a good average value of this cycle. The best mean values from 
1700 to 1920 were recently sent to Dinsmore Alter, who subjected them to 
periodogram analysis in a rapid and skilful mathematical attack. His 
correlation periodogram of this sequence is given in figure 16, which shows 
a strong crest at something over 20 years, with further conspicuous crests 
of about 40 and 80 and 120. A cyclogram made of the same data is repre- 
sented in figure 17, which shows at once that this 20-year cycle is composite 
during this 220-year interval. Through the first one-third, a cycle of about 
23 years is prominent and in the remainder a strong 19-year cycle appears. 
One recalls that 19 is approximately one-third of five times the sun-spot 
cycle. Thus this puzzling but strong cycle proves to be a composite of 
two taken in succession. 

A curve of Arizona tree growth plotted at ‘/; horizontal scale for 1200 
years shows approximately five times the sun-spot cycle from 1168 to 1512 
A. D. 

I have illustrated above a number of cycles but it will be seen that 
each one is a simple fraction of a small multiple of the 11-year sun-spot 
cycle. 

Sequoia Integrations.—The difficulty with the study of cycles in climate 
has been the limited lifetime or duration of any one cycle. One thinks he 
has a cycle well tied down only to find that it has vanished and another 
has taken its place. The possibility of finding a definite succession of 
cycles has long been recognized and the great sequoias have seemed the 
most promising location for the search. We have in our laboratory 50 
radials of these trees, of which about half show 2000 rings or more, and four 
extend back 3100 years; one reaches 3237 years, to 1305 B.C. They fall 
into two groups, the Grant Park sequoias and the Springville sequoias, 
some forty miles apart. It was discovered about 1921 that the double 
sun-spot cycle 22 or 23 years is more common in these trees than the 11-year 
cycle. The pine tree curve of the Sierra Nevadas resembles the sequoias 
in giving prominence to a cycle of about 23 years. A correlation periodo- 
gram of this series of data by Dr. Alter, reproduced in figure 18, shows a 
pronounced crest at about 221/2 years and its multiples, thus confirming 
these results. 

In 1927 a study was made of certain portions of the sequoia records in 
whith the presence of the 11-year cycle had been suspected. Integrations 
on a basis of 23 years gave immediate results in the form of a recurrence 
cycle of 260 to 280 years, in which 23/5 or 4.5 years repeated itself. At 
the present time the same form of integration has been extended so that 
now we have a series of curves each 23 years long, that are consecutive 
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PLATE 4 


13. Fifty-seven European trees and sun-spot cycle. Below: 13 
Eberswalde pines and 1 Swedish spruce. 

14. Cyclogram of Arizona tree record showing succession of cycles. 

15. Arizona 20-year cycle by integration at 18, 19, 20, 21 and 22 
years. 

16. Alter’s correlation periodogram of same, 0 to 130 years. 

17. Cyclogram of same showing succession of 23- and 19-year cycles. 

18. Alter’s correlation periodogram of Sierra Nevada pines showing 
22.5-year cycle at dots and its subdivision. 


and continuous through nearly 3000 years and each curve is a running or 
overlapping mean of three successive 23-year intervals. 

A large percentage of these 130 curves show clearly a 23-year cycle or 
some subdivision of it up to five equal parts. The division into 4 parts is 
illustrated in figure 19. Of course this work is far from finished but these 
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Hellmann cycle, (a quarter of 23 years) in California sequoias, 
by integration method, A. D. 155 to 146. 

Similar values drawn by lot and treated in same manner. 

Hellmann cycle in sequoias by same method, 1719 to 1914. 

Possible 275-year recurrence cycle in 2100 years of sequoia 
records (indicated in horizontal alignment). 

Average cycle succession in Grant Park sequoias in 2100 years. 
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results apparently agree in the two separated sequoia groves. Limited 
checking has also been carried successfully to other sequoia groupings that 
for various reasons have been segregated. As a matter of safety a couple 
of hundred years of sequoia values were ‘‘scrambled”’ and then treated by 
the same integration method, as shown in figure 20. This brought out 
the fact that the manner of smoothing eliminates two-year cycles which in 
turn leaves 3-year cycles, or 8 crests in the curve, most numerous, 7 crests 
in the next and soon. One scrambled curve has five crests, and one has 
four, showing that some of the fractionizing in the real curve is accidental, 
especially in isolated cases. In the scrambled curves the same fraction 
does not occur twice in succession. But in the natural curves the same 
occurs three to five or more times in succession and the chances are greatly 
against accident. As the length of the fraction increases, for example from 
41/, to 71/2 years, accident becomes less and less likely. So even three 
curves in succession showing triple subdivision, 71/2 years, are very proba- 
bly real. 

The Hellmann relation is fairly evident in upland sequoias for the same 
time, as appears in figure 21, which was developed like those above. It 
is interesting to note that in the large groups of sequoias it seems to extend 
back to about 1750. 

In these 130 curves, covering the age of the sequoias, a recurrence cycle 
of about 275 years seems to be a common character as shown in this pre- 
liminary study, three examples of which are given in the horizontal bands 
in figure 22. A cycle about 100 years long is common to them all and some- 
thing near 220 years is also there. All these are mentioned by Turner. 
Something near 300 years is discussed by Clough. Dr. Alter mentions 250 
years. The double value near 560 years also appears and it is possible 
that it will prove important. In order to visualize the possibilities, the 
275-year length is taken as a workable solution for the time being. A pro- 
visional pattern based on this 275-year succession is shown in figure 23. 
This diagram shows the sequence of different short cycles in the 275-year 
pattern derived from 2000 years of the Grant Park sequoias and is perhaps 
the most promising sign at the present time. In it we find that the per- 
sistence of the Hellmann relation since 1750 agrees with previous recur- 
rences of that cycle. If the various cycles continue as they have done we 
have some reason to expect that through the next 70 years this Hellmann 
relation should be less conspicuous or absent, with a probability of re- 
placement by a 23-year cycle or the same divided into five equal 
parts. 

This pattern therefore encourages us to hope that from the giant se- 
quoias, checked by historical material, we shall obtain in some form of 
cycle succession a real step forward toward a theory of climatic change and 
long range prediction. 
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III, PERIODICITY IN SOLAR VARIATION 
By C. G. ABBOT AND Mrs. A. M. Bonp 


SMITHSONIAN INSTITUTION, WASHINGTON, D. C. 


Long ago, Secretary Langley induced Congress to support the study of 
solar radiation at the Smithsonian Astrophysical Observatory. He pointed 
out that all life and all weather depend on it. He held out the possibility 
and hope that a sufficient knowledge of solar radiation, and of its behavior 
in our atmosphere, might even enable meteorologists to forecast long in 
advance the fat years and the lean years as Joseph is said to have done in 
Egypt. 

After 40 years of research, we have results which seem to us to justify in 
some degree Langley’s hope. We have not yet, it is true, tried the bold 
venture of long-range forecasting, but we have evidence to present to the 
Academy today that the sun’s output of radiation is variable; that its 
variation is periodic; that the United States weather departures from 
normal are periodic; and that a considerable part of the ranges of weather 
departures from normal are comprised in a series of periodicities which 
are identical with those found in the sun. 

For more than 25 years the staff of the Smithsonian Astrophysical 
Observatory has been measuring the intensity of solar radiation. At first, 
in Washington, we further developed the method devised by Langley and 
used by him about 50 years ago at Allegheny and at Mount Whitney. 
We devised the silver-dise pyrheliometer for ordinary daily measurements 
of the total intensity of solar radiation at the station. We also devised 
the water-flow and the water-stir standard pyrheliometers, whereby we 
reduced the scale of measurement to standard calories per square centimeter 
per minute. We improved the recording spectrobolometer of Langley 
so that in less than 10 minutes it could furnish an excellent record of the 
intensities in the solar spectrum between wave-lengths about 0.35 microns 
in the ultra-violet and 2.5 microns in the infra-red. ‘We devised graphical 
methods operated by instrumental appliances whereby we were able to 
compute the intensity and distribution in the spectrum of the solar radia- 
tion as it is outside the atmosphere, whenever we determined the intensity 
and distribution of it as it reached the surface observatory at different 
solar altitudes. We devised an instrument for measuring the brightness 
of the sky about the sun which we named the pyranometer. By its aid 
we have devised a brief empirical method for estimating the atmospheric 
transparency in all observed wave-lengths. We have also devised a 
spectroscopic method for estimating the quantity of precipitable water 
held in the form of vapor in the atmosphere. From our determinations 
of atmospheric transparency we have checked exactly with other methods 
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on the determination of the number of molecules per unit volume of gas. 

’ As the temperature of the earth and the fundamental factors of climate 
and weather depend on the intensity of solar radiation, we have made 
earnest efforts over a long period of years to secure accurate measurements 
of it. When we began this work in 1903, authorities were in doubt over 
the entire range as between Pouillet’s value of 1.76 calories, and Angstrém’s 
value of 4.0 calories for the solar constant of radiation. As a result of 
our work, carried on at all seasons—at stations ranging from sea-level to 
4500 meters elevation, checked by automatic apparatus exposed from 
sounding balloons at 25,000 meters elevation—there is now no doubt 











FIGURE 1 
Solar radiation station, Montezuma, Chile. 


that the true value lies certainly within one per cent of 1.94 calories per 
square centimeter per minute. 

We have discovered evidences of variability of the sun’s emission. 
Having devised a brief method for measuring the solar constant, we have 
applied it several times a day on all favorable days over quite a term of 
years. We have occupied mountain stations in desert lands in Arizona 
and Southern California, in Northern Chile and in South West Africa. 
Figure 1 shows a view of the apparatus at our station at Mount Montezuma 
in Northern Chile, 9000 feet above sea-level. The pyrheliometers and the 
pyranometer are exposed outside, and the solar altitude is measured with 
a theodolite. A beam of light is reflected into a cave observatory where 
the spectrobolometric work is done. Figure 2 shows the close accord at- 
tained in the monthly mean values of the solar constant at three widely 
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FIGURE 2 


Agreement of monthly mean solar constant values. Three stations. 
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FIGURE 3 


Agreement of solar constant values. 
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separated stations. It is clear that if the observations at the earth’s 
surface and the estimates of losses in the earth’s atmosphere could be 
correctly measured, then determinations of the solar constant (that is, the 
intensity of solar radiation outside the atmosphere) ought to agree exactly 
wherever made on the earth’s surface. In actual practice we have so far 
refined our determinations that our two best widely separated observa- 
tories, Montezuma, Chile and Table Mountain, California, do agree in 
their monthly mean values over a period of five years within an average 
difference of 0.08 per cent. The probable error of the mean curve shown in 
figure 2 is well below 0.1 per cent. 

It will be noted that the three stations not only agree closely, but unite 
to indicate fluctuations of the sun’s emission. The extreme range of varia- 
tion shown in figure 2 is 1.2 per cent. On an earlier occasion, in 1922, 
a range of the monthly mean values of nearly 3 per cent was observed, as 
indicated in figure 3, where values from Montezuma, Chile and Harqua 
Hala, Arizona, are compared. Taking the best results of the work as 
derived from the evidence of all stations, we find the variation of monthly 
mean solar constant values since 1920 as indicated in figure 4, curve A. 
We are able to reproduce it as the sum of seven regular periodicities of 
7, 8, 11, 21, 25, 45 and 68 months. The degree of approximation is shown 
by the smallness of the residuals in curve B. It is to be noted that larger 
residuals are found in the earlier years when the solar observing was less 
perfected than it became later. 

These results on periodicities have been obtained by Mrs. Bond with an 
instrument which we call the periodometer, which is on exhibition in the 
adjoining room.* It was constructed with the aid of a grant of $1000 from 
Research Corporation of New York. Its purpose is to discover and 
evaluate periodicities in long series of observations. It does not recognize 
the reality of any period until tested, and it evaluates its distribution in 
amplitude without regard to any assumed mathematical expression. It 
appears to us, for instance, that as the curve of sun-spot frequency is well 
known not to be of regular sine form, there is no reason to suppose that 
other solar periodicities should have a sine form. Hence, our instrument 
is designed to evaluate their forms as Nature fixed them, not according 
to the forms assumed in mathematical series and harmonic analyzing 
machines. 

Curves C, D, E, F, G, H, of figure 4, show the periodicities actually dis- 
covered in the solar radiation by aid of the periodometer. It will be seen 
that the 21-month period betrays also one of 7 months. In the cases of 
the shorter periods, we have been able to separate the data into several 
groups and thus to evaluate the periodicities independently at several 
epochs. These partial determinations are shown in curves Ci, C2, C3, 
* As shown at the Academy, April, 1932. 
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Dj, Dz, D3. In such cases we have been encouraged to find that the maxima 
and minima occur without change of phase in these independent epochs. 
Thus we regard the periodicities found as having reality and permanence. 
We ventured in November, 1930, to make a forecast for 1931 and 1932 
of the probable march of solar variation. It has been fairly well verified, * 


* . ” 





FIGURE 4 


Analysis of solar variation. A, original data. C, D, E, F, G, H, perio- 
dicities of 8, 11, 25, 45, 68 and 21 months. B, residual, after removing peri- 
odicities. C,, Co, C3, D,, De, Ds, partial analyses leading to curves C, D. 


although it called for solar constant values almost all the time since 1930 
about one per cent above the mean, notwithstanding that the values pre- 
ceding the date of forecast for several years had been prevailingly below 
the mean. 

It has been of great interest to us to note that several of the periodicities 


* Up to December, 1932, as shown by the dotted curve, figure 5. 
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found in solar variation are closely related to the sun-spot period of 111/, 
years or 135 months. Thus, 68 months is within its probable error one- 
half, 45 months one-third of 135 months. Again if we take a period ap- 
proximately three times as long, or 400 months, which is near the Bruckner 
period, 25 months is one-sixteenth, 21 months closely one-nineteenth, 
11 months closely one-thirty-sixth, 8 months is one-fiftieth and 7 months 
closely one-fifty-seventh of its duration. 

If we admit provisionally (subject to the findings of subsequent years) 
that the solar variation is made up of the seven periodicities named, it 
becomes of interest to see if these same periodicities are traceable in tem- 
perature departures of the weather. We have investigated this question 
for three widely separated United States stations, viz.: Clanton, Ala- 
bama, Washington, D. C., and Williston, North Dakota. We have taken 
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Prediction of march of solar variation made November, 1930, and fulfilment. 


our data from the climatological summaries of the United States Weather 
Bureau, 1918 to 1930. In order to eliminate the yearly march of tempera- 
ture, we have computed for each station the mean monthly temperatures 
1918-1930, and have subtracted them from the observed, thus giving 
monthly departures exactly suited to the epoch studied. Lest the in- 
fluences of shorter period changes should obscure the general march of 
events, we have smoothed the monthly temperature departures by taking 
five-month consecutive means of the form: 


1 + Op + Og + Oy t+ As Oy + Os + A + Gs + Oe 


et cetera. 





9) o 


With the data thus prepared, we have sought and evaluated with the 
periodometer all the periodicities which the curves disclosed. Our pro- 
cedure, as in the case of solar variation, is to subtract from the data the 
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effect of each periodicity as soon as determined, before proceeding to 
evaluate in the residual curve another periodicity. We continued the 
search and evaluation until no more periodicities could be perceived. 
The result obtained for Clanton is shown in figure 6. Periodicities of 
8, 91/2, 11, 16, 18, 21, 25, 33, 45 and 66 months were evaluated. These 
periodicities and their partial determining curves are indicated by letters 
C, D, E, F, G, H, I, J, K, L of figure 6. The residual shown in curve B 
plainly indicates the 11'/,-year period. We also note the large positive 





FIGURE 6 


Analysis of departures from normal temperature, Clanton, Alabama, and residual, B, 
after removing periodicities found from smoothed curve of temperature departure A. 


departure shown in the residual curve B for the year 1930, a year remark- 
able for the extraordinary drought and accompanying cloudlessness. A 
similar extraordinary positive departure for 1930 was found for Washing- 
ton, and also for Williston. It will be noted that strong periodicities 
of 8, 21, 25, 45 and 66 or 68 months found in Clanton temperatures are 
found also in solar variation. The total interval is too short to distinguish 
the last period exactly. The 11-month solar period is feeble in Clanton 
temperatures. The 135-month period is doubtless of solar origin, although 
it does not appear conspicuously in the solar variation between 1918 and 
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1930. The other Clanton temperature periodicities of 9'/2, 16, 18 and 33 
months were not found in the sun, but nevertheless 33 months is one-half 
of 66 months, which is conspicuously found as a periodicity in the sun, 
and 16 months is twice 8 months, a solar period. 

At Washington periodicities of 8, 9/2, 13'/2, 18, 25, 45 and 68 months 
are found. The residual curve shows clearly the 135-month periodicity 
in practically the same phase, though lesser amplitude, then Clanton. 
The extraordinary drought of 1930 produces its strong positive departure. 
Here again the strong periodicities of 8, 25, 45 and 68 months seem to 
reflect solar radiation changes. The periodicities of 91/2 and 18 months 
are found also at Clanton. The 13'/.-month periodicity is new. No 
appreciable influence of the 11-month solar periodicity is found. 

Williston, North Dakota, has a much wider range of departures than 
Clanton or Washington. Yet, this range is reduced nearly to zero by with- 
drawing the periodic departures discovered. They are, respectively, of 
7, 8, 11, 13/2, 18, 21, 25, 28, 45 and 68 months periods. The 111/,-year 
period is but indistinctly shown in the residual curve and has a very 
different phase from corresponding curves at Clanton and Washington. 

In order to fix our ideas of the relations between solar and terrestrial 
periodicities which we have discovered, we give in table 1 a summary of 


TABLE 1 
AMPLITUDES OF PERIODICITIES 
PERIODS OF 
SOLAR ORIGIN SUN CLANTON WASHINGTON WILLISTON 
Months Calories ° Fahr. ° Fahr. ° Fahr. 
7 0.005 ae ae 1.2 
8 0.005 0.7 0.7 1.0 
11 0.009 aM 1.3 
21 0.004 2.0 3 0.7 
25 0.010 1.6 WS 2.7 
45 0.013 1.8 2.2 4.0 
68 0.014 1.8 2.0 5.5 
135 2.0 1.5 0.7 
Sum of Amplitudes, Solar 9.9 8.1 i a | 
PERIODS OF 
TERRESTRIAL ORIGIN 
Months 
91/2 Ss 1.0 0.5 ps 

131/2 bie pay 1.3 3.0 
18 anit 2.1 LZ 2.4 
28 aa shag 3.1 

34 ie 2.2 

39 Ss 1.3 


Sum of Amplitudes, Terrestrial 
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FIGURE 7 
Curves of variation from normals. A of solar radiation, B,C,D, from normal 
temperatures of Clanton, Washington and Williston. Curves E,F,G,H, and 
I,J,K,L, are 25-and 45-month periodicities in these four curves A,B,C,D. 
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them.* We invite attention to the fact that a majority of the periodicities 
in terrestrial temperatures which we have found are identical in length 
with periodicities in solar variation. The sum of the maximum ampli- 
tudes of these periodicities which may be of solar causation is much greater 
than the sum of those which appear to be of terrestrial origin. In fact, 
if we make the hypothesis that a terrestrial temperature departure perio- 
dicity is of solar causation if it is of identical period with a solar periodicity, 
then we shall come to the result that a major part of the ranges of terres- 
trial temperature departures from normal used in our analyses are due 
to variations of solar radiation. 

We do not venture to claim this important conclusion as yet. We 
_ intend to carry on the research much longer. But at least the investiga- 
tion has decided promise. 

We have been especially interested to compare among themselves the 
curves giving the solar variation and the temperature departures at the 
three stations. We have also wished to compare the curves giving periodici- 
ties of 25 and 45 months in solar radiation with the corresponding tempera- 
ture periodicities at the three stations. These comparisons are shown in 
figure 7. 

It appears at sight that parts of the curves of temperature departures 
for Williston and Washington are very similar, but that the similarity is 
slight as between Williston and Clanton. On the other hand, there are 
many points of similarity between the temperature departures of Washing- 
ton and Clanton. The large departures of similar form found at Williston 
and Washington in the years 1918 to 1921 occur from one to two months 
later at Washington than at Williston. 

The reader’s attention is invited to a correlation of form between curves 
A and D, and to the close synchronization not only in that but in the cor- 
relations of curves Hand HandJand JL. Also there is a close synchroniza- 
tion of correlation between curves F and G and between curves J and K. 
The reader need not be reminded that all of these analyses are completely 
independent. 


* The results for Clanton in table 1 and in figure 7 are taken from an earlier analysis 
than that shown in figure 6. 
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IV. THE NATURE OF THE SOLAR CYCLE 
By W. S. ADAMS AND S. B. NICHOLSON 


MountT WILSON OBSERVATORY, CARNEGIE INSTITUTION OF WASHINGTON 


Correlations of sun-spots with the weather were attempted as early as 
1801 by Sir William Herschel. At that time the periodic nature of sun- 
spots had not yet been recognized, although spots had been observed tele- 
scopically for nearly two hundred years, and the available weather records 
were so poor that Herschel considered the price of wheat a better index of 
the weather than the records of the meteorologists. He found a slight 
correlation but pointed out its uncertainty and warned anyone from using 
it for more than it was worth. 

Since the days of Herschel our knowledge of sun-spots and the weather 
has been much improved, and it is the aim of this paper to call attention 
to the astronomical facts relating to the nature of the solar cycle which 
should be recognized when correlations with sun-spots are being investi- 
gated. 

Although sun-spots are cooler and faculae hotter than the normal solar 
surface, their areas are never large enough to affect directly the integrated 
solar radiation by more than a fraction of one per cent. Significant varia- 
tions in the total solar radiation must come from changes in the sun’s 
average effective temperature which are smaller than those observed in 
spots and faculae but which affect much larger areas, perhaps the whole 
solar surface. Such changes may be brought about by variations in the 
transparency of the solar atmosphere, as indicated by the observations 
of the Smithsonian Astrophysical Observatory. Sun-spots, faculae and 
prominences are therefore to be considered merely as indices of some more 
fundamental solar change. 

The spectrum occasionally reveals in the neighborhood of spots small 
areas which radiate abnormally, not as black bodies at a higher temperature 
but principally as regions in which the intensity of hydrogen emission is 
exceptional. The effect of the sun on the earth’s magnetic field and on the 
aurora is probably due in some way to such areas, and it is not inconceiv- 
able that the transparency of the earth’s atmosphere may also be changed 
sufficiently to affect the amount and quality of the transmitted solar radia- 
tion, although such an influence has not been established. 

The numbers and areas of sun-spots, faculae and prominences vary in 
a fundamental cycle of about eleven years. Both the period and the ampli- 
tude of this cycle vary independently, the period between nine and four- 
teen years, and the amplitude by about fifty per cent of its average value. 
The true cycle is from minimum to minimum, not from maximum to maxi- 
mum. This fact is shown most strikingly by the magnetic polarities of 
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sun-spots. At the last two minima the distribution of magnetic polarities 
in the spot-groups has been observed to reverse, spots throughout one cycle 
having polarities opposite to those of the preceding and following cycles. 
The unity of the cycle from minimum to minimum is also indicated by the 
movement of the zones of spot activity. After minimum, when the number 
of spots is increasing, they occur in two zones on either side of the solar 
equator in average latitude about twenty-five or thirty degrees. The lati- 
tude of these zones decreases continuously throughout the cycle to about 
ten degrees. The renewed activity of the next cycle begins again in the 
higher latitudes, a real discontinuity occurring at the minimum. 
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FIGURE 1 


Wolf Sun-Spot Numbers. The two upper curves give the yearly means 
of. the smoothed Wolf Numbers. The lower curves give the observed 
monthly values (light line) and the smoothed monthly values (heavy 
line) for the last two cycles. 


The reversal of polarity at the minimum indicates a fundamental period 
twice that of the eleven-year cycle. Further evidence of a fundamental 
difference in alternate cycles is shown in the last ten or twelve cycles by 
a difference in their intensity, but the difference is such a small fraction of 
the eleven-year amplitude that there seems to be no reason for considering 
the double period in searching for first-order correlations. A long period 
of between sixty and eighty years is also indicated by changes in the 
amplitude; but if this is real the present cycle should have been more ac- 
tive. The next cycle will certainly furnish a critical test of the reality of 
this long period. 

Sun-spots are often distributed non-uniformly in solar longitude so that 
the number visible at any time depends on the rotation period of the sun, 
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which approximates a month. The amplitude of this apparent varia- 
tion may be as great as the average amplitude of the eleven-year cycle. 
Since the synodic solar rotation period is about one month, the monthly 
means form a good index of the average activity on the whole solar surface. 
When these means are examined we see that the eleven-year cycle does not 
progress smoothly but in pulsations, the amplitudes of which are about one- 
fourth that of the main cycle, and the periods from six to fifteen months. 
These short-period pulsations often persist throughout a whole cycle with 
approximately the same period. They are equally apparent in the curves 
of faculae and other solar phenomena, and it is because of them that the 
exact time of any phase of the main cycle is so difficult to determine. 

In 1847 Dr. Rudolf Wolf began a series of solar observations in which he 
expressed the solar activity by a number, called the Wolf Relative Sun- 
Spot Number, which is equal to k(m + 10 g), where m is the number of indi- 
vidual spots, g the number of groups and k a factor depending on the 
observing conditions. The value of k for Mount Wilson is about 0.7. 
Wolf collected all the earlier observations of sun-spots and from them esti- 
mated the Wolf Numbers as far back as 1749. The first half of this series 
is therefore not so reliable as the last half, which has been made from daily 
observations taken for the express purpose. . This series of Sun-Spot 
Numbers, covering nearly two centuries, is the most homogeneous index 
of solar activity over such a long interval and is the most convenient to 
use in studying correlations. 

Many other criteria may be used as indices of solar activity. The areas 
of spots and faculae, measured at Greenwich since 1873, form an excellent 
series for this purpose; but large individual groups sometimes affect the 
mean areas to such an extent that most investigators consider the Wolf 
Numbers a better index of solar activity. Character figures based on 
bright hydrogen and calcium flocculi and on prominences are available 
for dates since 1923 and will be extended at least as far back as 1913. This 
series is still too short to be of much value in studying correlations with the 
eleven-year period. A comparison of the curves for these different criteria 
shows that on the whole they are very similar, although some minor differ- 
ences are evident. 

Dr. Abbot’s measures indicate that the observed changes in the solar 
radiation are principally in the violet part of the spectrum, and it was with 
this in mind that measures of the ratio of ultra-violet to green sunlight 
have been made at the Mount Wilson Observatory since 1924. These 
measures show a good correlation with sun-spots except in the year 1929, 
but the series is still too short to give much weight to the comparison. 

When only a short series of observations is available for studying correla- 
tions, especially at the time of maximum solar activity, the principal solar 
fluctuations are those of the short-period pulsations, which may differ funda- 
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mentally from those of the eleven-year cycle. That the relationship 
between sun-spots and the measures of total solar radiation is complicated 
may be seen from a casual comparison of curves representing their values. 
Whatever the relationship may be, it is obviously not adequately de- 
scribed by saying that high values of the solar radiation are associated 
with the sun-spot maximum and low values with the minimum. 
Correlations of sun-spots have been made with almost every conceiv- 
able variable, especially those pertaining to the weather. Many writers 


FIGURE 2 


Spectroheliograms of calcium flocculi overlying sun-spots showing changes in area and 
latitude at different phases of sun-spot cycle. 


have stated that the best correlation with the weather is that between sun- 
spots and tropical temperatures. This correlation is therefore chosen 
as an example of what may be expected from such studies. Larger cor- 
relation coefficients have been found in other comparisons where the time 
interval is short, but generally a longer series of observations shows that 
the material has been insufficient for such a statistical study. 

The mean annual temperature of six tropical stations, Batavia, Chocin, 
Colombo, Freetown, Nuwara Eliya and Trincomalee, was obtained from 
World Weather Records, Volume 79 of the Smithsonian Miscellaneous 
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Publications. This series, fifty-eight years in length, was compared with 
the Wolf Sun-Spot Numbers for the corresponding years. When combined 
in cycles equal to the solar cycle, a mean amplitude of 0°4 Fahrenheit was 
found, with the maximum temperature at the minimum spot phase. The 
correlation coefficient between the annual Wolf Number and the annual 
temperature of the six stations is —0.37. The maximum deviation for 
any year from the mean of the fifty-eight is 0°9, and the probable error made 
in assuming that the temperature for any year is equal to the mean is 0°31. 
If a curve with the period of the solar cycle is assumed, the mean temperature 
for any year can be assigned with a probable error of 0°26. This reduction 
of only 0°05 in the probable error illustrates the uncertainty in the mean 
amplitude of 0°4 which was obtained from the solar cycle. These prob- 
able errors can be expressed in another way. We may predict that the 
mean yearly temperature for these six stations will be equal to its mean 
for the last fifty-eight years and give one to one odds that the prediction 
will not differ from the observed value by more than three-tenths of a de- 
gree. On the other hand, if we utilize the solar cycle to improve the pre- 
diction, the odds are about three to two that the prediction will not be in 
error by more than three-tenths of a degree. The gain through the intro- 
duction of the sun-spot cycle is quite appreciable though not striking. 

Since mean temperatures from other stations also show a negative corre- 
lation with sun-spots, the correlation is probably real, in spite of its low 
coefficient, but we should follow Sir William Herschel’s example and warn 
the weather prophets not to use it for more than it is worth. 

In conclusion, it may be said that of the various correlations with solar 
activity which have been attempted that between variations in terrestrial 
magnetism and sun-spot activity is by far the best established. The paral- 
lelism of the curves is conclusive in spite of the fact that frequently mag- 
netic storms on the earth occur when no unusual activity is seen on the 
visible disc of the sun. The correlation of ultra-violet radiation with sun- 
spot frequency seems probable, although observations have extended over 
too short a period to give fully definitive results. The relationship of 
measures of solar radiation to solar activity is still uncertain, and correla- 
tions of rainfall and temperature with sun-spot numbers are at least as 
doubtful as are short-period climatic cycles themselves. That important 
connections between solar and terrestrial phenomena still remain to be dis- 
covered can hardly be questioned; but when we consider that very ac- 
curate measurements of solar radiation are available only for the past 
twelve or fifteen years, of ultra-violet radiation for eight years and suit- 
able weather records for not more than fifty to seventy-five years, it is not 
surprising that many of the attempts to establish correlations are still quite 
inconclusive. As regards all of these factors our successors will be much 
more fortunately situated. 
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V. CORRELATION OF SEDIMENTARY AND CLIMATIC RECORDS 
By IsatAaH BOWMAN 


AMERICAN GEOGRAPHICAL SOCIETY OF NEw YORK 


Hitchcock as early as 1841 suggested that the clays of the Connecticut 
Valley showed annual banding. Certain regular alternations in Cretaceous 
strata in Colorado were thought by G. K. Gilbert (1895) to have a possible 
correlation with the precessional period of about 21,000 years; and by 
using this period as a time unit he calculated that 3900 feet of sediments 
(Benton, Niobrara and Pierre) were deposited in 20,000,000 years. Ber- 
key (1905) provisionally interpreted the double (winter and summer) 
laminations in a glacial clay bed 35 feet thick at Grantsburg, Wisconsin, 
as a measure of the years, 1700 in all, in which the deposit was laid down. 

It will be observed that Gilbert and Berkey made studies in a single 
location. They were seeking measures of the periods of time in which 
given deposits were laid down—measures akin to those of Michelson on the 
speed of light, in that the end in view was a more accurate yardstick. 
The interpretation of sediments has now passed into two more advanced 
stages: (1) the cross-dating of sediments laid down in separate localities 
under conditions so nearly comparable as to provide, not identical, but at 
least closely similar records, and (2) tentative correlations between such 
equivalent records and other records of similar or greater precision such as 
the tree ring analyses of Douglass and the solar radiation results of Abbot. 

De Geer’s work on the “‘varved’’ clays of Sweden and Antevs’ critical 
field studies in New England and eastern Canada are among the most 
noteworthy of the attempts to cross-date sediments. De Geer, looking 
for a general geochronological scale, boldly seeks to correlate the ‘‘curves’”’ 
of banded clays (varves) from Asia, South America and North America 
with the varve curves of Sweden, while Antevs challenges these results and 
insists on a more rigid theoretical basis and on much wider field work, 
namely, the checking of such correlations by the study of the record of 
contemporaneous or at least parallel events in moraines, eskers, striae 
and the like. A number of observers have suggested the possible correla- 
tion of the varying thickness of varved clay bands, or ups-and-downs in 
varve “curves,” with changes in solar radiation; and Reeds! has essayed 
such a comparison. The theory runs as follows: changes in the degree 
of solar radiation affect the melt rate of ice through variations of rainfall, 
temperature and cloudiness. The yearly layers of transported sediment— 
each composed of a thicker and lighter summer layer and a thinner and 
darker winter layer—are reasonably assumed to vary in thickness accord- 
ingly. 

The so-called varved clays are of glacial-marginal origin. The lakes in 
which they were deposited were provided with sediment-bearing water 
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that flowed from the front of the continental ice sheet. A great deal of 
rock flour—fine silt and clay—was included in the discharge. While 
the coarser detritus came to rest near shore the fine silt was widely dis- 
tributed over the lake floors before it came to rest and the clay remained in 
suspension for a much longer period and upon deposition formed the char- 
acteristic dark winter band. The strongly marked banding in clays of 
glacial-marginal lakes reflects the strong contrasts of rigid winter and 
warmth-compelled summer discharge on the front of a retreating ice 
border. 

The fresh-water lakes in question were short-lived. With the retreat 
of the ice that impounded them, outlets were provided at lower levels and, 
in many cases, the entire lake floor was laid bare. Within distances of 
a few miles up to 100 or 200 miles such lakes may have existed in the same 
period or portions of the same period and thus have had corresponding 
experiences which have been found to be clearly identifiable. The rate 
of retreat of the ice—about a mile in 22 years in central New England— 
the extent of the common experiences which given lakes enjoyed, and even 
the record of exceptional catastrophic events in the life of a given lake, 
such as a change of outlet, are among the results of varve studies. All 
such datings are relative. No one has yet found the bottom of the oldest 
deposit nor the top of the youngest. 

When we come to the correlation of the varved clays with tree rings there 
is an obvious gap between even the oldest New England trees and the date 
of extinction of the much more remote lakes in which there were formed 
seasonally banded clays. One naturally turns to the West, where trees 
grow that are more than 3000 years old. The most promising localities 
in the Great Basin have been searched for varves and a few sequences 
have been identified by Antevs?; but they are short and he concludes: 
“Correlation between rate of growth of the Big Tree and fluctuations of the 
lakes in the arid regions in the western States should be avoided, until 
the conditions governing the growth of the sequoia are better known, 
and the datings of the lake fluctuations are more accurate.” 

There are two prime difficulties in a tree-lake correlation in the Great 
Basin. First, the Great Basin lakes are and have been mostly salt, with 
short fresh-water interruptions; and the salt flocculates the sediment 
so that it is deposited promptly, leaving little colloidal material to settle 
upon the off-shore reaches of the lake floor. Second, there is a gap between 
the complicated history of the Great Basin lakes and the oldest trees. 
Two steps seem indicated. The late sedimentary history of the Great 
Basin requires much more detailed study by modern techniques before 
we have a satisfactory interpretation. The sediments of mountain border 
lakes also require equally detailed study for there if anywhere one should 
be able to find the bridge between the more distant lakes of the Great 
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Basin on the one hand and the mountain glaciers and Big Trees on the 
other. 

Fortunately, in the case of the lakes of the Great Basin there is a marked 
recurring advance and retreat of lake waters, an exaggerated expression of 
climatic changes. In recent years all of the lakes, with the exception of a 
few of the largest, have either disappeared or almost disappeared. Their 
floors are now available for sectioning for the first or possibly the second 
time since white occupation of the West. As we have already noted, the 
salt content of the lakes of interior basins produces rapid deposition of 
even the finer sediments in contrast to the slower differential settling 
of fresh-water lakes. Yet the salt content was and is variable and some of 
the lakes now salt were once fresh through overflow at times of greater 
expansion, presumably in the Pleistocene though also in post-glacial time 
as lake has drained into lake. Moreover, it is worthwhile to look for 
seasonal characteristics possibly in the salt deposits themselves in the 
manner suggested by Gale’ in explanation of the potash deposits of Alsace. 

During 1931, 54 samples of lake sediments were collected in 13 different 
basins of the Great Basin province, each sample about a foot long, two 
inches wide and an inch deep. At the same time about 80 tree ring sec- 
tions were obtained from selected localities nearby. 

A comparative study is in progress to determine: 

(a) Vertical changes in physical and chemical character. 

(b) The degree to which the annual deposits are recognizable and vary 
in thickness and if such variations are cyclic or have any recognizable re- 
lation to the rainfall and evaporation records. 

(c) The presence or absence of correspondences between the varia- 
tions of thickness or constitution of lake sediments and the variation of 
tree ring thicknesses in sections cut from trees growing in the same 
watersheds. 

There are at least three main lines of further investigation in the Great 
Basin that are held in view: 

(a) The detailed study of the sedimentary process in selected sites which 
will show types of sedimentary alternations that may prove recognizable 
in sections now exposed; 

(6) Detailed study of tree ring records on lake floors now exposed 
after long submergence; 

(c) Comparative study of floor deposits in a north-south direction 
with the object of discovering the regional differences in sedimentation 
habit possibly comparable to the regional differences in rainfall habit and 
evaporation rate. 

Both theory and fact point to the conclusion that lakes of dry regions 
may register the larger climatic changes even if the smaller changes go 
unrecorded or prove not to be recognizable. Such lakes seem to accent 
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the degree of climatic change though with a much more pronounced lag in 
relati-n to their water-supply fluctuations in contrast to the close relation 
exhibited by fresh-water lakes. Those salt-water lakes that have no out- 
let spread widely following periods of exceptionally heavy rainfall. Follow- 
ing periods of drought they may dry up altogether and their total salt con- 
tent is then deposited upon the floor of the basin. The recurrent and 
possibly periodic expansion and contraction of arid-basin lakes produces 
an alternation of salt and sediment that bears on its back, so to speak, the 
fainter alternations of the seasons and of groups of drier and wetter years. 
As muds settle they seal up the salts deposited in an earlier and drier phase 
of the climatic cycle. Only the comparative and intensive study of a 
widely extended series of lakes in a climatic province whose present-day 
or recent habits can be determined instrumentally would seem to provide 
an adequate observational basis for the testing of hypotheses already sug- 
gested by a number of observers. 

By taking the total complexity of climate and hydrography into account 
we can hope to piece out the regional story which lakes seem able to tell. 
If their history could be stated with some degree of refinement of detail 
it might furnish us with the missing factors that we need to know to fill 
the gap between the larger and the shorter climatic oscillations. Their 
sediments go far back of the earliest tree records, just as the tree records 
go far back of the period of instrumental weather observations. Equally 
important is the diversity of climatic experience, region by region, in the 
same broad province. R. J. Russell* has recently represented the three- 
fold division of the Great Basin in terms of “‘desert years’’ and Meinzer® 
has called attention to the clear distinction between the northern and 
southern sections of the lake region of the Great Basin in the Pleistocene 
as well as today. We can hardly suppose that correlation of lacustrine 
evidence is possible unless account be taken of regional climatic differences 
such as are demonstrable in the case of adjacent climatic provinces today; 
precipitation, for example, being in opposite phase often within relatively 
short distances. 

The first great deficiency in correlation studies of lake sediments is the 
absence of detailed measurements of recent floor deposits in places where 
rainfall and evaporation records and tree rings are also available. Trees 
upon slopes adjacent to lake floors have not been systematically studied in 
relation to lake history. Tree trunks are known to be exposed directly 
upon basin lake floors that are periodically dry (yellow pine in Granite 
Lake, Spokane County, Washington®). So far as I know no one has 
attempted the correlation of tree ring records and sediments in critical 
situations nor has there yet been made a direct comparison of contem- 
poraneous sedimentary and precipitation records. 

The most neglected field appears to be essentially direct observation 
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upon the sedimentary process. The first requirement is a technique of 
observation. W. A. Johnston’ has come close to it and has shown the 
reality and degree of seasonal banding in the waters of a given lake, of 
known or knowable rainfall and snowfall and measurable stream discharge, 
evaporation rate, cloudiness and the like. And R. W. Sayles® has dis- 
cussed criteria, including the important and neglected point that ‘‘the 
ratios between the thicknesses of the coarse:and fine components should 
be studied rather than the actual thicknesses of the individual components, 
for the thicknesses of the layers are due to other factors in addition to 
the climatic factor.”” I would suggest that an analysis of the thicknesses 
of the winter layers might also show systematic variations of possibly 
diagnostic value. Bottom sampling through lake water has, when taken 
alone, quite limited value for the purpose in view. Deeper sections are 
required and an absence of marginal distortion due to the thrust of the 
sampling machine. 

Lake sediments are, of course, in process of deposition today. Every 
type of lake-floor section has its counterpart in a living lake or in one but 
recently vanished. Moreover, we have the means available for directly 
observing the growth of trees in basins of sedimentation. We can di- 
rectly observe the weather elements in their climatic or long-range com- 
positions in localities of critical importance. Spotty “‘correlations’’ have 
high suggestive value but they do not give us laws, only miscellaneous 
correspondences. Certain laws of solar occurrences may be well estab- 
lished; the sedimentary effects of such occurrences are another story, a 
complicated story. If correspondences and resemblances are the goal, 
we have reached the goal. If the actual regional operation of the forces 
is the objective we have only opened the book. The lesson is still to be 
read. The status of experimental work is suggested by the fact that in 
his excellent paper of 1915, Gale® turned to the work of T. M. Chatard, 
published in 1890 (based on data gathered in 1886) for the composition of 
crystalline residues obtained by evaporating the mother liquor of a salt 
lake (Mono Lake, California). 

The apparent conclusiveness of the first varve correlations has been 
misleading. Varved clays furnish a geochronological yardstick that has 
to be used with caution. The old question reappears: What is a cor- 
relation? There is no general agreement on the point. A wide gap has 
yet to be closed between the mathematics of curve analysis and the realities 
of the depositional process. For example, a recent paper points to solar 
radiation variations ‘‘not different in kind’ from those of varved clays. 
We are not told what constitutes “difference in kind.’’ The author is 
surprised to see how closely solar fluctuations ‘‘simulate” those of varved 
clays. What is simulation and how close must it be to be significant? 
A tentative correlation is offered because of “‘close agreement in the form 
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of the graphs.’’ Similarities of form are suggestive. How can they be 
conclusive or become the groundwork of correlations unless independent 
lines of evidence give support to a theory based on mere resemblance? 

The moment we detach a set of varve ‘‘curves’’ from the regional en- 
vironment and attempt to match them with a set from a far distant 
locality we are bound to require a degree of correspondence in the curves 
that is unattainable. As Antevs!! has pointed out: “If warm summers 
in one region corresponded to cold summers in another region, thick varves 
in the former area should be compared with thin varves in the latter. 
And if warm summers in one region corresponded, now with warm, now 
with cold summers in another region (as is the case between North America 
and northern Europe at present), thick varves in the former area should be 
compared, now with thick, now with thin varves in the latter. In other 
words, correlations between graphs can be made only for areas that are 
known, or with reasonable certainty can be assumed, to have experienced 
about the same annual fluctuations in the total summer temperature.” 

The essential point of this criticism is that connections between varve 
curves are not rigid and absolute but depend upon a set of independent 
supporting conditions of glacial or aqueo-glacial deposition. If the lo- 
calities are not far apart and if the late-glacial history of the region is known 
in some detail, a presumption is raised or a conclusion definitely reached 
that the two localities had a similar climatic history during a part of the 
period of clay accumulation. Under these circumstances there should be 
possible a matching of the curves with a degree of correspondence well 
above the level of accident. Thus Antevs finds an equal number of curves 
from Brattleboro and the Hudson, 80 miles apart, in agreement to the 
extent of 84 per cent; Concord and Albany, 120 miles, 75 per cent; New- 
burg and Hartford, 70 miles, 79 percent; and so on. Upon rigid com- 
parison, curves from the two sides of the Atlantic show correspondences in 
the main from 50 to 60 per cent only and thus have no meaning because 
‘‘Varve graphs are individually only moderately distinctive, and, in addi- 
tion, being records of summer heat, they present periodic more or less 
similar fluctuations.’’!2 Under these circumstances a close analysis of 
the curves and a high degree of correspondence are a rigid requirement if 
we are to get away from the mere opinion of an individual. 

The fundamental difference between Antevs on the one hand and De 
Geer on the other with respect to correlations of varves as between New 
England and Scandinavia and Patagonia and Scandinavia is illustrated 
by the following quotation from Caldenius: ‘“Two teleconnected diagrams 
from different territories of glaciation can hardly be expected to be similar 
to the same extent as two diagrams connected together from the same terri- 
tory of glaciation. That would be to set expectations to too high a pitch, 
when we know all the. matters which have a disturbing influence upon the 
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copying of the variations of the summer temperature in the thickness of 
the varves. The reflecting similarities of the teleconnecting diagrams 
become for that reason fewer in number, and they are confined to certain 
groups of varves, but the teleconnections are based upon the equal length 
of these repeated groups of varves with their different characteristic simi- 
larities and on the equal length of groups of varves occurring between them, 
where conformity is missing or where the curve, on which the teleconnec- 
tions are based, is indifferently formed. It is obvious that the conformity, 
stated in percentage indicates the similarities which are intended and are 
at hand, in an awkward and often misleading manner.’’* 

This carries the comparison far outside the requirements of rigid correla- 
tion and into the field of simple comparison of form—and not the details 
of form as in the case of comparisons between varves from a given lake and 
its neighbor, but to comparisons of ‘‘certain groups of varves,”’ that is, 
constellations of forms that simply look alike. This is suggestion, not proof; 
and it is valuable only if kept in the field of suggestion and not combined 
with sun-spot-varve correlations that have their own modes of error. 

With every refinement of technique in varve studies there is an advance 
away from the position that varves are self-contained entities and toward 
the position that varves are parts of a record of which the whole is glacial 
history in all its complexity. One of the most puzzling features is the 
present-day regional variation in climatic habit. Climatic provinces— 
even those that are closely adjacent—do not have a uniform history. 
Their cycles are not commonly in the same phase if we take groups of 
adjacent regions into account though their phases may in time reveal 
causal relationships. We can only suppose that these differences. have 
always existed. The consequence is that a single idea or a single technique 
of study cannot be applied wholesale. The problem requires to be broken 
up into its regional components. 

The studies of Clayton'‘ have shown what some of these regional com- 
ponents are. His most comprehensive generalization is that, within the 
tropics, the temperature increases and decreases with the solar radiation, 
but with a slight lag, while from 30° to 60° N. or S. there is an inversion. 
A correlative effect of increase of solar radiation is to intensify the belts 
of low pressure near the equator and near the 60th parallel and displace 
the intervening high-pressure belts with the consequence that there is an 
increase in the intensity of atmospheric circulation. 

Clayton then goes on to show how this broad terrestrial generalization 
is affected by the well-known changes of pressure over continents and 
oceans with the progression of the seasons. He also carries his analysis 
forward to the point of including both the Arctic and Antarctic regions 
in the scheme of concurrent change in atmospheric pressure and solar 
radiation. He recognizes, however, three variables which influence ‘‘the 
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pressure and other weather conditions at any one place on the earth sur- 
face.” First is the variation of pressure intensity; second, an annual 
period in which pressure conditions are reversed between summer and 
winter, as over continents and oceans of temperature latitudes; and third, 
a movement of the centers of action north and south of their mean position 
with the varying intensities of solar radiation. He concludes, ‘‘As a re- 
sult of these three variations the weather at very few places on the earth 
shows a simple direct influence of solar variation.’”’ By this he means 
that ‘‘even when the annual period is eliminated by taking the same month 
in succeeding years, the two remaining variables cause the weather at any 
given place to follow the solar radiation, sometimes directly and sometimes 
inversely.’’ These complexities call for elaborate interpretations from 
the observations at any one station in terms of or in relation to stations 
elsewhere in the world. He believes that up to the present time the 
studies have not been conclusive and do not show the exact manner 
in which atmospheric waves are formed in response to changes in the 
atmosphere brought about by solar influence. He seems to occupy 
a neutral ground with respect to the changes of longer range. While 
believing that the day by day weather changes are directly related to solar 
radiation and while pointing to the close analogy between such changes 
and those of longer period, he concludes that the long-range changes have 
not yet been shown to be due to solar changes except possibly in part.% 

In the same way that the lakes of enclosed basins exaggerate climatic 
changes by wide expansions and contractions, so streams exhibit cycles 
of run-off much more clearly than do rainfall records. This is because 
the run-off is cumulative and increases three times as fast as the rainfall. 
Correlations seem to be established between curves of stream flow, the sun- 
spot cycle, twice that cycle, and the Briickner cycle, which is three times 
the sunspot cycle. The matter of stream flow has been studied in relation 
to the cyclical changes in lake levels, particularly the Great Lakes." 
Thus the cumulative as well as the amplifying nature of stream flow may 
prove to be combinable with the records of lake-level changes and provide 
the sedimentation cycle of lake floors with an accented rhythm of great im- 
portance in determining the reality of solar effects. I should like to 
underscore this point because our attention has been focused for so long 
upon a special kind of lake deposit, namely, the varved clays. These 
were formed in glacial-marginal lakes under conditions of seasonal changes 
in water discharge from melting ice. They represent local drainages and 
special conditions of drainage that have been assumed to be exceptionally 
favorable to a type of lamination that records the seasons. It remains to 
be seen from the direct study of the sedimentation process in existing lakes 
of various types whether periodicities of greater amplitude may not be 
determinable, as for example upon lakes of high latitude or altitude whose 
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tributary streams have strong seasonal differences of discharge and some 
of which are affected by ice melting and others not. 

Three main lines of inquiry suggest themselves. The first is the study 
of the sedimentation process in lakes lying in different climatic zones; 
second, the study of sediments now exposed in widely contrasted climatic 
zones; and third, the study of changes in river discharge and lake levels 
in widely different climatic zones. On theoretical grounds one would 
conclude that the relationships between rainfall, evaporation rate and 
sun-spots, on the one hand, and changes in river and lake levels on the 
other hand, if strongly suggested or established, would vary in type and in 
degree of correlation as between the tropics and temperate zones and as 
between the humid and arid regions. Biel’s recent map of rainfall varia- 
bility gives us a comprehensive world view of the fact long known in 
meteorology that, in general, the drier the climate the greater the rainfall 
variability from year to year. Each type of climatic province has its 
special regimen of rainfall, stream discharge, change in lake level, dif- 
ference in evaporation rate and in character of the sedimentary process. 
The result will necessarily be a marked difference in the clearness of register 
of these facts in the sedimentary record. It is equally certain that some 
types of records will necessarily show more clearly than others those re- 
lationships that are thought to be real and that as yet are expressed in 
correlation coefficients too low in value to be conclusive. To take a 
specific instance, Lake Victoria!’ in Africa, where the rainfall seems pri- 
marily responsible for variations in level of the lake, the correlation co- 
efficient of 0.915 has been found between the rainfall of the plateau region 
and the level of the lake, a correlation coefficient of 0.64 being found be- 
tween sun-spots and rainfall in the same region. Similarly, for Lake 
Nyasa, Dixey'’® has pointed out the cyclical nature of the variations in 
lake level and their correspondence with sun-spot numbers, the lake rising 
at sun-spot maxima. A much lower degree of correspondence has been ob- 
served between sun-spots and the successive levels of Lake George, Aus- 
tralia. An 11-year periodicity seems to be recognizable and a less defi- 
nite 19-year periodicity. C. E. P. Brooks in commenting upon these low 
correlation values observes that ‘‘weather cycles are treacherous things, 
and it would not be safe to base a forecaston them.’’ His conservatism is 
based upon the short-range changes in Lake George. Other influences 
than rainfall evidently have a predominating effect and the lake may prove 
to be one of those having a regimen not susceptible of analysis in terms that 
will yield results for any of the correlations here held in view. 

The phenomena that we examine in the case of climatic change come to 
us in diluted form. The saying is that we have to take the sunlight ‘‘as 
it comes to us,” that is, through the atmosphere, just as we have to take 
the sedimentary record as it stands. And as it stands it is fragmentary, 
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besides being the expression of a number of forces some of which are 
known to be without regularity of occurrence or with a regularity of such 
long periodicity as to have no meaning in the short-term studies with which 
we are here mainly concerned. Sediments, like some tree rings, may 
express the complex result of a number of agencies. In an attempt to 
correlate sediments, climatic changes, water supply, tree rings or any 
other set of elements, we are confronted with the difficulty of dilute effects. 
Advance in correlations that involve sediments seems to depend upon a 
new technique of study of the sedimentation process and the search for 
lakes, whether existing or extinct, whose bottom deposits have been in 
delicate relation to specific climatic conditions. The lakes of the Great 
Basin are among those which seem to promise data of critical value. 
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In response to the request of the Chairman for discussion of the above 
contributed papers, the following remarks were made from the floor by 
Messrs. A. E. Kennelly, E. W. Brown, H. N. Russell and E. B. Wilson. 


We seem to be justified in concluding from the interesting papers, just read, that 
correlations are found between variations in certain solar phenomena (sun-spots and 
radiation) and variations in certain terrestrial phenomena (magnetic, thermal and 
vegetative). These correlations are not sufficiently clear to warrant definite conclusions 
at present; nevertheless, they seem to be sufficiently evident to make more extended 
studies well worthwhile, especially when we consider how recent are the observations 
on which these correlations have been obtained. Quantitative sun-spot records are 
not much more than one century old, and quantitative radiation records only a few 
decades. These intervals are but droplets in the stream of solar and terrestrial time. 

It would be most interesting to unravel, if possible, the history and causes of the sun- 
spot cycles. How far back can they be traced in vegetation? Are they discernible in 
the sections of fossil trees? Are these spot cycles caused by actions below the solar sur- 
face, or in external space? These cosmic phenomena are surely eminently worthy of 
study by the National Academy of Sciences. Even if we assume that their complexity 
will baffle all attempts at future long-range weather forecasting on the earth, they may 
be expected to lead to knowledge of the past history of the earth and sun. We are, 
therefore, justified in believing that further research will establish historical connections 
between the sciences of geology, heliology and meteorology. 

A. E. KENNELLY 


The analysis of a continuous series of observations into cycles, usually called harmonic 
analysis, is useful: 

(a) When physical theory indicates the presence of a cycle and we wish to know its 
amplitude and phase or to correct its period; 

(b) When it is desired to discover a cycle either for purposes of prediction or to give 
an indication of physical theory. 

Experience has, however, shown that the discovery of a cycle which was not obvious 
from the graph has little value. A hidden cycle means that it is so overlaid with other 
effects that it cannot be used for prediction or for the discovery of the principal physical 
causes of the oscillating phenomena. 

Analysis into many hidden cycles has usually no physical meaning. Mathematical 
theory shows that we can always represent a graph by cycles. Since each harmonic 
term introduces three constants, the more cycles we take the better should be the 
representation. 

A partial test of the reality of a cycle is furnished by calculation of the probable errors 
of its constants. This should always be done whether the cycle is hidden or not. 

The question of correlation between two sets of oscillating phenomena is to some ex- 
tent mixed up with that of harmonic analysis. The mere fact that the two sets appear 
to give a cycle which has nearly the same period in each has usually no physical meaning 
outside of a known physical cause. Where several cycles are found, the laws of chance 
alone may give us one in one set which nearly corresponds with one in the other set. 
This is particularly the case with the sub-harmonics. If we find a cycle for example, 
of twelve years, we are almost sure to find one. or more of six, four, three, etc., years. 
If we find only an occasional multiple, the cycle should be regarded as highly suspicious 
as to its reality. 
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The question as to the degree of correlation between two sets of phenomena needed 
to assert the existence of a common physical connection, is to some extent an open one. 
Usually, however, unless a correlation is obvious it is scarcely worth detailed study for 
the same reasons as those given for the existence of a cycle. 

I am in full agreement with Humphreys when he states that there are only two known 
climatic cycles, the day and the year and these result from geometrical changes of posi- 
tion of the earth’s surface with regard to the sun. The nature of the motion of a fluid 
like the atmosphere is such that irregularity is almost the necessary consequence—that 
is, that it is impossible to predict its motions by mathematical cycles alone. The 
more promising line of attack is that which has already been successful, namely, corre- 
lation of phenomena which follow one another in time though not necessarily in place. 
The physical basis for this is the fact that the motions of the atmosphere and their effects 
depend to a greater extent on phenomena which have immediately preceded them than 
on those which are separated from them by a longer interval of time. 

‘Trigger’? action is sometimes invoked. Such action demands a peculiar set of cir- 
cumstances which will only rarely occur and it involves the existence of instability. 
When instability is present predictions as to the time when future events will take place 
are rarely possible. 

E. W. Brown 


The conspicuous periodicities which appear in solar phenomena (and the more doubt- 
ful ones which are found in the weather) differ radically from those familiar to the 
astronomer. The latter depend, in one way or another, on simple dynamical motions 
(such as the Earth’s rotation and orbital revolution) and are precise timekeepers. A 
similar oscillation in a body of the Sun’s mass and density should have a period of two 
hours or less. The sun-spot cycle is more than fifty thousand times longer, and is 
evidently based on some fluctuation of thermal or hydrodynamic character. Sun-spots 
are periodic in space (heliographic latitude) as well as in time; and this, and the re- 
versal of polarities in successive cycles, shows clearly that each cycle is a unit—an out- 
burst which affords a release for some sort of internal strain, each beginning before the 
last died out. One would not expect such cycles to be more than very roughly periodic, 
like the eruptions of one of the more regular geysers. 

This makes an analytical treatment of the periodicity very difficult, for almost all 
the available methods assume the existence of fixed underlying periods. The formulae 
derived from such analyses habitually fail to predict the future course of the phenomena, 
not on account of any errors in the reckoning, but because the processes themselves are 
not accurately predictable. Michelson, whose judgment in physical matters was 
excellent, after studying the sun-spot curve, concluded that it should be represented 
by a function of variable period, amplitude and phase.! 

In the case of the weather, we have conspicuous, but not exactly repeating, cycles of 
diurnal and annual period, weak, but probably real, evidences of the main solar cycle, 
and evidence also of other cycles, which may be determined by terrestrial conditions, 
for example, in the atmosphere or the ocean. These, like the solar cycle, are probably 
inherently irregular in period and amplitude. Superposed on them are the much larger 
variations connected, for example, with cyclonic storms. 

To disentangle the various cycles from one another and from the other fluctuations, 
is a task of excessive difficulty, and I fear that even generalized prediction of weather 
conditions by the extrapolation of empirical cycles derived from observation, if attain- 
able, is far in the future. 

This opinion, though unfavorable to certain methods of prediction, is very far from 


1 Astrophys. Jour., 28, 273 (1918). 
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being adverse to the continuance of such admirable investigations as Dr. Abbot and Dr. 
Douglass have described. . Each of these researches has already far more than justified 
its existence and its continuance. The one has provided a precise chronology of funda- 
mental value to the western archaeologist; the other has proved the existence of fluctua- 
tions in the sun’s radiation, enduring for months or years, whose explanation offers im- 
portant problems which the astrophysicist is as yet hardly ready even to attack. Fur- 
ther investigation along these and similar lines is sure to be fruitful and may lead to 
important practical results. 
H. N. RUSSELL 


Previous speakers have referred to empirical statistical analysis on the one hand and 
to analytical mathematical analysis on the other. I wish to mention a point which is 
of equal importance to both. I have worked very little with climatic cycles. The only 
work of that sort which I have done is in attempting to find some reasonably high cor- 
relations between the death rate from penumonia and various aspects of the weather such 
as temperature, the diurnal variation of temperature, or the suddenness of temperature 
changes, or humidity or barometric pressure, and I have found that the correlations are 
low. Even if one treats the problem by partial correlation one has an outstanding varia- 
tion which is large. 

Important as such a treatment is whether for the empirical descriptive representation 
of a part of the phenomenon or as a suggestion for an analytical and mathematical 
discussion of the rational basis of the phenomenon, it is scientifically of the utmost 
importance not to give up the exploration for other possible variables which may be much 
more highly correlated with the phenomenon. It is quite possible that in respect to 
the mortality from respiratory diseases or the growth of trees some variable not now 
tabulated in the meteorological records is of more crucial importance than any of the 
variables which are tabulated. For example, there is the ionization of the air, which is 


decidedly variable, on which certain preliminary experiments and calculations (which 
are not ready for release because the series is still too short to be dependable) indicate 
the possibility that it may be even more highly correlated with the respiratory mortality 
than any other variable yet tried. The search for new and better variables is fundamen- 
tal to both the empirical and the analytical aspects of science. 


E. B. WILson 











